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Superconducting properties of metallic nanowires can be entirely different from those of bulk 
superconductors because of the dominating role played by thermal and quantum fluctuations of 
the order parameter. For superconducting wires with diameters below ~ 50 nm quantum phase 
slippage is an important process which can yield a non-vanishing wire resistance down to very low 
temperatures. Further decrease of the wire diameter, for typical material parameters down to ~ 10 
nm, results in proliferation of quantum phase slips causing a sharp crossover from superconducting 
to normal behavior even at T = 0. A number of interesting phenomena associated both with 
quantum phase slips and with the parity effect occur in superconducting nanorings. We review 
recent theoretical and experimental activities in the field and demonstrate dramatic progress in 
understanding of the phenomenon of superconductivity in quasi-one-dimensional nanostructures. 
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I. INTRODUCTION 



The phenomenon of superconductivity was discovered 
[l[ as a sudden drop of resistance to immeasurably small 
value. With the development of the topic it was realized 
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that the superconducting phase transition is frequently 
not at all " sudden" and the measured dependence of the 
sample resistance R{T) in the vicinity of the critical tem- 
perature T(7 may have a finite width. One possible rea- 
son for this behavior - and frequently the dominating 
factor - is the sample inhomogeneity, i.e. the sample 
might simply consist of regions with different local crit- 
ical temperatures. However, with improving fabrication 
technologies it became clear that even for highly homoge- 
neous samples the superconducting phase transition may 
remain broadened. This effect is usually very small in 
bulk samples and becomes more pronounced in systems 
with reduced dimensions. A fundamental physical reason 
behind such smearing of the transition is superconducting 
fluctuations. 

An important role of fluctuations in reduced dimen- 
sion is well known. Above Tc such fluctuationsyield an 
enhanced conductivity of metallic systems 0, 0, Q ■ 
For instance, the so-called Aslamazov-Larkin fluctuation 
correction to conductivity SaAL ^ {T — Tc)^*-^"^/^' be- 
comes large in the vicinity of Tc and this effect increases 
with decreasing dimensionality D. Below Tc - according 
to the general theorem [1, 0| - fluctuations should de- 
stroy the long-range order in low dimensional supercon- 
ductors. Thus, it could naively be concluded that low 
dimensional conductors cannot exhibit superconducting 
properties because of strong phase fluctuation effects. 

This conclusion, however, turns out to be some- 
what premature. For instance, 2D structures undergo 
Berezinskii-Kosterlitz-Thouless (BKT) phase transition 
[1, H, as a result of which the decay of correlations 
in space changes from exponential at high enough T to 
power law at low temperatures. This result implies that 
at low T long range phase coherence essentially survives 
in samples of a finite size and, hence, 2D films can well 
exhibit superconducting properties. 

Can superconductivity survive also in (quasi)-lD sys- 
tems or do fluctuations suppress phase coherence thus 
disrupting any supercurrent? The answer to this ques- 
tion would clearly be of both fundamental interest and 
practical importance. On one hand, investigations of this 
subject definitely help to encover novel physics and shed 
more light on the crucial role of superconducting fluctu- 
ations in ID wires. On the other hand, rapidly progress- 
ing miniaturization of nanodevices opens new horizons 
for applications of superconducting nanocircuits and re- 
quires better understanding of fundamental limitations 
for the phenomenon of superconductivity in reduced di- 
mension. A detailed review of the present status of this 
field is the main purpose of this paper. 

It was first pointed out by Little [11] that quasi-one- 
dimensional wires made of a superconducting material 
can acquire a finite resistance below Tc of a bulk material 
due to the mechanism of thermally activated phase slips 
(TAPS). Within the Ginzburg-Landau theory one can 
describe a superconducting wire by means of a complex 
order parameter ^'(a;) = I'i' {x)\e^'^^^\ Thermal fiuctua- 
tions cause deviations of both the modulus and the phase 
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FIG. 1: Scliematics of the phase slip process. Spatial variation 
of the amplitude of the order parameter l^*] (left axis, dashed 
line) and phase ip (right axis solid line) at various moments of 
time: (a) before, (b) during and (c) after the phase slippage. 



of this order parameter from their equilibrium values. A 
non-trivial fluctuation corresponds to temporal suppres- 
sion of |5'(a;)| down to zero in some point (e.g., a; = 0) 
inside the wire, see Fig. 1. As soon as the modulus of 
the order parameter |4'(0)| vanishes the phase ip{0) be- 
comes unrestricted and can jump by the value 27rn, where 
n is any integer number. After this process the modu- 
lus |vE'(0)| gets restored, the phase becomes single valued 
again and the system returns to its initial state accumu- 
lating the net phase shift 2iTn. Provided such phase slip 
events are sufficiently rare, one can restrict n by n = ±1 
and totally disregard fluctuations with |n| > 2. 

According to the Josephson relation V = hip/2e each 
such phase slip event causes a nonzero voltage drop V 
across the wire. In the absence of any bias current the net 
average numbers of "positive" (n = -1-1) and "negative" 
(n = —1) phase slips are equal, thus the net voltage drop 
remains zero. Applying the current / oc| 4* ^ Vip one 
creates nonzero phase gradient along the wire and makes 
"positive" phase slips more likely than "negative" ones. 
Hence, the net voltage drop V due to TAPS differs from 
zero, i.e. thermal fluctuations cause non-zero resistance 
i? = V/I oi superconducting wires even below Tc- We 
would also like to emphasize that, in contrast to the so- 
called phase slip centers [H, [H, [13, [lB| produced by a 
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large current above the critical one I > Ic, here we are 
dealing with fluctuation-induced phase slips which can 
occur at arbitrarily small values /. 

A quantitative theory of the TAPS phenomenon was 
first proposed by Langer and Ambegaokar (16j and then 
completed by McCumber and Halperin jl7| • This LAMH 
theory predicts that the TAPS creation rate and, hence, 
resistance of a superconducting wire R below Tc are de- 
termined by the activation exponent 

R{T)^exp{-U/T), U^^^sar), (1) 

where U{T) is the effective potential barrier for TAPS 
proportional to the superconducting condensation energy 
{No is the metallic density of states at the Fermi energy 
and Ao(T) is the BCS order parameter) for a part of the 
wire of a volume where superconductivity is destroyed 
by thermal fluctuations (s is the wire cross section and 
5(T) is the superconducting coherence length). At tem- 
peratures very close to Tc eq. ([T]) yields appreciable resis- 
tivity which was indeed detected experimentally [H, [l^ . 
Close to Tc the experimental results fully confirm the 
activation behavior of R{T) expected from eq. ((1]). How- 
ever, as the temperature is lowered further below Tc the 
number of TAPS inside the wire decreases exponentially 
and no measurable wire resistance is predicted by the 
LAMH theory [l^, [l3| except in the immediate vicinity 
of the critical temperature. 

Experiments O [l^ were done on small whiskers of 
typical diameters ~ 0.5 /im. Recent progress in nano- 
lithographic technique allowed to fabricate samples with 
much smaller diameters down to - and even below - 10 
nm. In such systems one can consider a possibility for 
phase slips to occur not only due to thermal, but also 
due to quantum fluctuations of the superconducting or- 
der parameter. The physical picture of quantum phase 
slippage is qualitatively similar to that of TAPS (see Fig. 
1) except the order parameter |^'(a;)| gets virtually sup- 
pressed due the process of quantum tunneling. 

Following the standard quantum mechanical argu- 
ments one can expect that the probability of such tun- 
neling process should be controlled by the exponent 
~ exp(— [//fiwo), i-e. instead of temperature in the ac- 
tivation exponent ^ one should just substitute Hujo, 
where uq is an effective attempt frequency. This is be- 
cause the order parameter field ^'(x) now tunnels under 
the barrier U rather than overcomes it by thermal ac- 
tivation. Since such tunneling process should obviously 
persist down to T = one arrives at a fundamentally 
important conclusion that in nanowires superconductiv- 
ity can be destroyed by quantum fluctuations at any tem- 
perature including T = 0. Accordingly, such nanowires 
should demonstrate a non-vanishing resistivity down to 
zero temperature. Assuming that hioo ~ Aq(T) one 
would expect that at Ao(r) < T < Tc the TAPS de- 
pendence ([1]) applies while at lower T < Aq (T) quantum 
phase slips (QPS) take over, eventually leading to satu- 
ration of the temperature dependence R{T) to a non-zero 



value in the limit T ^ Aq. 

This behavior was indeed observed: Giordano (20| per- 
formed experiments which clearly demonstrated a no- 
table resistivity of ultra-thin superconducting wires far 
below Tc- These observations could not be adequately 
interpreted within the TAPS theory and were attributed 
to QPS. Later other groups also reported noticeable de- 
viations from the LAMH theory in thin (quasi-) ID wires. 
These experiments will be discussed in Chapter 6. 

It should be noted, however, that despite these de- 
velopments the idea that in realistic samples supercon- 
ductivity can be destroyed by quantum fluctuations was 
initially received with a large portion of scepticism. On 
one hand, this was due to a number of unsuccessful at- 
tempts to experimentally observe the QPS phenomenon. 
On the other hand, some early theoretical efforts have 
led to the results strongly underestimating the actual 
QPS rate. Also, unambiguous interpretation of the ob- 
servations j20|] in terms of QPS was questioned because 
of possible granularity of the samples used in those ex- 
periments. If that was indeed the case, QPS could eas- 
ily be created inside weak links connecting neighboring 
grains. Also in this case superconducting fluctuations 
play a very important role |21[ however - in contrast to 
the case of uniform wires - the superconducting order 
parameter needs not to be destroyed during a QPS event. 

First attempts to theoretically analyze the QPS ef- 
fects ^2, 23, 24] - as well as a number of later studies 
- were based on the so-called time-dependent Ginzburg- 
Landau (TDGL) equations. Unfortunately the TDGL 
approach is by far insufficient for the problem in question 
for a number of reasons: (i) A trivial reason is that the 
Ginzburg-Landau (GL) expansion applies only at tem- 
peratures close to Tc whereas in order to describe QPS 
one usually needs to go to lower temperatures down to 
T -> 0. (u) More importantly, also at Tc - T < Tc 
TDGL equation remains applicable only in a special limit 
of gapless superconductors, while it fails in a general sit- 
uation considered here, (iii) TDGL approach does not 
account for dissipation effects due to quasiparticles in- 
side the QPS core (in certain cases also outside this 
core) which are expected to reduce the probability of 
QPS events similarly to the standa rd p roblem of quan- 
tum tunneling with dissipation [H, Hall^]. (iv) TDGL 
approach is not fully adequate to properly describe exci- 
tation of electromagnetic modes around the wire during 
a QPS event (this effect turns out to be particularly im- 
portant for sufficiently long wires). Thus, TDGL-based 
description of QPS effects simply cannot be trusted, and 
a much more elaborate theory is highly desirable in this 
situation. 

A microscopic theory of QPS processes in supercon- 
ducting nanowires was developed .29, .30J with the aid 
of the imaginary time effective action technique [s^, |3l[ . 
This theory remains applicable down to T = and prop- 
erly accounts for non-equilibrium, dissipative and electro- 
magnetic effects during a QPS event. One of the main 
conclusions of this theory is that in sufficiently dirty su- 
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perconducting nanowires with diameters in the 10 nm 
range QPS probability can already be large enough to 
yield experimentally observable phenomena. Also, fur- 
ther interesting effects including quantum phase transi- 
tions caused by interactions between quantum phase slips 
were predicted [ISIli- 

An important parameter of this theory is the QPS fu- 
gacity 

y ^ "Score 6Xp( »S'coi o)5 "Score ^ 5? J 

where is the dimensionless conductance of the wire 
segment of length ^. Provided is very large, typically 
9i ^ 100, the fugacity y remains vanishingly small, QPS 
events are very rare and in many cases can be totally 
neglected. In such case the standard BCS mean field 
description should apply and a finite (though possibly 
sufficiently long) wire remains essentially superconduct- 
ing outside an immediate vicinity of Tc- For smaller 
9s, ^ 10-;- 20 QPS effects already become important down 
to r — > 0. Finally, at even smaller ~ 1 strong fluctu- 
ations should wipe out superconductivity everywhere in 
the wire. We also point out that in the case of nanowires 
considered here the parameter is related to the well 

known Ginzburg number as GiiD ^ ^Is'^'^ ^ i-e. the con- 
dition ~ 1 also implies that the fluctuation region 
becomes of order Tq- 

Another important parameter is the ratio between 
the ("superconducting") quantum resistance unit Rq ~ 
Trh/2e'^ = 6.453 kf2 and the wire impedance Z^, — 

/i = Rq/2Zw, 

where C is the wire capacitance per unit length, L = 
47rA|_/s is the wire kinetic inductance and is the Lon- 
don penetration depth. Provided this parameter becomes 
of order one, /i ^ 1, superconductivity in sufficiently long 
wires gets fully suppressed due to intensive fluctuations 
of the phase of the superconducting order parameter. 
We note that both and /i scale with the wire cross 
section s respectively as oc s and ^ oc ^/s. It fol- 
lows immediately that with decreasing the cross section 
below a certain value the wire inevitably looses intrinsic 
superconducting properties due to strong fluctuation ef- 
fects. For generic parameters both conditions g^ ^ 1 -f- 10 
and fj, ^ I are typically met for wire diameters in the 
range y/s < 10 nm. 

A number of recent experimental observations are 
clearly consistent with the above theoretical conclusions. 
Perhaps the first unambiguous evidence for QPS effects in 
quasi-lD wires was reported by Bezryadin, Lau and Tin- 
kham [32] who fabricated sufficiently uniform supercon- 
ducting wires with thicknesses down to 3 5 nm and ob- 
served that several samples showed no signs of supercon- 
ductivity even at temperatures well below the bulk criti- 
cal temperature. Those results were later confirmed and 
substantially extended by different experimental groups. 



At present there exists an overwhelming experimental 
evidence for QPS effects in superconducting nanowires 
fabricated to be sufficiently uniform and homogeneous. 
Below we will analyze the main experimental results and 
compare them with theoretical predictions. 

Yet another interesting issue is related to persistent 
currents (PC) in superconducting nanorings. It was 
demonstrated [1^ that QPS effects can significantly mod- 
ify PC in such systems and even lead to exponential sup- 
pression of supercurrent for sufficiently large ring perime- 
ters. Another important factor that can substantially af- 
fect PC in isolated superconducting nanorings at low T 
is the electron parity number. Of particular interest is 
the behavior of rings with odd number of electrons which 
can develop spontaneous supercurrent in the ground state 
without any externally applied magnetic fiux [34] . 

The structure of our Review is as follows. A theory 
of superconducting fiuctuations in nanowires will be ad- 
dressed in Chapters 2-5. In Chapter 2 we discuss a gen- 
eral derivation of the real time effective action of a su- 
perconductor suitable for further investigations of ffuc- 
tuation effects at temperatures below Tc- We also for- 
mulate the Langevin equations and analyze their relation 
to TDGL-type of equations frequently used in the litera- 
ture. In Chapter 3 we adopt our general formalism to the 
case of superconducting nanowires and demonstrate the 
importance of superconducting fluctuations in such struc- 
tures. In Chapter 4 we will briefly review LAMH theory 
of thermally activated phase slips. Quantum phase slip 
effects will be analyzed in details in Chapter 5. Chapter 6 
is devoted to an elaborate discussion of key experiments 
in the field and their interpretation in terms of existing 
theories. In Chapter 7 we will analyze persistent cur- 
rents in superconducting nanorings. Chapter 8 contains 
a brief summary of our main observations and conclu- 
sions. Some technical details are presented in Appendix. 

II. EFFECTIVE ACTION, LANGEVIN AND 
GINZBURG-LANDAU EQUATIONS 

A. General formulation 

The starting point of our analysis is the formal expres- 
sion for the quantum evolution operator on the Keldysh 
contour or the so-called Keldysh "partition function" . As 
usually 21], the kernel of the evolution operator J' can 
be expressed in the form of a path integral over quantum 
fields defined both at the forward (below denoted by the 
subscript f) and the backward (denoted by the subscript 
b) branches of the Keldysh contour. We have 

X giS[4>pAF,VF,AF\-iS[tt>BAByB,AB\ ^2) 

Here i^p/B^i^F/B electron Grasmann fields, a =t, i is 
the spin index, A^/s are superconducting complex order 
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parameter fields (which emerge as a result of the standard 
Hubbard-Stratonovich decoupling of the BCS coupling 
term [lH), Vp/B and Ap/B are respectively the scalar 
and vector potentials. The action S is defined as follows 



iS = 



U{r) + eVjTp°' - AiP^^ 
No\A\^ 



Svr 



A 



(3) 



Here ef is the Fermi energy, U{r) is the disorder poten- 
tial, A'o is the density of states per unit spin at the Fermi 
level, and A is the BCS coupling constant. In our nota- 
tions the electron charge is e* = — e, i.e. we define e > 0. 
Here and below we employ the notation d'^X = dtd'^r. 

As the above action is quadratic in the electron fields 
one can integrate them out exactly. After that one arrives 
at the path integral 



J = / VAf/b VVf/b VAp/B e 



iSctt[A 



F/ B 1 ^ F/ E 



'■F/l 



where the effective action reads 



d^X 



A 



Tr In Q 



(4) 



(5) 



The inverse 4x4 matrix Green-Keldysh function Q ^ can 
be split into 2x2 sub-blocks (indicated by a hat): 



'22 



(6) 



where 

Qii = Cfz 

O22 = 

Here _ffo = " 
Hamiltonian. 



2{V,i)s} 



2{V,*s} 



raVc. 



mvt 



(7) 



-VV2m + U{r) 



^ J_( Vipp + ^Ap 
^ 2m\ V(^B + f A 



ep is the single electron 
(8) 



2e 









Vb 



2e 



(9) 



define the gauge invariant combiantions of the supercon- 
ducting order parameter phase and the electromagnetic 
potentials. The matrix A is constructed analogously. It 
reads 



A = 



Af 

As 



(10) 



Note that both fields A^ and A^ are now real since we 
have already decoupled the phase factors <Pf/b by the 
gauge transformation. 



B. Perturbation theory 

Although the above expression for the effective action 
^cff ® is exact it remains too complicated for practical 
calculations. In order to proceed let us - analo gous ly 
to the derivation in the Matsubara technique (sol. ISlj - 
restrict our analysis to quadratic fiuctuations. For this 
purpose we split the order parameter into the BCS mean 
field term Aq and the fluctuating part. Performing a shift 
^F/B ~* Aq + ^F/B "we re-define the order parameter 
field in a way to describe fluctuations by the fields A^/b. 
Expanding the effective action ([5]) in powers of A, vs and 
$ to the second order we obtain the action describing the 
quadratic fluctuations in the system 



iS'^^g = iSi + iSa + iScm + iS^A, (11) 



where 



iSi — Tr 



a-^e$(G'- G) - a^A{F + F) 



(12) 



defines the first order (in 'I>,'Us,A) contribution to the 
action, 



No 



J d^X{\Ao + Ap\^-\Ao + AB\') 
i Tr {(F<T,A)2 + {F&^Af + 2Ga,AGa,A} , 



(13) 



accounts for fluctuations of the absolute value of the or- 
der parameter field, 



i d^X 



ip2 Ei2 tr2 I it2 



Stt 



- -Tr^m*^(G-G) + 



{ -{V,'Os} + e$ ) G 



1 2 



<^z { ^{V, vs} + e$ j Fa, f ^{V, vs} - e<S> ] F 



^z ( -{V,t}s}-e$)G 



(14) 



describes electromagnetic fields and their coupling to the 
phase of the order parameter field and 



z5^A = Tr ^ a, ( -{V,-i)s} + e$ ) i^a.AG 



^z ( -{V,vs} + e.^\Gu,AF 
^z ( '-{V,Vs}~e<i>]Ga,AF 
'^z ( '-{^,vs}~e<^]Fa,AG\. 



(15) 
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is responsible for coupling of electromagnetic and phase 
fluctuations to the absolute value of the order parameter 
field. 

All the above contributions to the action are ex- 
pressed in terms of equilibrium non-perturbed normal 
and anomalous 2x2 Green-Keldysh matrices which enter 

as sub-blocks G, F, G, F into the 4x4 matrix 



^0 = 




(16) 



obtained by inverting the expression ^ a,t <t = 0, A = Q 
and Ap — i!^B — Aq. The formal expressions for the 
sub-blocks G, F, G, F are defined in Appendix A.l. 

One can demonstrate that the matrix Green-Keldysh 
functions satisfy the Ward identities I|B1|) and (IB2p spec- 
ified in Appendix A. 2. Making use of eq. (|Bip we can 
rewrite the first order contribution to the action in the 
form 



iSi 



-Tr 

ie 
2toc 



a,A{F + F) 



iAo Tr 



a,0{F~F) 



■Tr 



a,{V,A}(G + G) 



e Tr 



a,V{G-G) 



Observing that F — F and introducing the current j'o 
and the charge density po in the non-perturbed state we 
further transform the action iSi and get 



iSi = -2Tr 



a,AF 



-I J d^XpoiVF-Vs). 



(17) 



In addition, assuming the non-perturbed system state to 
be in thermodynamic equlibrium, we set both the current 
and the charge density equal zero. Then we obtain 



iSi = -2Tr 



(18) 



Finally, we assume that the equilibrium value of the order 
parameter Ao(r) satisfies the standard BCS gap equa- 
tion 



tanh 



2T 



(19) 



where is the Debye frequency. In this case iSi is 
cancelled by the first order contribution coming from iS^ 
and the action does not contain the first order terms in 
SA any more. 

The Ward identities (|Bip and ()B2p also allow one to 
transform the contribution iScm (|14l) and cast it to the 
form 



iSc 



iSj + iSL + iSo, 



where the terms iSj, iSl and iSo define the terms of 
a different physical origin which we will identify respec- 
tively as Josephson, London and Drude contributions to 
the effective action. They read 



iSj = 2e^ Tr 



a.^Fa.^F 



(20) 



iSr 



tt2 tt2 

d^X 



2Tr 



iSD=i I d'^X 
1 



<T.{V,-0s}F(7,{V,t)5}i^ , (21) 



Ei 



Sir 



■Tr 



a,{\7,u}G]\(a,{\7,u}dY 



2a^{V, u}Fa^{V, u}F ~ 4ma^u^{G - G) 



where 



u = 



2m 



dt' WV 



to 



,(22) 



(23) 



At low frequencies and temperatures the Josephson con- 
tribution iSj can be large, thus suppressing fluctuations 
of the gauge invariant potential In this case one can 

set $ = and get (j) = 2eV, which is just the well 
known Josephson relation between the phase and the 
electric potential. Note that for ultra-thin superconduct- 
ing wires the Josephson relation can be violated, as it 
will be demonstrated below. 

The London contribution ISl is responsible for the 
screening of the magnetic fleld penetrating inside the su- 
perconductor. Finally, the Drude contribution iSd re- 
mains non-zero in the normal state where it accounts for 
Ohmic dissipation due to flowing electric currents. Since 
the correction to the magnetic susceptibility in normal 
metals is usually small, one can ignore the vector po- 
tential in the expression for u ([23|) . Afterwards one can 
again apply the Ward identities (|B2p and rewrite iSo in 
the form 



d^X^, 



El 



^Tr 
2 



,VG 



,VG 



2a^VFa^VF 



(24) 



At last, let us consider the cross term iS^pA- Again 
applying the identity (jBl|l we cast this term to the form 
similar to eq. fTT]) : 



iS^A = - Tr 



a, A [SF-SF 



d^X 6MAf -AB)-i d^X SpoiVp -Vb), 
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where Sj,5po oc A^, define the first order corrections 
to the current and the charge density due to fluctuating 
order parameter fields. One can verify that in the absence 
of both particle-hole asymmetry and charge imbalance 
these corrections vanish. Likewise, in this case we have 
SF — SF. Thus, we conclude that 



iSy:,A — 0. 



(25) 



For clarity, we now summarize again the results derived 
this section. The complete expression for the effective ac- 
tion describing quadratic fluctuations in a superconduc- 
tor reads 



iS, 



(2) _ 
cff — 



iSi + iSA + iSj + iSl + iSo, 



(26) 



where the terms iSi, iSa, "iSj, iSl and iSo are defined 
respectively in eqs. dH]), (EOl), (ED) and 

C. Gaussian fluctuations in dirty superconductors 

Below we will mainly be interested in the limit of so- 
called dirty superconductors, i.e. we assume that the 
concentration of defects in the system is sufficiently high 
and the electron motion is diffusive. In the case of ultra- 
thin superconducting channels only this limit appears to 
be of practical interest, since usually the electron elastic 
mean free path / does not exceed the diameter of the 
wire. Since we will mainly focus our attention on wires 
with diameters in the 10 nm range, realistic values of 
I should be typically in the same range, i.e. we have 
I ^ (,Q ^ Hvf/^o ~ 1 Mm and I 'ti Xl- 

In order to account for processes with characteristic 
length scales exceeding the electron mean free path I it 
will be convenient for us to perform disorder averaging 
directly in the effective action. To this end we substitute 
explicit expressions for the Green functions (|A3|) into the 
effective action derived in the previous section and then 
apply the standard rules of averaging for the electron 
wave functions. In the diffusion approximation we have 

E ('^(^" - - ^2)Xn{r)xn{r')xm{r')Xm{r)) 



No 



ReI?(Ci-6,r-,0 



mn 

Xn(n)Xn(r2)Xm(r'3)Xm(r4)) \^^^ri=r, r3=r,=r' 

= -^^5^0 Re [z(6 - - 6, r-, r')] . 



(27) 



Here D = vpl/^ is the diffusion constant and ^{uj, r, r') 
is the diffuson defined as a solution of the diffusion equa- 
tion 



{-iu - DVl) V{lo, r, r') = S{r - r'). 



(28) 



In the following we will mainly consider spatially ex- 
tended systems in which case one has 



P(w,r,r') 



^3 giq(r-r') 



(27r)3 -iuj + Dq^' 



(29) 



Employing eqs. (|27p we arrive at the following expres- 
sion for the effective action 



E-E+ H 
An 



+ I j d'^Xd*X'|A-(X)i^f'^'A+(X') 



E-{X)K^'^' E+{X'] 



- j ct^xd'^x'i^^-{x)kl'^' ^-{x') 

+ ^-{x)kf'''^-{x')-v-s{x)k^''''vs{x') 

+ E-{X)k^^'''E-{X')^, (30) 

where we introduced "classical" A+ = {/^f + ^b)/'2 and 
"quantum" A~ = Ap — A^ components of the order 
parameter field and used analogous definitions for other 
fluctuating variables v'^ ^ and E'^ . The four kernels 

K^'^ {j = A, L, J, D) are defined as follows 



x.x' 



K 



kf-' 



(2^)4 



xm, 



d'^Q iu(t-t')+tq(r-r') Im [Xj{Q)] 



{2nY 



coth 



2r' 

(31) 



where we denote d^Q — duod^q. 

Explicit expressions for the functions xa, XJj Xl and 
Xd are rather cumbersome. They are presented in Ap- 
pendix A. 3 respectively in eqs. (ICip . (|C2p . (jC3P and 
(|C4[) . Here we provide simple analytical expressions valid 
in some limiting cases. 

Let us first concentrate on the low temperature limit 
T <C Aq. In this limit at small frequencies and wave 
vectors w, Dq^ <C Aq we find 



XA = -2iVo 1 



XJ = 2e'No 1 + 









^ 8Ao , 




TrDq^\ 


6A2 


8Ao J 


■0, 





XD 



4Ao ' 

while in the limit \uj\,Dq^ > Ao(0) = Trg-Tc (7 



(32) 
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0.577 is the Euler constant) one finds 



XA 



XJ 



-27Voln- 



-ibj + Dq^ 
Ao(0) 



XL 



XD 



uj In : — — iVq In 

Ao(0) ^ 

m'^NoDAliT) 
tj2 + 1)2^4 

In 



Ao(0) 



io; In 



2Ao(0) 



2Ao(0) 



-iuj + Dq^ ' 



(33) 



where a = 2e^NQD is the normal state Drude conduc- 
tivity. Here we expUcitly indicated the temperature de- 
pendence of the superconducting gap Aq(T) in order to 
emphasize that these asymptotic expressions are valid at 
all temperatures rather than only in the limit T <C Aq . 

Let us now consider higher temperatures |T — Tc\ <C 
Tc- At T > Tc our general expression for XA reduces to 
the standard result 



, T /I -tLO + Dq^ 

In h * - ^ — 

Tc \2 AttT 



(34) 



Xd is again defined by the Drude formula (j33p . while two 
other ^-functions vanish identically in this temperature 
interval, XJ — XL — 0. The latter observation implies 
that phase fluctuations remain unrestricted in this case. 
Hence, no Taylor expansion of the action in the phase 
Lp can be performed. In this case it is more convenient 
to undo the gauge transformation restoring the initial 
dependence of the action on the complex order parameter 
field and then to expand the action in this field. For 
simplicity ignoring electromagnetic fields and expanding 
the action to the second order in we find 

^S^:^ = \ I d^Xd4x'{(A-(X))*xf-^'A+(X') 

+ A-(X)i^f'^' (A+(X'))*} 

- J d'^Xd^X' {A-{X))* K^^^' A-{X'). (35) 

In the limit of small frequencies and wave vectors 
Dq^ ^ SttTc one recovers the standard expression 



XA 



-2iVoln 



T TTiVn 



Tc 



AT 



i-tUJ + Dq^), 



(36) 



which usually serves as a starting point for the deriva- 
tion of the so-called time dependent Ginzburg-Landau 
equation (TDGL) which we will address shortly below. 

Turning to temperatures below the critical one, T < 
Tc, and expanding the functions (|Cip - (|C4|) in powers 



of Aq we obtain 

7C(3) NoAl 



XA 



XJ 



27r2 
ttNo 
4T 



rp2 

{-iuj + Dq^) 



ttNo 



^0 



2T 



XL 



XD 



7C(3) e^NpAl 
27r2 r2 

irm'^NQDAl 
f 

2e^NQD 
—iLU + Dq'^ \ 
ire^NoDAl 



1 



^0 



Dq^' 



-2iuj + Dq^ 



iuj + Dq^ 

-2iuj + Dq^ 
-iuj + Dq^ 

7C(3) Al 
27r2 r2 

1 



{-iuj + Dq^y 

A^ ^ 
i-iUJ + Dq^f I 



1 



(37) 



2T {-iuj + Dq^y 

These results apply for Aq ^ \ix}\,Dq^ and in the limit 
of small wave vectors and frequencies, Dq^ <C 2ttTc. 
Here Aq (T) obeys the standard ECS self-consistency gap 
equation at T ~ Tc'. 



In 



Tc _ 7C(3) Al 
87r2 T2 



T 



O 



A4 
y4 



(38) 



For |w|, Dq^ <^ Aq we obtain non-analytic expressions. 
For example, XJ in this limit reads 



XJ 



^jVpAo f 7C(3) Aq 
2T \ 7r2 T 



D'^q'^ 



In 



yfj^Tly^+Dq^ 
^tj2 _^ £,2^4 _ J^q2 



UJ 



(39) 



More accurate expressions for the kernels xa, XJ 
and XL valid at temeperatures close to Tc and at any 
|w|,Dg2 <^ 2TrTc are given in Appendix, see Eqs. l|C7t - 
IC12p . Only in the hmit \oj\ < Dq"^ < Aq those ex- 
pressions match with the well known results for the co- 
efficients of the linearized (time-independent) Ginzburg- 
Landau equation: 



XA = 



XL 



7C(3) jVoA2 
27r2 T2 ' 
nrn'^NoDAl 



XJ = 



7C(3) e^NpAl 
27r2 T2 ' 



(40) 



At frequencies Dq^ ^ |w| <C Ag the functions xa and 
Xj turn out to be parametrically different taking much 
higher values: 



XA 



ttNqAq 
AT 



XJ 



ttc^NoAq 



(41) 



while XL is still given by Eq. P0)l . The Drude suscep- 
tibility Xd may be taken in the usual form p3p in both 
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cases. Thus, already at small frequencies (well below the 
gap Aq <C Tc) microscopic results can strongly deviate 
from those frequently used within semi-phenomenological 
TDGL approach. At higher frequencies and/or wave vec- 
tors \uj\ + Dq^ ^ Ao this difference becomes even more 
pronounced, cf. eqs. ([37]) . 

One can demonstrate that these kernels are not inde- 
pendent and obey the following exact identity 



XJ 



4e2A^ 



XL, 



(42) 



which directly follows from the Ward identities (|B2p . In 
addition, in the diffusive limit the kernels xj and xl are 
related to each other as 

XL(cj,g) xl(w,0) 2 X./(w,g). (43) 

This latter relation applies only for dirty superconduc- 
tors. 

For clarity, it is worthwhile to display the relation be- 
tween the x-kernels derived here and some other quan- 
tities analyzed in the literature. For example, one can 
introduce the complex conductivity cr((5) of a supercon- 
ductor (sal 



(44) 



where jq and Eq are the Fourier components of respec- 
tively the current density and the electric field. In Ref. 
[36t the function Q(Q) was analyzed which expresses the 
current via the vector potential: 

3q = -Q{Q)Aq. (45) 

Both Q(ff) and a(ff) are related to the kernels xl and 
X£) as follows 

o2 



2(0) = —^(0) = - 

c c 



^xl{,Q)-^"xd[Q) 



(46) 



and averaging is performed over three different stochastic 
variables ^a, defined by the pair correlators 



X Im 



e-^'^'+'i'- coth — 
Im [xa(0)], 

d^^-^.t+^,r ^^^^ ^ 

27r)4 2T 
X lm[xjiQ)+q\DiQ)], 



e"^"*+*'f coth — 
(27r)4 2T 



'Xl{Q) 



-xd{Q) 



(a(t,r)O(0,0)) ^-J ^e--*+^'''- coth^ 



x^lm[x.(0)]. 



(49) 



All other cross correlators of the above stochastic vari- 
ables are equal to zero. 

The representation (IT71) is just the result of the stan- 
dard Hubbard-Stratonovich decoupling transformation 
in the effective action ([5Dt . We have also used the iden- 
tity 

^---VF--lf .-V*--™^. (50) 
c ot e at 

Let us now find the least action paths for 5^. Setting 
the variational derivatives of the action (|48| with respect 
to quantum fields A", (^~, and A~ equal to zero we 
arrive at four different equations for the fields A+, 
V'^ and which provide the minimum for the action 
((48|) . The first equation describes fluctuations of the 
absolute value of the order parameter. It reads 



D. Langevin equations 



d^X' K^'^' A+(X') = U{X)- 



(51) 



Let us now rewrite our results in a slightly different 
manner. The effective action S^^ can be equivalently 
deflned by means of the following formula 



where 



iSe = i I d*X 



E-E+ H 
47r 



+ i 



j d'^Xd'^X'i^l^-{X)Kl'^' A+{X') 



<^-{X)Kf^' <i>+{X') + E-{X)K^-^ E+{X') 



v-{X)K^-'''v+{X') 



i J d^X [CaA- + + ^LVg] ., 



(48) 



The second one is the continuity equation for the super- 
current. We obtain 

^+Vjs = — 77- + — VCl, (52) 
at at m 

where we introduced the superconducting density ps and 
the superconducting current density jg 

Ps = - I d^X' Kf'^'^+iX'), 



3s 



e 

m 



d^x' K^^^'v+{X'). 



(53) 



The remaining two saddle point equations take the form 



47r 



d^X'K^^^ VE+{X')-ps{X)^UX) 

(54) 
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and 



V X ^^^ 



An 



1 dE+ ^jN_ _ 

Aire dt c c m 



Here 



3n 



dX' K 



dt' 



(55) 



(56) 



is the normal quasiparticle current. 

Eqs. dni]), dni]), dnH) ana dSS]) together with noise 
correlators (|49)) represent the set of Langevin equations 
fully describing quantum dynamics of the order param- 
eter and electromagnetic fields for dirty superconduc- 
tors within the Gaussian approximation. As it is clear 
from our derivation, these equations remain valid pro- 
vided the electron distribution function is not driven far 
from equilibrium. Generalization of our approach to non- 
equilibrium situations is also possible but will not be dis- 
cussed here. 



E. Time dependent Ginzburg-Landau equation 

Now let us establish the relation between our results 
and the approach based on the so-called time dependent 
Ginzburg-Landau equation (TDGL) which is widely used 
to model various non-stationary effects in superconduc- 
tors at temperatures close to Tc- For example, above 
the critical temperature this TDGL approach allows to 
correctly evaluate the so-called Aslamazov-Larkin fluctu- 
ation correction to the conductivity of the system. Be- 
low Tc it enables one to describe formation of phase slip 
centers and the resistive state of current biased supercon- 
ducting wires. Relative simplicity of the TDGL approach 
makes it possible to apply powerful numerical methods 
thus making this technique particularly appealing. In 
many cases the TDGL-based analysis was employed even 
far beyond its applicability range, e.g., in order to de- 
scribe quantum phase slips in superconducting nanowires 
at T ^ 0. 

The TDGL equation is usually written in the following 
simple form 0] 



^ - 2ieV + — -D 

Ot TGL 



7C(3),^,2' 



A, 



where 



TGL 



8|T-Tc| 

is the so-called Ginzburg-Landau time and 

16T2 



(57) 



(58) 



{CA{t,r)U{t' ,r')) 



5{t~t')6{r~r'). (59) 



Although this form can be justified for gapless supercon- 
ductors at high concentration of magnetic impurities, in 
general no consistent microscopic derivation of eq. ([57|) 
was ever performed. Nevertheless, it is usually believed 
that eq. ((57)) is microscopically justified at temperatures 
above Tc where the average value of the BCS order pa- 
rameter is zero, Aq = 0. Unfortunately, our present mi- 
croscopic derivation (as well as earlier imaginary time 
analysis 0, HH) does not fully support this statement. 
An independent analysis based on the real-time non- 
linear g-model was recently performed by Levchenko and 
Kamenev [37] who also noticed that even at T > Tc 
eq. ([57]) is not quite correct. These authors formulated 
a more accurate real-time TDGL equation in the form 
convenient for a comparison with eq. (|57|) : 



d „ dK,(X) 1 ^ / 2e 
— - 2^e^^ + Div + t—A{X) 

ot dt TGL V C 



+ ^|A(X)|^' 



A(X)-^A(X) I d'X'V{X,X') 



(60) 



where 1C{X) = / dX'V{X, X') [V{X') - DVA{X')\ . 

It is obvious that eq. (|60p does not in general coin- 
cide with eq. ([57]) . At the same time, it is satisfactory 
to observe that eq. ((60)) agrees with our results up to 



terms ^ A+ , . In order to demonstrate this fact it is 
necessary to identify A = (Aq -I- A'*')e*''' and consider 
terms linear in A'^ and Lp'^ . In this way one arrives at our 
Langevin eqs. ([5T|) - (|52|) with Ag =const, where the func- 
tions xA) XJi Xl coincide with those given by Eq. ([57|) 
in the leading order in Aq. Thus, we conclude that in 
the limit Aq — > our Langevin equations are equivalent 
to TDGL- type of equation [s^ within the order ~ Aq. 
Some differences, however, arise for higher order terms, 
namely for terms ~ Aq originating from the functions 
Xj and xl as well as for terms ^ Aq emerging from ^A- 
This observation indicates that eq. (|60|) is still not fully 
justified at temperatures T < Tc- 

In fact, it is not quite clear to us whether it would 
be of any practical importance to pursue the GL expan- 
sion up to terms ^ A'' in the TDGL-type equations. Of 
course, a regular expansion in powers of A* can be per- 
formed in the initial effective action ([5]). In the order 
~ A* this expansion generates many complicated non- 
local (both in space and in time) terms containing the 
quantum field (A~)" with n ranging from 1 to 4. In or- 
der to recover terms ~ A'^ in the TDGL equation one 
should disregard all terms in the action ~ (A~)" with 
n > 2. In certain situations this approximation might be 
difficult to justify. In addition, the remaining terms ~ A"^ 
are hardly tractable except in the zero frequency limit. 
Finally, the whole approach remains restricted to temper- 
atures T Tc- In view of all these problems it appears 
more appealing to perform the expansion of the effective 
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action in superconducting fluctuations around the mean 
field value Ag. This strategy was pursued in the bulk of 
this chapter. Restricting these expansion to second order 
terms in we arrive at the Langevin equations (j5ip - 
([52)1 with XA, XJ and xl defined in Eqs. (|CHIC3p . This 
approach remains applicable at all temperatures down to 
T = and is sufficient for practical calculations in a large 
number of situations. 



III. THIN METALLIC WIRES 

We now turn to the specific case of sufficiently long and 
very thin superconducting wires which will be of partic- 
ular interest for us here. For such systems the terms 
describing the action of free electro-magnetic field can be 
rewritten in the form 



wire take the form 



E-E+ 



H H+ 



Air 

I I dtdx ( CV+V- 



A+A- 



(61) 



Here we have defined the coordinate along the wire x, 
the capacitance per unit length of the wire C and the 
inductance times unit length L. A stands for the com- 
ponent of the vector potential parallel to the wire. For 
a cylindric wire with radius tq embedded in a dielectric 
environment with susceptibility e, the capacitance C and 
inductance L are 



1 92A+ ttD d^A+ 



12AI 



2e 



a<2 dt 

as d^V+ 



ttDAq 



8Ao dx^ 

2edA+\ 



dx J 



8Ao dx^ 
A 



+ 2e^Nos V+ 



_as_ / 



8Aoc \dtdx ' 



2ec 



dx 



c 

1 d'f+ 
lcPA+_ 

c a<2 

+ —A+ 
c 






= 

= 0. 



For dirty metallic wires with the diameter of the order 
of superconducting coherence length ^ = D/ Aq one 
finds 1/L » ttctAos/c^ and C <^ 2e'^NQS. In this case 
the last equation gives A^ — > 0, while the second and 
third equations describe the propagation of the plasmon 
Mooij-Schon mode [s^ with dispersion 



= caq, 

where the velocity of this mode cq is 



Co 



ttctAqs 



c 



(64) 



(65) 



and L — An\\/s — l/2ne'^NQAQDs is the kinetic induc- 
tance of a superconducting wire. 



B. Gaussian fluctuations of tlie order parameter 



21n(i?o/ro)' 



L«21n(i?oAo), 



(62) 



where i?o is the distance from the center of the wire and 
the bulk metallic electrode. 

In order to transform other terms one should apply a 
simple rule 



d^X 



dt dx, 



d^Q 



(2^)4 

where s is the wire cross section. 



1 r diodq 



(2^)2 



(63) 



A. Propagating modes 

In the low temperature limit T <C Aq all ^-functions 
((32)) are real and, hence, the noise terms in all four 
Langevin equations (|51ll55p vanish. This enables prop- 
agation of electromagnetic modes along a quasi-lD su- 
perconducting wire. The equations of motion for such a 



The effective action ([5(H) fully accounts for Gaussian 
ffuctuations in diffusive superconducting structures. For 
instance, from Eqs. (|5HI55p one readily establishes the 
correlation functions for all fiuctuating variables in our 
problem. For the order parameter fields A* in a quasi- 
lD wire we have 



(A+(t,r)A+(0,0))=-i / -^.t+^,r 
s J (27r)^ 



X Im 



1 



Xa(Q) 



coth 



2T' 



(66) 



Correlation functions v'^ , $^ and are defined analo- 
gously via the corresponding kernels xj, Xl and xd- 

Consider Gaussian fluctuations of the order parameter 
in thin one-dimensional wires. The simplest possible av- 
erage (A+(0, 0)A+(0,0))/Aq is divergent since the func- 
tion XA grows very slowly at large tu and q. Let us define 
and analyze a slightly different object 



(A+(1/Ao,0 A+(0,0)) 



A2 



(67) 



One can verify that, for example, the non-local kernel 
K^'^ significantly decays as long as \t\ exceeds I/Aq 
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and \r\ becomes bigger than the coherence length ^. 
Therefore, the parameter TZ provides a qualitative mea- 
sure of the ratio between the fluctuation correction to the 
current and its mean field value. 

In order to estimate the parameter TZ we note that at 
low temperatures T ^ Aq the kernel xa (IC1|) can be 
expressed in the form 

XA - ~NoFo{Lu/Ao,^q), 

where F{x,y) is a certain dimensionless function. Then 
we obtain 



n 



where 



2^^ JO 



is a numerical prefactor. Below we will demonstrate that 
the probability for quantum tunneling of the order pa- 
rameter field A in superconducting nanowires is propor- 
tional to exp[— S'qps], where Sqps ''^ 1/TZ is the action 
of single quantum phase slips (QPS). Hence, for 7^ <C I 
Gaussian fluctuations are small and QPS events are rare, 
which are important pre-conditions for the BCS mean 
field theory. On the other hand, at 7?. > 1 one enters the 
regime of strong non-Gaussian fluctuations which fully 
suppress the mean field order parameter thus driving the 
wire to a normal state. The concept of QPS also becomes 
ill-defined in this regime of strong quantum fiuctuations. 

The same conclusions can be extracted from the so- 
called Ginzburg-Levanyuk criterion. Let us consider the 
Ginzburg number Gi defined as the value (Tc — T)/Tc 
at which the fiuctuation correction to the specific heat 
becomes equal to the specific heat jump at the phase 
transition point. In the case of quasi-lD wires this num- 
ber reads M- 



Gi 



1.3 



ID 



(68) 



where Te = l/vp the elastic mean free time. Typically 
in thick wires one finds Gim ^ 1 and fiuctuations be- 
come strong only in a very narrow region close to Tc, i.e. 
at \Tc — T\/Tc < GiiB- One can also rewrite eq. (p5)) as 



Gi 



0.15 



ID 



or simply Giin ~ where is the dimensionless 

conductance of the wire segment of length ^. Thus, for 
5^ ~ 1 the width of the fluctuation region ST is compa- 
rable to Tc and the BCS mean field approach becomes 
obsolete down to T = 0. 



^2/3 



1 



^2/3 ' 



(69) 



C. Matsubara effective action 

To complete our analysis we will briefiy address the 
imaginary time (Matsubara) version of the effective ac- 
tion. Technically it is more convenient to deal with this 
form of the action provided one needs to account for 
quantum tunneling processes. This is precisely what we 
will do below when we describe quantum phase slips in 
superconducting nanowires. The calculation is described 
in details in Refs. [s^, [3l[ and is completely analogous 
to one carried out above in real time. 

Our starting point is the path integral representation 
of the grand partition function 



Z = VAVVVA 



(70) 



where Se is the Euclidean version of the effective action, 
the fluctuating order parameter fleld A as well as scalar 
and vector potentials V and A depend on coordinate x 
along the wire and imaginary time r restricted to the in- 
terval Q < T < (3 = 1/T. As before, assuming that devi- 
ations of the amplitude of the order parameter field from 
its equilibrium value Aq are relatively small we expand 
the effective action in powers of 5A{x, r) = A(a;, r) — Aq 
and in the electromagnetic fields up to the second or- 
der terms. The next step is to average over the random 
potential of impurities. After that the effective action 
becomes translationally invariant both in space and in 
time. Performing the Fourier transformation we obtain 

Sim 



+ X.J 



dujdq 



ILU 



Ls 
2 



XD 



XL 

Aw? 



c 



qV 
2 



-A 

c 



Xa\SA\' 



(71) 



The functions xa, Xj, Xl and xd, which depend both 
on the frequencies and the wave vectors, are expressed 
in terms of the averaged products of the Matsubara 
Green functions [s^, l3ll |. These functions represent the 
imaginary time version of the analogous real time func- 
tions Xj already encountered above. In order to recover 
the expressions for Xji'-^jk) one just needs to substitute 
-iuj |w| in eqs. (IUT|l - ((r^ for Xji^, k). The action Se 
(j7ip represents the imaginary time analogue of the real 
time effective action ([3D]). 

Note that the action Se ||7T|) is quadratic both in the 
voltage V and the vector potential A. Hence, these vari- 
ables can be integrated out exactly. Performing this in- 
tegration one arrives at the effective action which only 
depends on ip and SA. We obtain 



diudq 
(2^ 



{H^,qM'+XA\SA\'}. (72) 



The general expression for JF(w, q) and the saddle point 
relations between the electromagnetic potentials and the 
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fluctuating phase Lp are presented in Appendix A. 4 (eqs. 
(|DT|) - ((D3l) ) for completeness. 

As we already discussed, usually the wire geometric in- 
ductance remains unimportant. Therefore here and be- 
low we put L = 0. Then eqs. (|Dip - (|D3|) get simphfied 
and read 



where Fbcs is the standard BCS free energy. The or- 
der parameter is defined by the saddle point equation 
dF/dA = and can be written in the form A = 
Ao — SAq, where Aq is the solution of the BCS self- 
consistency equation dFBcs/dAo — (fTO]) and the fluc- 
tuation correction 6Aq has the form 



XL 



q [f + Xoq 



X.lXL „2 



T + X.7 + Xoq^ 



and 



V = 



XJ 



T+XJ 



xsq" 



2e 



A = 0. 



(73) 



(74) 



Note that according to eq. (|7^ the Josephson relation 
V — ip/2e is in general not satisfied. This relation may 
approximately hold only in the limit xj ^ C/s + xoq^- 
Making use of the results presented in Appendix one eas- 
ily observes that in the important limit of small elastic 
mean free paths I the latter condition is obeyed only 
at low frequencies and wave vectors uj/Aq <^ 1 and 
Dqy Ao < 1. 

Let us now perform yet one more approximation and 
expand the action in powers of to and q^. Keeping the 
terms of the order q'^ and uj'^q'^ we find 



Se — - 



8Ao 



uj q 



2iVo( 1 



UAl 



(75) 



The term 



turns out to be equal to zero. At even 



smaller wave vectors, Dq^ /2A <^ 2C /ttc^Nqs <C 1, we 
get 



Sf, — — 



dujdq 



'Cu? -I- TTcrAosg^ 



■sxA\&Af 
(76) 



Here, as before, we have assumed C/2e'^NQS <^ 1. The 
first two terms in this action correspond to the effective 
Hamiltonian of the form 



dx[{dr'pf/2+{d^^f/2] 



which again defines the Mooij-Schon plasma modes prop- 
agating along the wire with the velocity cq (|65p . 

The effective action (|72p allows to directly evaluate the 
fiuctuation correction to the order parameter in super- 
conducting nanowires. Performing Gaussian integration 
over both Lp and 6 A we arrive at the wire free energy 
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Making use of the above expressions for the functions 
F{lu, q) and xa('j-', q) and having in mind that for a wire 
of length X one has dF/dAo = 2NosX, at T ^ we 
obtain 



5Ao 



(79) 



In eq. ((79)) fluctuations of both the phase and the abso- 
lute value of the order parameter give contributions of the 
same order. The estimate ([75]) again demonstrates that 
at low temperatures suppression of the order parameter 
in superconducting nanowires due to Gaussian fluctua- 
tions remains weak as long as 3> 1 and it becomes 
important only for extremely thin wires with Giiu ^ 1. 

Finally let us return to the action ([7T|) which we will 
use to illustrate a deficiency of the TDGL approach in the 
imaginary time. Considering the superconducting part of 
the action only and assuming that temperature is close 
to Tc we can identify ^ A with A and set Aq = in all 
the x-kernels. Exactly as in the real-time approach one 
then has Xj = XL = and the phase fluctuations become 
unrestricted. For this reason one should again undo the 
gauge transformation and return to the complex order 
parameter field. 

For simplicity let us ignore both the scalar and the 
vector potentials. The TDGL action for the wire is then 
usually written in the form 

^TDGL =^oTs^ J dx^^\A\^ + Nos J drdx 



( ttD 


dA 


[w 


dx 



T-T( 



Tc 



^'A\^ + 0(\Af) 



(80) 



This form can be obtained from the action (j7ip by for- 
mally expanding the kernel XA in Matsubara frequencies 
and wave vectors w^, Dq^ <C AttT, cf. eqs. ([M]) and (|36p . 
Note, however, that since the validity of the GL expan- 
sion is restricted to temperatures T ~ Tc, the Matsubara 
frequencies a;„ = 27rnT are never really smaller than AttT. 
Hence, the expansion *(l/2 + |w„|/47rr) - ^'(1/2) 
TTUJn/8T which yields ([80]) is never correct except in the 
stationary case uj^ = 0. Already these simple arguments 
illustrate the failure of the TDGL action ([^ in the Mat- 
subara technique. Further problems with this TDGL ap- 
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proach arise in the presence of the electromagnetic po- 
tentials V and A. We refer the reader to the paper [31[ 
for the corresponding analysis. 



IV. THERMALLY ACTIVATED PHASE SLIPS 

As we already discussed, sufficiently thin supercon- 
ducting wires can acquire non-zero resistance even be- 
low Tc due to fluctuations of the superconducting order 
parameter. In this section we will address thermal fluctu- 
ations which are particularly important in the immediate 
vicinity of Tc- 

The theory of thermally activated phase slip s (TAPS) 
was developed by Langer and Ambegaoka r llal and then 
completed by McCumber and Halperin [i7|. Here we 
will briefly review this LAMH theory with minor mod- 
ifications related to the fact that the TDGL-based ap- 
proach is not sufficiently accurate to correctly determine 
the pre-exponent in the expression for the TAPS rate. 

This rate "/taps is defined by the standard activation 
dependence 



n) 



where 5F is the free energy difference which determines 
an effective potential barrier which the system should 
overcome in order to create a phase slip. 



A. Activation exponent 

In order to evaluate 5F we make use of the general ex- 
pression for the Ginzburg-Landau free energy functional 
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The saddle point paths for this functional are determined 
by the standard GL equation 



ttD d^A T~Tc 



8T dx^ 



Tc 



A 



7C(3) 
87r2T2 



ApA = 0. (83) 



In the case of a long quasi- ID wire this equation has two 
solutions, a trivial one 



Ao 



I8tt'^T{Tc -T) 



7C(3) 



(84) 



providing the minimum for the free energy and a 
metastable one 



Am (x) = Ao tanh [ ^ ^^'^^^ x 



(85) 



The potential barrier 5F in eq. ((8T|) is set by the differ- 
ence 



167r2 



SF = F[Am{x)] - F[Ao] = ^^f^oV^iTc - Tf'^ ■ 



(86) 



Note that this result applies as long as the transport cur- 
rent / across the wire is sufficiently small. With increas- 
ing / the height of the potential barrier 5F decreases and 
finally vanishes as / approaches the critical (depairing) 
current of the wire. The corresponding expression for 
SF{I) can be found in Refs. [3, Recently a mi- 

croscopic calculation of SF in the case of a clean single 
channel superconducting wire was reported in Ref. [39| . 



B. Pre-exponent 

Now let us turn to the pre-exponent B in the expres- 
sion for the TAPS rate (|8ip. In order to evaluate B one 
should go beyond the stationary free energy functional 
([82| and include time-dependent fluctuations of the or- 
der parameter field A(x,t). In Ref. [l3] this task was 
accomplished within the framework of a TDGL-based 
analysis. Employing TDGL equation it is possible to 
re-formulate the problem in terms of the corresponding 
Fokker-Planck equation 40] which can be conveniently 
solved for the problem in question. Since the important 
time scale within the TDGL approach is the Ginzburg- 
Landau time TcL (|58p . this time also enters the expres- 
sion for the pre-exponent B derived in [l7| . 

Unfortunately, as it was demonstrated, e.g., in Chap- 
ters 2 and 3, the TDGL approach fails below Tc- Hence, 
one should employ a more accurate effective action anal- 
ysis. The microscopic effective action for superconduct- 
ing wires is rather complicated and it cannot be easily 
reduced to any Fokker-Planck-type of equation. For this 
reason, below we will take a somewhat different route [4l| 
and combine our effective action approach with the well 
known general formula for the decay rate of a metastable 
state (see, e.g., (2l[27t) 



7(r) = 2ImF(T). 



(87) 



As our effective action does not contain the parameter 
TcL we expect that our final result for the pre-exponent 
will not contain this parameter either. 

Following the standard procedure [4^ we expand the 
general expression for the effective action around both 
saddle point solutions (|84|) and ((85|) up to quadratic 
terms in both the phase if and 6 A. Neglecting the contri- 
butions from fluctuating electromagnetic fields we obtain 



S, 



F[Ao,m]+5^S, 



O/M, 



(88) 
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where 



dxdx' 



X [SA{Un, x)x^^\\uJn\; X - x')dA{uJn, x') 
+ (p{uJn, x)k^^\\uJn\; X - x' )(p{uJn, x')] , 



S^Sm 



sT 



2 ^ 



dxdx' 



X [5A{uJn,x)ji-^'\\uJn\;X, x')5 A{uJn, x') 

+ ^{0Jn,x)k'^^'\\uJn\]X,x')ip(u:n,x')\. (89) 

Here w„ = 2TrTn are Bose Matsubara frequencies. The 
functions and k!p^ are expressed in terms of the ker- 
nels XA, Xj and Xi as follows: 



2^ 
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(90) 



The functions Xa^' and A:^^^' describe the fluctua- 
tions around the coordinate dependent metastable state 
Aj\/(x), and, therefore, cannot be straitforwardly related 

to XA, Xj and xl- Fortunately, the explicit form of Xa^'' 
and k^'' is not important for us here. 



The pre-exponent B in eq. ([8T|l is obtained by integrat- 
ing over fluctuations 6 A in the expression for the grand 
partition. One arrives at a formally diverging expression 
which signals decay of a metastable state. After a proper 
analytic continuation one arrives at the decay rate in the 
form dSIl) with 



B = 2TIm Yl 



\ 



detxA''(^«) detfc^°^(a;n) 
detxA''^''(w„) detfc^"'*^^(w„) 



(91) 



Here it is necessary to take an imaginary part since one 
of the eigenvalues of the operator fcM(O) is negative. 

The key point is to observe that at T Tc all Matsub- 
ara frequencies |aj„| = 27rr|n| - except for one with n — 
- strongly exceed the order parameter, \u>„\ 3> Ao(T). 
Hence, for all such values the function xa(*|w„|) ap- 
proaches the asymptotic form (|34p which is not sensi- 
tive to superconductivity and we have detxA''(w„) ~ 



detfc^*^^(a;„). The cor- 



detxA^''(^ra) and det k^\ujn) 
responding determinants in eq. ([9T|) cancel out and only 
the contribution from LUn = remains. It yields 



B ~ 2TIm 



\ 



detxi"^(0) det4"^(0) 



detxi^''^(0) detk'J">{0) 



(M), 



(92) 



The ratio of these determinants can be evaluated with 
the aid of the GL free energy functional ([5^ with the 
result pL7j] 
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(93) 



where SF is the free energy barrier (j86p . X is the wire 
length and ^(T) = y/TTD/4{Tc - T) is the superconduct- 
ing coherence length in the vicinity of Tc- 

Combining all the above results we arrive at the TAPS 
rate 



"fTAPS 




exp 



5F_ 

T 



(94) 



As we expected, this result (|94l) does not contain the 
Ginzburg-Landau time tgl and exceeds the correspond- 
ing expression [l3| by the factor ^ Ttgl- 

Eq. ([94)1 is applicable sATc — T <^Tc and as long as 
5F » T. Combining these two inequalities with eq. ([86|) 
we arrive at the condition 



(95) 



Gh-D < < 1, 



where the Ginzburg number GiiD is defined in eq. 
The double inequality ([M]) is standard for the GL theory. 
Obviously, it also restricts the applicability range of the 
LAMH theory. 



C. Temperature-dependent resistance and noise 

Every phase slip event implies changing of the super- 
conducting phase in time in such a way that the total 
phase difference values along the wire before and after 
this event differ by ±2tt. Since the average voltage is 
linked to the time derivative of the phase by means of 
the Josephson relation, {V) — ((p/2e), for the net voltage 
drop across the wire we obtain 



v^- -[r2.(/)-r_2,(/)] 

e 



(96) 



where r±2ir are the TAPS rates corresponding to the 
phase changes by ±27r. In the absence of any bias current 
/ ^ both rates are equal r±2ir = Itaps and the net 
voltage drop V vanishes. In the presence of a non-zero 
bias current the symmetry between these two rates is 
lifted since - in complete analogy to the case of Joseph- 
son junctions (cf., e.g., [2l|) - the potential barrier for 
these two processes differ. As long as the bias current / 
is sufficiently small, we obtain 

^±7r//2eT 



r±27r(^) = iTAPSd 



(97) 



Thus, at such values of / and at temperatures slightly 
below Tc the I — V curve for quasi- ID superconducting 
wires takes the form 



V = -■^TAPS sinh^, 
e Zel 



(98) 
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wit h It A PR defined in eq. (|94|) . This important result 
[161 [TtI impHes that thermal fluctuations effectively de- 
stroy long range phase coherence in the system and the 
wire acquires non-vanishing resistance R = V/I even be- 
low Tc- This resistance demonstrates strong (exponen- 
tial) dependence on temperature and the wire cross sec- 
tion 



R{T) cx exp(-5F(r)/T), 



(99) 



leading to effective fluctuation-induced broadening of the 
superconducting phase transition which can be detected 
experimentally. The corresponding discussion is pre- 
sented in Chapter 6. 

To complete our description of thermal fluctuations 
in superconducting wires we point out that in addition 
to non-zero resistance (jHH) TAPS also cause the volt- 
age noise below Tc- Treating TAPS as independent 
events one immediately concludes that they should obey 
Poissonian statistics. Hence, the voltage noise power 
Sv = {VV) is proportional to the TAPS rate ^taps- 
The contributions from TAPS changing the phase by ±27r 
add up and we obtain 



47r2 tt/ 
= — IT APS '^o^'i^^- 



(100) 



Similarly to the wire resistance the voltage noise rapidly 
decreases as one lowers the temperature away from Tc- 
Only in the vicinity of the critical temperature this noise 
remains appreciable and can be detected in experiments. 



V. THEORY OF QUANTUM PHASE SLIPS IN 
SUPERCONDUCTING NANOWIRES 

As temperature goes down thermal fluctuations de- 
crease and, hence, TAPS become progressively less im- 
portant and eventually die out in the limit T — > 0. At low 
enough temperatures quantum fluctuations of the order 
parameter field A take over and essentially determine the 
behavior of ultra-thin superconducting wires. As we have 
already discussed, the most important quantum fluctua- 
tions in such wires are Quantum Phase Slips (QPS). Each 
QPS event involves suppression of the order parameter 
in the phase slip core and a winding of the superconduct- 
ing phase around this core. This conflguration describes 
quantum tunneling of the order parameter fleld through 
an effective potential barrier and can be conveniently 
described within the imaginary time formalism. Below 
we will elaborate on the microscopic theory of quantum 
phase slips in superconducting nanowires. In doing so, 
to a large extent we will follow the papers [H, [2^, [sO] • 



A. QPS action 

Let us denote the typical size of the QPS core as xq 
and the typical (imaginary time) duration of the QPS 



event as tq. At this stage both these parameters are 
not yet known and remain to be determined from our 
subsequent analysis. It is instructive to separate the total 
action of a single QPS Sqps into a core part Score around 
the phase slip center for which the condensation energy 
and dissipation by normal currents are important (scales 
X < xq, t < To), and a hydrodynamic part outside the 
core S'out which depends on the hydrodynamics of the 
electromagnetic fields, i.e. 



S, 
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(101) 



Let us first evaluate the hydrodynamic part 5out- This 
task is simplified by the fact that outside the core the 
absolute value of the order parameter field remains equal 
to its mean field value Aq, and only its phase ip{x,T) 
changes in space and time. Without loss of generality we 
can assume that the absolute value of the order parame- 
ter is equal to zero at r = and a; = 0. For sufficiently 
long wires and outside the QPS core the saddle point so- 
lution corresponding to a single QPS event should satisfy 
the identity 



dxdr^ — drdx^ — 2'kS{t, x) 



(102) 



which follows from the fact that after a wind around the 
QPS center the phase should change by 2tt. In a way QPS 
is just a vortex in space-time with the phase distribution 
(P{x,t) described by the saddle point solution 



(p{x,T) ~ — arctan(x/coT). 



(103) 



Substitutung the solution (|103p into the action (|75p we 
obtain 

•Sout = M ln[min(co/3, X)/ max(coTo, xq)] , (104) 
where the parameter 




(105) 



sets the scale for the the hydrodynamic contribution to 
the QPS action. Here and below a — e^/hc ~ 1/137 is 
the fine structure constant. We also note that at T ^ 
the contribution ^out (|104p diverges logarithmically for 
infinitely long wires thus making single QPS events un- 
likely in this limit. 

Let us now turn to the core contribution to the ac- 
tion of a single QPS. In order to exactly evaluate this 
contribution it is necessary to explicitly find the QPS 
saddle point of the full non-linear effective action. This 
is a formidable task which can hardly be accomplished in 
practice. On the other hand, this task is greatly simpli- 
fied if one is aiming at estimating the term ^coro up to a 
numerical prefactor of order one. Below we will recover 
the full microscopic expression for the core contribution 
•Score leaving only this numerical prefactor undetermined. 
In this way we fully capture all essential physics of QPS. 



17 



The dimensionless prefactor can be regarded as a fit pa- 
rameter which can be extracted, e.g., from the compari- 
son with available experimental data. 

The above strategy allows us to approximate the com- 
plex order parameter field inside the QPS core by two 
simple functions which should satisfy several require- 
ments. The absolute value of the order parameter 
|A(a::,T)| should vanish at a; = and r = and coin- 
cide with the mean field value Aq outside the QPS core. 
The phase ip{x, r) should flip at a: = and r = in a way 
to provide the change of the net phase difference across 
the wire by 27r. On top of that, in a short wire and out- 
side the QPS core the phase (p should not depend on the 
spatial coordinate in the zero bias limit. All sufficiently 
smooth functions obeying these requirements can be used 
to estimate S'coro- For concreteness, in what follows we 
will choose 

\dA{x, t)\ = Ao exp(-xV2a;^ - T^f^To). (106) 
for the amplitude of the order parameter field and 



— tanh ( — - 
2 \xqt 



(107) 



for its phase. Rewriting the action ((75)) in the space-time 
domain 



32e2Ao ydxdTj 



sNo I dx dr 



(where we dropped unimportant terms oc {Dq^Y) and 
substituting the trial functions (|106[) . (|107[) into the ac- 
tion (jlOSp one finds 
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(109) 



where aj are numerical factors of order one which depend 
on the precise form of the trial functions, C = CX is the 
total capacitance of the wire and X is the wire length. 
Note that fictitious divergencies emerging from a singular 
behavior of the function (jl07p at x = xq and t — tq 
are eliminated since the order parameter vanishes in this 
space-time point. 

Let us first disregard capacitive effects neglecting the 
last term in eq. (jl09p . Minimizing the remaining action 
with respect to the core parameters xq, tq and making 
use of the inequality C /e^N^s <^ 1, we obtain 
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0205 2e27VoA^ 



0305 2e-^NoDA^ ' 



(110) 



These values provide the minimum for the QPS action, 
and we find 



5coro = 2sNQ{yJa2a^DAo + ^ aia^(TAo/2e-^No). (Ill) 

Substituting the Drude expression for the normal conduc- 
tance of our wire a = 2e^NaD into eqs. (|110p and (jllip 
we arrive at the final results for the core parameters 




Xo 



and for the core action 



To 



\ 1/4 1 

a2a4 \ 1 



(112) 
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^corc = ttANos^DAo = A-^- = -g^. (113) 

Here A — 2(^0203 -I- y/a^a^)/?: is the numerical prefac- 
tor Rj^ is the total normal state wire resistance, Rq — 
TTh/2e'^ = 6.453 kfi is the "superconducting" resistance 
quantum, ^ — y^D/Ao is the superconducting coherence 
length and = 4:{Rq/ Rn){X/(^) is the dimensionless 
normal conductance of a wire segment of length ^. 

As it was already pointed out, the results pi2p and 
(I113P hold provided the capacitive effects are small. This 
is the case for relatively short wires 



e^Nos 

c ■ 



(114) 



In the opposite limit the same minimization procedure of 
the action (|109p yields 



Xo 



Aro - ^XC/^e-'Nos » 1. 



(115) 



The QPS core action then takes a somewhat more com- 
plicated form 



A 



Rn 



3/2 
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(116) 



where A' is again a numerical prefactor. Eqs. (|10ip . 
(|104p . (|113p and (|116p provide complete information 
about the action for single QPS in diffusive supercon- 
ducting nanowires. 

Let us analyze the above expressions. Introducing the 
number of conducting channels in the wire A/" — pps/Air, 
setting C ^ 1 and making use of the condition /hvp ^ 
1 satisfied for typical metals, one can rewrite the inequal- 
ity (jll4p in a very simple form 



X<^iAf. 



(117) 



To give an idea about the relevant length scales, for typ- 
ical values ^ ^ 10 nm and N ^ 10^ 10'^ - according 
to eq. (|117p - the wire can be considered short pro- 
vided its length does not exceed 1 10 /xm. This condi- 
tion is satisfied in a number of experiments, e.g., in Refs. 
32 , IS [4^ lisl. [4^ . On the other hand, in experiments 
471 . |48| . I49I |50| much longer wires with lengths up to 
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X ^ 100 ixm were studied. Apparently such samples are 
effectively in the long wire regime X ^ . 

In the case of short wires we observe a clear separation 
between different fluctuation effects contributing to the 
QPS action: Fluctuations of the order parameter field 
and dissipative currents determine the core part (|113|) 
while electromagnetic fluctuations are responsible for the 
hydrodynamic term (|104p . In the case of longer wires ca- 
pacitive effects also contribute to the core part (|116p . We 
also observe that in the short wire limit different contri- 
butions to the QPS action depend differently on the wire 
thickness (or the number of conducting channels N) : The 
core part (|113p decreases linearly with the wire cross sec- 
tion, 5'core oc A^, whcrcas the hydrodynamic contribution 
shows a weaker dependence /i (x ^/JI . In the long wire 
limit the dependence of the core part (|116p on M also 
becomes weaker, Score VA/", due to capacitive effects. 

Yet another important observation is that in the inter- 
esting range of wires thicknesses rg = y/s > 5-^10 nm the 
core part Score usually exceeds the hydrodynamic term 
'-^ fi. E.g. for C ~ 1 we obtain fi w 30(ro/Ai). Setting 
ro ~ 5 ^ 10 nm and estimating /i and Score for typical 
system parameters Pp^ ~ 0.2 nm, I ^ 7 nm, ^ 10 nm, 

XI ~ 100 nm, vp = 10^ m/s and Aq ~ 1 10 K we find 
/i 1 -f- 3 and Score ^ The latter inequality becomes 
even stronger for thicker wires. Note that the condition 
Score 3> /i allowed us to ignore the hydrodynamic part 
of the QPS action while minimizing the core part with 
respect to xq and tq. 

At the first sight the result pi3p derived in the short 
wire limit could create an illusion that our microscopic 
description would not be needed in order to recover the 
correct form of the core part Score- Indeed, the same form 
could be guessed, e.g., from an oversimplified TDGL- 
based approach or just from the "condensation energy" 
term (proportional to xa) without taking into account 
dissipative effects. For instance, minimization of the con- 
tribution ^ |<5Ap (the last three terms in eq. (|108p ) 
is formally sufficient to arrive at the correct estimate 
Score ~ NoS\/DAo. The same equation demonstrates, 
however, that not only the amplitude but also the phase 
fluctuations of the order parameter field provide impor- 
tant contributions to the QPS action. If the latter fluc- 
tuations were taken into account without including dissi- 
pative effects (this would correspond to formally setting 
CT — > in eq. (jlOSp ) minimization of the core action 
would immediately yield the meaningless result xq —^ 0, 
tq (cf., eq. (jllOp ') implying that the hydrodynamic 
contribution Sout could not be neglected in that case. 
Minimization of the total action Sgps would then yield 
the estimate for the action parametrically different from 
that of eq. (|113p . On the other hand, in the strong damp- 
ing limit (approached by formally setting cr ^ cxd in eq. 
(jlOSp ) the core size would become very large and the 
action Score (|llip would diverge meaning that no QPS 
would be possible at all. These observations clearly illus- 
trate crucial importance of dissipative effects. Under the 
condition C/c^Nqs <C 1 (usually well satisfied in metal- 



lic nanowires) dissipation plays a dominant role during 
the phase slip event, and the correct QPS core action 
cannot be obtained without an adequate microscopic de- 
scription of dissipative currents fiowing inside the wire. 
Only employing the Drude formula for the wire conduc- 
tivity a — 2e^NQD enables one to recover the correct 
result (|113p whereas for some other models of dissipation 
different results for the core action would follow, cf., e.g., 
Ref. [13. 

To complete our discussion of the QPS action let us 
recall that in the course of our derivation we employed 
two approximations: (i) we expanded the action up to 
the second order in SA{x,t) = A{x,t) — Aq and (ii) 
in eq. ((75|) we expanded the action ((72)) in powers of 
uj/Aq and Dq^ /Aq. The approximation (i) is sufficient 
everywhere except inside the QPS core where A{x,t) is 
small. In these space- and time-restricted regions one can 
expand already in A{x,t) and again arrive at eq. (j72p 
with SA{x,t) — > A{x,t) and with all the ^-functions 
defined in Appendix A3 with Aq = 0. Both expansions 
match smoothly at the scale of the core size xq ~ ^, 
To ^ I/Aq. Hence, the approximation (i) is sufficient to 
derive the correct QPS action up to a numerical prefactor 
of order one. 

The approximation (ii) is sufflcient within the same ac- 
curacy. One can actually avoid this approximation and 
substitute the trial functions (|106p . (|107p directly into 
the action ([7^ . Neglecting capacitive effects in the limit 
pi4p one can rewrite the QPS action as a function of the 
dimensionless parameters xq/S, and AqTq only. Making 
use of the general results for the x-functions collected in 
Appendix A3 and minimizing the QPS action with re- 
spect to Xq and To one again arrives at the result l|113p 
with A ^ 1. If the inequality (|114p is violated, the accu- 
racy of our expansion in powers of cj/ Aq can only become 
better (cf. eq. pT5)l V 

B. QPS rate 

We now proceed further and evaluate the QPS rate 
IQPS- Provided the QPS action is sufficiently large 
Sqps ^ 1 this rate can be expressed in the form 

IQPS = Bexp{-SQPs)- (118) 

The results for the QPS action derived above allow to 
determine the rate Jqps with the exponential accuracy. 
Here we evaluate the pre-exponential factor B in eq. 
(|118p . For this purpose we will make use of the stan- 
dard instanton technique [4^ . 

Consider the grand partition function of the wire Z. 
As we already discussed this function can be expressed 
via the path integral 




VAVipcxp{-S), (119) 



which will be evaluated the saddle point approximation. 
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The least action paths 



,55/|^A| = 0, 5S/5lp = 



(120) 



determine all possible QPS configurations. Integrating 
over small fluctuations around all QPS trajectories one 
represents the grand partition function Z in terms of in- 
finite series where each term describes the contribution 
of one particular QPS saddle point. Provided interac- 
tion between different quantum phase slips is sufficiently 
weak one can perform a summation of these series in a 
straightforward manner with the result 



Z = eM-F/T), 
where F defines the wire free energy 



(121) 



i?5ycxp(-<525i[<5r]) 



^]v5Ye^p{-5-^So[5Y])' 

= F - ^^^^ 
" 2 ■ 



(122) 



Here is the free energy without quantum phase slips, 
5Y = {SA, 6(p) describe fluctuations of relevant coordi- 
nates (fields), (S^S'o,i[(5y] are the quadratic in 6Y parts 
of the action, and the subscripts "0" and "1" denote the 
action respectively without and with one QPS. 

The integrals over fluctuations in eq. (|122p can be 
evaluated exactly only in simple cases. Technically such 
a calculation can be quite complicated even if explicit an- 
alytical expressions for the saddle point trajectories are 
avalilable. In our case such expressions for the QPS tra- 
jectories are not even known. Hence, an exact evaluation 
of the path integrals in eq. (|122p is not possible. Fur- 
thermore, any attempt to find an explicit value for such 
a prefactor would make little sense simply because the 
numerical value of A in eq. (|113p is not known exactly. 

What to do in this situation? Below we will present 
a simple approach which allows to establish the correct 
general expression for the pre-exponent B up to an unim- 
portant numerical prefactor. Our approach may be use- 
ful not only in the case of superconducting wires but for 
various other situations since numerical prefactors in the 
pre-exponent are usually of little importance. 

In order to evaluate the ratio of the path integrals in 
eq. (|122p let us introduce the basis in the functional 
space 5'fc(z) in which the second variation of the action 
around the instanton S^Si [SY] is diagonal. Here the basis 
functions depend on a general vector coordinate z which 
is simply z = (t, x) in our case. The first N functions 
^ffe are the so-called "zero modes" reflecting the instan- 
ton action invariance under shifts in certain directions 
in the functional space. In our case the problem has 
two zero modes corresponding to shifts of the QPS posi- 
tion along the wire and in imaginary time, i.e. N — 2. 
Obviously such shifts do not cause any changes in the 
instanton energy. The eigenfunctions corresponding to 
these zero modes are: ^'fc(X) = dY/dzk, where Y{z) the 
instanton (or QPS) trajectory, k < N and the number of 



zero modes N coincides with the dimension of the vector 
z. An arbitrary fluctuation SY{z) can be represented in 
terms of the Fourier expansion 



N 



SY{z) = ^'^k 



k=l 



dY{z) 
dzk 



J2 Uk^k{z). (123) 



k=N+l 



Then we get 



oo 

5'^So[SY] = 2 X! AknUkU^, 

k,n—l 
^ oo 



(124) 



k=N+l 



where for k < N the Fourier coefficients Uk = 5zk are 
just the shifts of the instanton position along the fc— th 
axis and Xk are the eigenvalues of 5^S'i[(5y]. Integrating 
over the Fourier coefficients one arrives at the standard 
formula for the ratio of determinants with excluded zero 
modes 1431 



JV6Yexp{~S^Si[SY]) 
JVSYexp{-6^So[SY]) 



d6xi.. j dSxjq 





det Akn 



(125) 



(27r)^ n Afe 

A;=Ar+l 



where Lk is the system size in the fc— th dimension. Now 
we will argue that with a sufficient accuracy in the lat- 
ter formula one can keep the contribution of only first 
N eigenvalues. Indeed, the contribution of the "fast" 
eigenmodes (corresponding to frequencies and wave vec- 
tors much larger than the inverse instanton size in the 
corresponding dimension) is insensitive to the presence 
of an instanton. Hence, the corresponding eigenvalues 
are the same for both and 5^ Si and just cancel out 
from eq. (|125p . In addition to the fast modes there are 
several eigenmodes with frequencies (wave vectors) of or- 
der of the inverse instanton size. The ratio between the 
product of all such modes for 5"^ Si and the product of 
eigenvalues for 5'^Sq with the same numbers is dimen- 
sionless and may only affect a numerical prefactor which 
is not interesting for us here. Dropping the contribution 
of all such eigenvalues one gets 



^V5Y eyiY>{-5^Si[5Y]) 
^V5Y eyiY>{-5^So[5Y]) 



dSxi.. j dSxN 





det^ 



kn\k,n<N 
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(126) 



What remains is to estimate the parameters Akk for 
k < N. For this purpose let us observe that the sec- 
ond variation of the action becomes approximately equal 
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to the instanton action, = ^A^kZ^^ ~ Sgps, when 
the shift in the A;— th direction becomes equal to the in- 
stanton size in the same direction Szk = zq^- Then we 
find Akk ~ 2SqpsIzI^ and 



N 



det A 



kn\k,n<N 



k=l 



kk 



N 
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N 



(127) 



Finally, combining eqs. pT5)l . and (fTTTl) we 

obtain 



N/2 



(128) 



Here b is an unimportant numerical prefactor. This for- 
mula demonstrates that the functional dependence of the 
pre-exponent can be figured out practically without any 
calculation. It is sufficient to know just the instanton ac- 
tion, the number of the zero modes N and the instanton 
effective size z^k for each of these modes. 

In fact, a similar observation has already been made 
(4^ in the case of local Lagrangians equal to the sum 
of kinetic and potential energies. Here we demonstrated 
that the result ()128|) can be applied to even more general 
effective actions, including nonlocal ones. 

Turning to the interesting for us case of QPS in super- 
conducting nanowires we set Li = 1/T, L2 = X. Then 
eq. p28| yields 



B = b 



SqpsX 



(129) 



This equation provides an accurate expression for the 
pre-exponent _B up to a numerical factor h of order one. 
This result is parametrically different from the resuls 
derived within the TDGL-type of analysis |24| or sug- 
gested phenomenologically in Ref. [23|. The inequal- 
ity Sco-cc 3> /i allows to substitute (Score (|113p instead 
of Sqps in eq. (I129|) . Substituting also xq ^ and 
To ~ I/Aq, for the QPS rate we finally obtain 



R X 

IQPS ~ An— 2---2- exp(-S'QPs) ~ g^Ao^ exp(-yl55/4). 
Hn t, 



X 
1 



(130) 

This result concludes our calculation of the tunneling rate 
for single quantum phase slips. 

Finally we would like to emphasize that the method 
employed here works successfully in various other prob- 
lems described both by local and nonlocal in time La- 
grangians. Several examples of such problems are dis- 
cussed in Ref. (3d| . Here we mention only one such 
example. It is the well known problem quantum tunnel- 
ing with dissipation [2^. In the limit of strong dissipa- 
tion quantum decay of a metastable state was treated 
by Larkin and Ovchinnikov [sH who found the exact 
eigenvalues and, evaluating the ratio of the determinants, 
obtained the prefactor in expression for the decay rate 



B (X 77''/^ /m^, where rj is an effective friction constant 
and m is the particle mass. This result would imp ly that 
the pre-exponential factor in the decay rate [51[ should 
be very large and may even diverge if one formally sets 
TO ^ 0. Later it was realized [52,] that this divergence 
should be regularized by means of proper renormaliza- 
tion of the bare parameters in the effective action. After 
that the high frequency contribution to the pre-exponent 
is eliminated and one arrives at the result [S^l B <x 1/ y/r] 
which does not contain the particle mass m at all. This 
result allowed to fully resolve a discrepancy between the- 
ory and experiments 53]. Note that the result (s^ is 
trivially reproduced from eq. (|128p : The pre-exponent 
B can also be expressed in the form B ^ y/Sb/ro, where 
Sb oc r] is the instanton (bounce) action and tq oc 77 is 
its typical size. It is remarkable that in this particular 
case our approach correctly reproduces even the numer- 
ical prefactor. 



C. QPS interactions and quantum phase 
transitions 

Although typically the hydrodynamic part of the QPS 
action (|104p can be smaller than its core part (|113p 
the former also plays an important role since it deter- 
mines interactions between different quantum phase slips. 
Consider two such phase slips (two vortices in space- 
time) with the corresponding core coordinates {xi,ti) 
and (x2, T2). Provided the cores do not overlap, i.e. pro- 
vided \x2 — xi\ > xq and \t2 — nj > tq, the core contri- 
butions are independent and simply add up. In order to 
evaluate the hydrodynamic part we substitute the super- 
position of two solutions (p{x—xi, T—Ti)+ip(x—X2, T—T2) 
satisfying the identities 

dxdr(p - drdx^P = 27ri^i^2^(T - Ti,2tX- Xi,2) (131) 

(where vi^2 — ±1 are topological charges of two QPS 
fixing the phase change after a wind around the QPS 
center to be ±27r) into the action ((75)) and obtain 



-(2) 

'qps 



IJLV1V2 In 



{xi - X2f + c^(ri - T2) 



2 1 



(132) 

i.e. different quantum phase slips interact logarithmically 
in space-time. QPSs with opposite (equal) topological 
charges attract (repel) each other. 

The next step is to consider a gas of n quantum phase 
slips. Again assuming that the QPS cores do not overlap 
we can substitute a simple superposition of the saddle 
point solutions for n quantum phase slips ip = X)" '^(^ ~ 
Xi, T — Ti) into the action and find 



(n) 
QPS 



nSr, 



An) 



where 



\ xo 



(133) 



(134) 
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Here pij = [c^iri — tj)"^ + {xi — Xj)"^)^/"^ defines the dis- 
tance between the i-th and j-th QPS in the {x, r) plane, 
Vi = ±1 are the QPS topological charges and $o = hc/2e 
is the flux quantum. In eq. (jl34p we also included an ad- 
ditional term which keeps track of the applied current 
/ flowing through the wire. This term trivially follows 
from the standard contribution to the action [2l| 



The grand partition function of the wire is represented as 
a sum over all possible configurations of quantum phase 
slips (topological charges): 



n=0 ''^0 -^0 Jxo 



dX2n 
Xo 



dri 

To 



dT: 



2n \ ^ 



To 



(2«)^ 
/ , ^'"-fV "--int 
i',=±l 



exp(-5LT^)(135) 



where an effective fugacity y of these charges is related 
to the QPS rate as 

y = exp(-S'corc) ~ S'corc exp(-S'corc) • (136) 

We also note that only neutral QPS configurations with 

n 

^ = (137) 

i 

(and hence n even) contribute to the partition function 
(|135p . This fact is a direct consequence of the boundary 
condition (p{x, r) = (fi(x, r -|- (3) in the path integral for 
the partition function [21| . 

It is easy to observe that for / = Eqs. (HSJ), 
define the standard model for a 2D gas of logarithmi- 
cally interacting charges ui. The only specific feature of 
our present model as compared to the standard situa- 
tion is that here the space and time coordinates are not 
equivalent and one can consider different limiting cases 
of "long" and "short" wires. 

Let us first consider the limit of very long wires and 
assume that T ^ 0. Following the standard anal ysis 
of logarithmically interacting 2D Coulomb gas [sl. [ol. llCll| 
which is based on the renormalization group (RG) equa- 
tions both for the interaction parameter fi and the charge 
fugacity y. Defining the scahng parameter £ = ln(p/^) 
we have d, [l3| 



^222 

-47r ^ y , 



diy = (2 - ^j.)y. 



(138) 



Following the standard line of reasoning we immediately 
conclude that a quantum phase transition for phase slips 
occurs in a long superconducting wire at T — > and 



/.! = fi* = 2 + Airy w 2. 



(139) 



This is essentially a Berezinskii-Kosterlitz-Thouless 
(BKT) phase transtion [1, [l3| for charges Vi in space- 
time. The difference from the standard BKT transition 



in 2D superconducting films is only that in our case the 
transition is driven by the wire thickness ^/s (which en- 
ters into fi) and not by temperature. In other words, 
for thicker wires with /j, > fj,* quantum phase slips with 
opposite topological charges are bound in pairs (dipols) 
and the linear resistance of a superconducting wire is 
strongly suppressed and T-dependent. This resistance 
tends to vanish in the limit T 0. Thus, we arrive at an 
important conclusion: at T = a long quasi-lD super- 
conducting wire remains in a superconducting state, with 
vanishing linear resistance, provided its thickness is suffi- 
ciently large and, hence, the electromagnetic interaction 
between phase slips is sufficiently strong, i.e. > /i*. 

On the other hand, for /i < /i* the density of free 
(unbound) quantum phase slips in the wire always re- 
mains finite, such fluctuations destroy the phase coher- 
ence (and, hence, superconductivity) and bring the wire 
into the normal state with non- vanishing resistance even 
at T = 0. Thus, another important conclusion is that 
superconductivity in sufficiently thin wires is always de- 
stroyed by quantum fluctuations. 

The above analysis is valid for sufficiently long wires. 
For typical experimental parameters, however, X < cq/T 
(or even X cq/T), and the finite wire size needs to be 
accounted for. For this purpose we modify the above 
RG treatment in the following manner. Starting from 
small scales p ~ we increase the scaling parameter £ 
and renormalize both and y according to RG equations 
(|138p . Solving these equations up to ^ = = 
we obtain the renormalized fugacity y — y[£x)- For 
larger scales £ > Ix only the time coordinate matters 
and we arrive at the partition function fully equivalent 
to one for a (OD) superconducting weak link (Josephson 
junction) in the presence of quantum fluctuations of the 
phase. These systems are described in details in Ref. 
[2l| , therefore an extended discussion of this issue can be 
avoided here. We only point out that our renormalized 
fugacity y is equivalent to the tunneling amplitude of the 
Josephson phase (normalized by the Josephson plasma 
frequency), i.e. 



y 



^Soexp(-S'o) 



(140) 



where 5*0 = ^J^Ej/Ec is the action of an instanton 
(kink) in a Josephson junction, Ej and Eq are respec- 
tively the Josephson and the charging energies. 

The subsequent analysis essentially depends on the 
presence (or absence) of additional dissipation in our sys- 
tem. In the absence of dissipation at T ^ the Joseph- 
son junction always remains in the normal (i.e. non- 
superconducting) state since superconductivity is sup- 
pressed by quantum fluctuations of the phase and tunnel- 
ing of Cooper pairs is prohibited [21]. Due to eq. (|140p 
exactly the same conclusion applies for sufficiently short 
superconducting wires studied here. In the presence of 
additional dissipation, however, quantum fluctuations of 
the phase can be suppressed and, hence, superconductiv- 
ity can be restored |21J . 
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Of particular importance is the case of Ohmic dissipa- 
tion which is realized either provided the wire is shunted 
by a normal resistor Rs, or the wire itself has a non- 
vanishing normal conductivity a — aqp even far from the 
QPS core. The latter situation can occur, e.g., at finite 
(and not too low) temperatures due to the presence of a 
sufficient number of quasiparticles above the gap or pos- 
sibly also due to nonequilibrium effects. Bearing in mind 
that in a number of experiments with ultra-thin wires (to 
be discussed below) QPS effects were observed already at 
sufficiently high temperatures T > Ao(T) it is worth to 
briefly address the model with Ohmic dissipation here. 

In this case our consideration should be modified em- 
ploying a two stage scaling procedure [2§|. At it was al- 
ready explained, we first proceed with 2D RG equations 
(fnS)) up to the scale £ = ix = ln{X/£_). For simplicity 
we assume that Ohmic dissipation does not significantly 
affect RG eqs. (|I38p at such small scales. Eventually we 
arrive at the renormalized fugacity y = y(ix)- At larger 
scales the space coordinate is irrelevant and the problem 
reduces to that of a ID Coulomb gas with logarithmic 
interaction. Therefore, (for y <C 1) further scaling is de- 
fined by [Hi m, m Hi, [13 



9iy = (1 - 7)y, da = o. 



(141) 



where the interaction parameter 7 depends on the dissi- 
pation strength being equal to either (i) 7 = Rg/Rs or 
to (ii) 7 — ■KS(jqp/2e^X. For 7 > 1 the fugacity scales 
down to zero, which again corresponds to a supercon- 
ducting phase, whereas for 7 < 1 it increases indicating 
a resistive phase in complete analogy to a single Joseph- 
son junction with Ohmic dissipation. In the case (ii) the 
phase transition point again depends on the wire cross 
section s as well as on its total length X and the value 

Thus, two different quantum phase transitions can 
occur in superconducting nanowires. One of them is 
the BKT-like phase transition which is controlled by 
the strength /i of inter-QPS electromagnetic interactions 
and eventually by the wire thickness. Another one is 
the Schmid phase transition occurring in the presence 
of Ohmic dissipation [H, 111] provided, e.g., by an ex- 
ternal shunt resistance Rs- Note that this situation is 
somewhat reminiscent of that occurring in chains of re- 
sistively shunted Josephson jun ctions or granular arrays 

[5i[63,[Ml,[63,[6llMl6i,|6i,[63,[6i,[6i,[m0 where, 
however, Schmid-like QPT is driven by local (rather than 
global) shunt resistance. 



D. Wire resistance at low temperatures 

Let us now turn to the calculation of the wire resis- 
tance in the presence of quantum phase slips. We first 
consider the limit of long wires. At any nonzero T such 
wires have a nonzero resistance i?(T, /) even in the "or- 
dered" phase n > fj,* . In order to evaluate R{T) in this 
regime we proceed perturbatively in the QPS fugacity y. 



Since for fi > fi* quantum phase slips form closed pairs 
(dipols) and, hence, interactions between different dipols 
can be neglected. For this reason it suffices to evalu- 
ate the correction SF to the wire free energy due to one 
bound pair of quantum phase slips with opposite topo- 
logical charges. This procedure is completely analogous 
to that described in details in Ref. f2l| (see Chapter 5.3 
of that paper for an extended discussion). Taking into 
account only logarithmic interactions within bound pairs 
of quantum phase slips we can easily sum up the series 
in cq. (|135p and arrive at the result 
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where p 



TO "^0 Jxo ^0 

(142) 

It is convenient to first in- 



,2Nl/2 



tegrate over the spatial coordinate x and take the wire 
length X ^ 00. For nonzero / the expression in Eq. (|142p 
is formally divergent for — > 00 and acquires an imag- 
inary part Im SF after analytic continuation of the in- 
tegral over the temporal coordinate r [21]. This indi- 
cates a QPS-induced instability of the superconducting 
state of the wire: the state with a zero phase difference 
S(p{X) = (p{X) — ip{0) = decays into a lower energy 
state with 6ip{X) = 27r. The corresponding decay rate 
r27r is defined by eq. ([87)1 . The rate for the opposite 
transition from 27r to (or from to —2tt) - which is 
nonzero as long as T > - is defined analogously with 
I ~L The average voltage drop across the wire is 
then given by the difference between these rates 



t/ = ($o/c)[r2x(/)-r2.(-/)] 



(143) 



Evaluating the imaginary part of the free energy ImSF 
from dST]), P15)) we finally obtain 

CToXo T{p)T{2fi-l) \2cT) 
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where r(a;) is the Euler gamma-function. For the wire 
resistance R{T,I) = V/I this expression yields 



for T > $0^ and 
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(146) 



for T <^ $0-^- Thus, for sufficiently thick wires with fi > 
fi* , we expect a strong temperature dependence of the 
linear resistance which eventually vanishes in the limit 
T — !■ indicating the superconducting behavior of the 
wire. 

Unfortunately for thinner wires with fi < fj,* and at 
low temperatures the above simple analysis becomes in- 
sufficient because of the presence of unbound QPS and 
the necessity to account for many-body effects in the 
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gas of quantum phase slips. We expect that the tem- 
perature dependence of the wire resistivity should be- 
come linear at the transition to the disordered (i.e. non- 
superconducting) phase. At T ^ ^ol/c we expect a 
strongly nonlinear I — V characteristics V ^ I"^ in thick 
wires, and the universal power a(/i*) = 2 in thin wires at 
the transition into the resistive state with V ^ I. Note 
that in contrast to the BKT transition in 2D supercon- 
ducting films, in the case of wires the jump is expected 
to be from a = 2 to 1 rather than from a = 3 to 1. 

In the case of short wires X <^ cq/T we again proceed 
in two steps. A 2D scaling analysis yields the "global" 
fugacity y. In analogy with resistively shunted Josephson 
junctions [2l[, the voltage drop from the imaginary part 
of the free energy reads 
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where tq ~ XC/e^^. 
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(147) 

This equation gives R oc y^T^^ ^ 



and R oc y'^P'<~'^ respectively at high and low T . The 
above result is valid for 7 > 1 and also for smaller 7 at 
not very small T [2l| . In the low temperature limit in 
the metallic phase the linear resistance becomes [2lj 



R = Sa^p/X 



(148) 



or R — Rs, i.e. R is just equal either to the quasiparti- 
cle resistance of the wire or to the shunt resistance. The 
physics of this result is exactly the same as in the case 
of resistively shunted Josephson junctions: the super- 
conducting channel turns out to be completely blocked 
by quantum fluctuations and the current can only flow 
through the normal resistance. 



E. Discussion 

It is worthwhile to emphasize that the above results 
were derived for sufhciently uniform nanowires with (al- 
most) constant cross-section s. As we already mentioned 
above, the very existence of quantum phase slips and/or 
the possibility to experimentally observe QPS effects in 
such wires were debated at an early stage of the research. 
At present there already exist a large number of indepen- 
dent experiments by different groups providing an over- 
whelming evidence that QPS effects essentially determine 
low temperature properties of ultra-thin superconducting 
wires. A detailed comparison between these experimen- 
tal results and the above theoretical predictions will be 
carried out in the next chapter. 

It is quite obvious that our QPS theory does not re- 
quire strict uniformity of the system, it will also apply 
provided there exist relatively small diameter variations 
along the wire which are practically unavoidable in any 
realistic experimental situation. The probability for QPS 
to occur will be higher in thinner parts of the wire than in 
thicker ones, otherwise all the physics remains essentially 
the same. This situation should be contrasted to that of 



chains of Josephson junctions and ID granular wires in 
which case quantum phase slips occur inside the junctions 
(tunnel barriers) and do not require any suppression of 
the superconducting order parameter. In other words, 
such quantum phase slips have no core at all. Accord- 
ingly, no condensation energy is lost during such phase 
slip events and also dissipative effects remain insignifi- 
cant (unless Josephson junctions are shunted by external 
normal resistors). 

Note that, although the QPS core physics for uni- 
form and granular quasi- ID wires is entirely different, 
the behavior of these systems at larger scales can be 
qualitatively similar. For instance, it is well known 
that the BKT-like quantum phase transition (driven by 
the ratio E j/Eg) occurs in ID superconducting granu- 
lar arrays (sol . [60 . [73j similarly to the case of uniform 
wires [28|. Also, Schmid-like quantum dissipative phase 
transitions in arrays and chains of resistively shunted 
Josephs on j unctions were studied in a great detail, see 
[H, m, m, H m, M, M, IH, IZO, IZII, lli and further 
references therein. One of the main conclusions reached 
there was that dissipative phase transitions yield local 
phase coherence in a long chain of junctions provided 
each of such junctions is shunted by an effective Ohmic 
resistance below some value of order the quantum resis- 
tance unit Rq. In contrast, at T — > BKT-like quan- 
tum phase transition yields global phase coherence in 
Josephson chains provided Ej > Eq- As a result, one 
arrives at a non-trivial phase diagram with all possi- 
ble combinations of local and global order [H, [lOl . Let 
us also mention that more recently the models of re- 
sistively shunted superconducting granular chains were 
studied [t^ , FtsI [tgI aiming to interpret some recent ex- 
perimental data [H, 13] seemingly indicating that the 
superconducting phase transition can be driven by the 
global wire resistance. Within this approach dissipation 
is introduced phenomenologically employing an effective 
two-fluid model. It remains unclear, however, what could 
be the physical origin of Ohmic dissipation in sufhcicntly 
long and nominally uniform superconducting wires at 
temperatures well below the gap Aq. An alternative pro- 
posal aimed to explain recent exp erimental observations 
in relatively short nanowires [13 13 discussed below 
in Sec. 6.3 

We should also stress that even if the diameter remains 
strictly constant the wire can be considered uniform only 
at length scales exceeding the elastic electron mean free 
path I whereas at the scale ~ I this uniformity is natu- 
rally broken by the presence of impurities and defects. 
Scattering of electrons at such impurities and defects 
plays a very important role since it breaks Galilei in- 
variance thus preventing from momentum conservation 
for electrons propagating along the wire 31]. Momen- 
tum non-conservation is in turn necessary for QPS to 
occur since quantum phase slips "unwind" the supercur- 
rent thus generating momentum transfer. At T = (i.e. 
in the absence of phonons) this extra momentum can- 
not be absorbed in perfectly clean and uniform wires in 
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FIG. 2: Superconducting nanowire containing a thinner part 
(constriction) of length Xc and cross section Sc- 



which case QPS would be strictly prohibited. On a for- 
mal level, this implies that the QPS action Sqps should 
become large in the clean limit [73|- fact, this effect 
can be observed within our analysis too, one just needs 
to formally take the limit cr —> oo in eq. (jllip . 

For completeness, let us mention about several propos- 
als concerning the nature of QPT in long superconduct- 
ing nanowires with impurities. Khlcbnikov and Pryadko 
|78| argued that this transition should occur at /i = 1 as 
a dissipative QPT ^ rather than the BKT-like QPT 
discussed here. Accordingly, for the wire resistance they 
suggested 

R cx T^/^-^ (149) 

instead of i? oc T^^-a^ (|145p . Sachdev et al. 

[t^, on the contrary, found QPT of the BKT universal- 
ity class but claimed it to be superconductor-insulator 
transition (SIT) rather than superconductor-metal tran- 
sition (SMT) as originally suggested in Ref . ^28] . 

In their analysis Khlebnikov and Pryadko [78] used a 
simplified Gross-Pitaevski action which is in many re- 
spects different from the microscopic effective action for 
superconducting wires derived in Chapter 2. Within this 
phenomenological approach the authors 78] first evalu- 
ated the QPS rate for a given disorder configuration and 
only then performed averaging over disorder. In contrast, 
within our analysis averaging over disorder is performed 
in the cflFective action before introducing QPS. This latter 
way is appropriate for diffusive superconducting wires, 
e.g., because the QPS core size ~ ^ greatly exceeds the 
electron mean free path I. The condition ^ ^ I allows 
to treat the wire as effectively homogeneous for QPS and 
the problem acquires an extra zero mode corresponding 
arbitrary shifts of the QPS core along the wire. This zero 
mode is lacking within the analysis [781 ]. 

Although for homogeneous (at scales > /) and uni- 
form wires QPT is of the BKT universality class and the 
temperature dependence of the resistance in the ordered 
(superconducting) phase is given by eq. p45p . there ex- 
ist other situations in which the dependence (|149p may 
apply at not very low T . Consider, e.g., a long supercon- 
ducting wire with cross-section s and length X which is 
uniform everywhere except a small part (constriction) of 
length Xc X and cross-section Sc < s, as shown in Fig. 
2. QPS would then occur easier inside the constriction 
than in the rest of the wire. Assuming for simplicity that 
Xc is of order of (or just slightly larger than) ^, for the 



corresponding QPS fugacity inside the constriction area 
we obtain 

Uc ^ S'coro(Sc/s)exp(-S'coro(Sc/s)), (150) 

where S'core is the QPS core action (|113p for the wire. As 
before, we can express the partition function Z in terms 
of series in powers of the QPS fugacity (cf., eq. (|119p ) 
which is now equal to tJc inside the constriction and to 
y otherwise. In other words, we should now deal with 
two different types of logarithmically interacting QPS. 
Quantum phase slips in the wire are treated exactly as 
before, and we again arrive at the BKT-like QPT at 
/i = /I* w 2. QPS inside the constriction interact log- 
arithmically only in time, thus giving rise to the Schmid- 
likc QPT at /J, = 1. In this case the quasi-lD supercon- 
ducting wire with Mooij-Schon plasmons as elementary 
excitations plays the role of an effective dissipative envi- 
ronment for the constriction. Such QPT was discussed 
by Hekking and Glazman [si?] in the case of a narrow 
superconducting ring interrupted by the Josephson junc- 
tion. In the situation considered here this QPT appears 
to be of little importance since //*>!, i.e. at T ^ 
BKT-like QPT occurs when quantum phase slips inside 
the constriction are still bound in pairs. 

On the other hand, at not very low temperatures and 
provided j/c S> y quantum phase slips inside the con- 
striction may give an additional contribution to the wire 
resistance with the temperature dependence (I149p . In 
order to demonstrate that we modify eq. (|142p by taking 
into account the difference between QPS fugacities in the 
wire and in the constriction. This is trivially handled by 
splitting the space integrals into those over the wire and 
constriction areas. As a result, the linear resistance of 
the system becomes 

RiT) = R.^{T) + RciT), 

where Rw{T) is determined by QPS inside the wire (|145p . 
i.e. Rw{T) cx y^T'^^~^, while the additional contribution 
Rc{T) comes from QPS inside the constriction. This con- 
tribution has the form Rc{T) cx y^T^^"^. Hence, at suf- 
ficiently high temperatures the constriction contribution 
Rc(T) can dominate over Ryj{T), in which case the de- 
pendence (|149p applies. However, at lower temperatures 
Rc{T) becomes irrelevant and R{T) again is determined 
by eq. (|145p . Obviously, the same arguments hold in the 
case of several constrictions or, more generally, for some 
wires with strongly fluctuating cross-section. A similar 
scenario was recently discussed in Ref. [8l| where Gaus- 
sian fluctuations of the wire thickness were considered. 

Finally, let us briefly address the non-superconducting 
phase In our opinion, frequently used over- 

simplified approaches, like TDGL or effective Luttinger 
liquid models, are by far insufficient in order to judge 
whether this disordered phase is actually metallic or in- 
sulating, i.e. whether we are dealing with SMT or SIT 
at r = and /i = /i*. Also RG equations (|138p cannot 
unambiguously resolve this issue. To illustrate this point 
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let us compare the behavior of Josephson chains and ho- 
mogeneous wires. The QPS fugacities y of both systems 
obey RG equations (|138p and, hence, for < fi* these 
fugacities grow upon renormahzation in both cases. In 
the case of Josephson chains it only implies strong fluc- 
tuations of the phases across barriers with no suppres- 
sion of superconductivity inside metallic grains. In this 
case inter-grain tunneling of Cooper pairs is prohibited 
by electron-electron interactions and the system is in the 
insulating (Coulomb blockade) phase. 

On the contrary, large renormalized QPS fugacity y 
in the case of homogeneous wires implies that the QPS 
gas becomes dense, the QPS cores overlap and, hence, 
for fi < fi* and T ^ superconductivity gets eventu- 
ally destroyed everywhere in the wire. In other words, 
in this regime quantum fluctuations drive homogeneous 
ultra-thin superconducting wires normal. This impor- 
tant difference between homogeneous and granular wires 
was emphasized in Ref. [28]. What remains is to figure 
out whether such normal diffusive wires stay metallic or 
turn insulating at low enough T. 

Few years after the work [2^ it was realized that, 
similarly to highly conducting tunnel barriers [13, [H, 
[S^ l disordered metallic wires with large conductances 
also demonstrate the feature of weak Coulomb blockade 
[H, S US S m caused by non-trivial interplay be- 
tween scattering and electron-electron interactions. The 
presence of an (exponentially small) Coulomb gap im- 
plies that disordered wires most likely turn insulating at 
T = 0. However, in generic metallic wires with many 
conducting channels such a state can only be observed at 
exponentially low T, while for any experimentally attain- 
able temperature one expects a metallic behavior with 
a weak trend for the wire resistance to grow with de- 
creasing T due to Coulomb blockade corrections [s^ . In- 
terestingly, such a trend was indeed observed in recent 
experiments [13, [13, [HI and the measured dependence 
R{T) was found to closely follow theoretical predictions 
[Sq l . Further theoretical analysis of the role of electron- 
electron interactions in quasi- ID structures [3, [13, [Ml 
would be highly desirable [9^ . 



VI. EXPERIMENTS ON SUPERCONDUCTING 
NANOWIRES 

A. General considerations 

It follows from our theoretical analysis that both TAPS 
and - in particular- QPS can be observed in sufficiently 
thin superconducting wires, i.e. experiments aimed to 
observe QPS-related phenomena should be performed on 
samples with smallest possible cross section values s. On 
top of that, it is highly desirable to deal with structurally 
(and chemically) homogeneous samples in order to elim- 
inate spatial variations of Tc along the wire and rule out 
the effect of (possibly existing) constrictions and tunnel 
barriers. It has been shown that structural imperfections 



of real quasi-lD structures (non-uniform cross section, 
existence of probes and finite length effects) might effec- 
tively mask the phenomena related to thermal or quan- 
tum fluctuations [o^. 

In addition to the above requirements, in order to pro- 
vide optimal conditions for observation of fluctuation ef- 
fects it is crucial to properly choose the sample material. 
This is obvious from the fact that, for instance, the QPS 
core action (I113P at low enough T can be represented in 
the form 

•Score ^ s/ ppf, 

where pjv = l/c is the wire resistivity in the normal 
state. Hence, it is desirable to select superconducting ma- 
terials with smaller values of Tc (or Ao(0)) and perform 
experiments on sufficiently dirty nanowires with higher 
values oi pn- The latter requirement is in line with the 
well known general observation that fluctuation effects 
are more pronounced in dirtier systems. This require- 
ment also implies that the electron elastic mean free path 
should obey the "dirty limit" condition I < ^ ^ Vl^ 
where ~ fivp/^o is the BCS "clean" coherence length. 
The parameters for various conventional superconductors 
are listed in Table I. Judging from these numbers the 
most suitable materials for experimental investigations 
of QPS effects are those with higher resistivity values, 
i.e. MoGe and a:InO. On the other hand, a certain dis- 
advantage of these materials is that then can be strongly 
inhomogeneous. In this respect, Zr or Ti can be advanta- 
geous for experimental studies of quantum fluctuations. 

Though measuring the temperature dependence of the 
system resistance R(T) in the vicinity of the transition 
point might seem a routine experimental task, in the case 
of nanostructures more care is required. Typically, R{T) 
dependencies are measured in the current-biased regime. 
A standard requirement is to keep the bias current / 
much smaller than the critical (depairing) current Ic 
in order to avoid hysteresis effects due to overheating. 
Additionally, in order to stay within the linear response 
regime, the measuring current / should remain smaller 
than the charactristic scale Iq = ksTc/'^o 0, [13 equal 
to few tens of nA for the majority of materials (see Table 
I) . The characteristic normal state resistance of a typical 
metallic nanowire with diameter ^/s ~ 10 nm and length 
X ~ 1 /im is of order few kO. Hence, in the normal state 
the expected voltage is at least few p, V, which is certainly 
not a problem to measure. However, with decreasing 
temperature only slightly below Tq the resistance R{T) 
drops exponentially and the measured signal quickly re- 
duces to the nV range. Employing room-temperature 
electronics, it is preferable to use ac lock-in technique in 
order to increase the signal-to-noise ratio. This can be 
associated with some hidden problems. One of them is 
that even a tiny fraction of dc component (e.g. from the 
ground loop) adds a parasitic signal ~ dV/dI{T) to the 
"valuable" one i?(T) = V{T)/I. For this reason it is ad- 
visable to decouple the ac current source from the sample 
using a low-noise transformer. 
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Material 






Co 


Sc(0) 


A(0) 




Pjv 




K 


nm 




mT 


me V 


lO" m/s 


X cm 


W 


0.015 


- 


165* 


0.1 


0.002* 


1.8* 


5.7 


Ir 


0.11 


- 


22.5* 


1.6 


0.017* 


1.8* 


5.1 


Hf 


0.13 


- 


19.0* 


1.3 


0.020* 


1.8* 


35 


Ti 


0.40 


- 


6.2* 


5.5 


0.061* 


1.8* 


42.0 


Zr 


0.61 


- 


4.0* 


4.7 


0.093* 


1.8* 


44.0 


Zn 


0.86 


28 


1.8 


5.5 


0.120 


1.82 


5.9 


Al 


1.19 


16 


1.6 


10 


0.175 


2.02 


2.7 


Q:InO 


2.8 


- 


0.330 


- 


0.430* 


- 


3000 


In 


3.41 


65 


0.360 


28 


0.520 


1.74 


8.4 


Sn 


3.73 


50 


0.230 


30.5 


0.557 


1.88 


11.0 


MoGc (film) 


5.5 


720 


0.005 


66 


1.1 




200 


Pb 


7.2 


40 


0.090 


80.3 


1.365 


1.82 


20.7 


Nb 


9.25 


85 


0.030 


206 


1.520 


1.37 


12.5 



Table-I. Material parameters of several superconductors. 
The data are taken from Refs. [H, [H, [H, [oj; for MoGe 
from Refs. gl, H, HI] and for a:InO from Refs. p^.[lOll]. 
The numbers marked with * are estimated values. 

Another problem might originate from the presence 
of rf filters which are mandatory to protect nano-sized 
samples from noisy electromagnetic environments. Very 
often such filters are just RLC circuits shunting the rf 
component to the ground through a capacitor as high as 
few ^F. Such configuration might provide good results 
with high-Ohmic systems (e.g. tunnel structures). In a 
superconducting nanowire measured in 4-probe configu- 
ration each electrode contacting the "body" of the sample 
is typically made of the same material. Hence, the resis- 
tance of these probes also varies from few kfi down to 
zero over the same temperature range within the transi- 
tion. Depending on particular configuration, at a certain 
T the ac sample impedance might become comparable 
to that of the current leads through the ground, caus- 
ing the current re-arrangement throughout the sample. 
Even at rather low measuring frequencies ~ 10 Hz the 
parasitic effect might manifest itself as a non- monotonous 
i?(r) dependence with unusual "bumps" or "foot" at the 
bottom of the transition. Additional complication might 
arize from the non-negligible dependence of the gain of 
the nanovolt pre-amplifier on the total impedance of the 
load. 

Concluding this part, an experimentalist should be ex- 
tremely cautious in designing the measuring set-up for 
unambiguous interpretation of the data on superconduct- 
ing nanowires. 

Structural and geometrical homogeneity of quasi- ID 
samples is the central question in interpretation of ex- 
perimental data related to contribution of superconduct- 
ing fiuctuations. Theoretically one usually assumes that 
(i) the critical temperature Tc of quasi-lD wires under 
consideration remains spatially constant, i.e. it does not 
vary along the wire, (ii) the cross section s does not vary 
along the wire either and (iii) the measuring probes are 
non-invasive. Unfortunately, in realistic nanowires none 
of these conditions is usually well satisfied. To which 
extent can these imperfections be neglected while inter- 
preting the experimental data using models developed 



under the assumptions (i)-(iii)? 

If structural imperfections (e.g. non-uniform chemical 
composition) can alter the critical temperature, then the 
shape of the experimentally observed R(T) dependence is 
determined by the sequence of transitions of various parts 
of the wire with different local critical temperature Tc{x) 
. If the degree of such inhomogeneity is not too strong, 
a step-like R{T) transition may not be observed as the 
variations of the critical temperature are averaged on the 
scale of the coherence length ^ resulting in a relatively 
wide " smooth- looking" i?(T) dependence. 

One might naively expect that working with super- 
conducting samples fabricated of initially pure material 
can help to eliminate the problem. Unfortunately this 
is not the case. First, properties of low-dimensional su- 
perconductors are known depend on the fabrication pro- 
cess details, such as thin film deposition rate, residual 
pressure in the vacuum chamber, material of the sub- 
strate etc. Second, even if to make an effort to keep the 
fabrication parameters constant, it is hard to get rid of 
size-dependent effects. 

It is a well-known experimental fact that the critical 
temperature Tc of thin superconducting films frequently 
differs from one for bulk samples. The same tendency 
is observed in metallic nanostructures. In indium, alu- 
minum and zink T(7 increases wi th d ecreasing character- 
istic dimension [Ifl SI, SI, S ^Q^- On the contrary, 
in lead, niobiu m and MoGe a n opposite tendency is ob- 
served [H H [lOl [lOl [lOi]. No noticeable variations 
of Tc were detected in tin nanowires fl06l [iOtI [lOSl [l09l | . 
The origin of this phenomenon is no t clea r. There exist 
model s pr edicting both suppression jlld l and enhance- 
ment of the mean field critical temperature for low 
dimensional superconductors. One can take the empir- 
ical fact of variation of Tc with effective diameter of a 
superconducting nanowire as granted. Using Al lift-off 
fabricated nanowires as a representative example (Figs. 
|3] and |4]) it has been shown that size-dependent effects 
are extremely important for interpretation of fluctuation 
phenomena in quasi- ID sj^stems [93|. The conclusion of 
a quantitative analysis [93] is rather disappointing: even 
chemically pure nanostructures cannot be considered as 
sufficiently homogeneous as soon as the size dependence 
of the critical temperature Tc comes into play. One 
should study atomically homogeneous systems as single 
crystalline whiskers [l^, [l^ . Unfortunately, modern nan- 
otechnology does not enable growth of high-quality quasi- 
lD single crystalls of arbitrary diameter made of all ma- 
terials of interest. The lithographic processes results in 
much lower quality samples. The only exception are ma- 
terials, such as A/, Zn and In, where the critical temper- 
ature decreases with reduction of the effective dimension 
(e.g., nanowire diameter). Under certain conditions, ex- 
periments on lift-off fabricated nanowires made of these 
materials can be interpreted using fluctuation models de- 
veloped for perfect ID channels. Provided broadening of 
the R{T) dependencies is detected below the bulk critical 
temperature T^""^', the size effects cannot account for the 
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FIG. 3: Empirical dependence of the critical temperature Tc 
on cross section a = s for aluminum nanowires. Line is a 
guide for the eye. Inset: SPM measured distribution of cross 
sections for a typical lift-off fabricated nanowire: length X 
— 10 fim and effective diameter -^s « 75 nm [93|. 
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FIG. 4: Bullets (•) correspond to simulated i?(r) transition 
without taking into consideration the nodes. Diamonds (o) 
represent the results of similar calculations with contribution 
of the node regions. Open circles (o) denote the experimen- 
tally measured i?(T) dependence. Simulations and experi- 
mental data are for the same wire as in the inset from the 
previous figure [osl ]. 



phenomenon, and presumably "real" physics is observed. 



B. Experiments on TAPS 

As we have already discussed in Chapter 4 the TAPS 
mechanism provides non-zero voltage drop across the su- 
perconducting wire at temperatures below Tq. For ref- 
erence purposes let us rewrite again the expression for 
the I — V curve derived in the presence of TAPS at suffi- 



ciently low bias currents /. Combining eqs. ([94)) and ([98]) 
and restoring some fundamental constants (set equal to 
unity in our theory analysis) we obtain 
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where ^(T) is the temperature-dependent coherence 
length of a dirty wire, /cg is the Boltzmann constant. 
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is the the potential barrier for TAPS and C(3) ~ 1.202. 

Very quickly after development of the LAMH theory 
[IB, [l3| two experimental groups reported experiments 
aimed at verification of the TAPS model. At 70-s micro- 
fabrication technique was not much developed. In these 
early studies of ID superconductivity [3 [3 metallic 
whiskers (Fig. [5]) with characteristic diameter ~ 1 ^m 
and lengths up to 1 mm were utilized. Growth and basic 
properties of these highly anisotr opic objects have been 
widely described in the literature |ll2l | . In the particular 
case of tin whiskers " s queez e" method was typically ap- 
plied for their growing [ll3j |. Then the crystals were lit- 
erally hand picked from the ingots and positioned on sub- 
strates. The ele ctro des were made either by conducting 
paste or epoxy fTTi|, soldered by Wood's metal [H, [Hi 
or Sn- Pb al loy 19 ], or squeezed by a soft metal (e.g. in- 
dium) (ll5l |. An obvious disadvantage of the procedure 
is inevitable damage of the crystal within the locus of 
the electric probes. An attempt to overcome these dif- 
ficulties was made by combining planarization technique 
with electron beam lithography [ll6l Ill7j | (Figs, [5] and 
[S]), The technique appears promising, but extremely time 
consuming and resulting in a rather low yield of "work- 
ing" samples. 

An example of experimental R{T) dependence mea- 
sured on tin whisker is shown in Fig, [T] One observes 
a very goo d q uantitative agreement between the exper- 
iment 0, [ig and the TAPS model [H, [T^. The su- 
perconducting transition R{T) is very steep: resistance 
of the sample drops five orders of magnitude within the 
temperature range 5T ~ 1 mK below Tc- This is a conse- 
quence of rather large effective diameter values ^/s ~ 0,5 
/im and extremely high homogeneity of single crystals. In 
the case of thinner single-crystalline structures the width 
of the i?(T)-curve would be larger, cf, eq, ([M]) . Unfortu- 
nately, the dimensions of whiskers made of superconduct- 
ing materials do not vary much. Manual manipulation of 
sub-1 ^.m is extremely time consuming and results in a 
low yield of suitable samples, while the dimensions of 
whiskers are still too large to use scanning probe (SPM) 
technique well developed for manipulation of nano-sized 
objects, such as carbon nanotubes. 

An alternative to whiskers is superconducting mi- 
crocylinders encapsulated in a dielectric substrate (e,g , 
fab ricated u sing Taylor-Ulitovski method jllSl . 
[ml [12^ , The quality of these wires is not 
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FIG. 5: SEM image of a typical tin whisker on the surface 
of spin-on-glass III6II . Inset shows the magnified view of the 
crystal surface [llTl ]. 




FIG. 6: SEM image of a whisker-based microstructure on 
Si/SiO substrate covered with Spin- On- Glass. Slightly bent 
bright horizontal line is the tin whisker. Electrodes are made 
on top of the crystal by e-beam lithography followed by evap- 
oration of copper contacts [lid ] . 



as high as that of atomically-piire whiskers. However, 
an important advantage of this approach is the possi- 
bility to fabricate very long samples. Fabrication of a 
wire with metal core of about few /im and length up 
to a km (!) is a routine work, while fine adjustment of 
parameters enables fabricat ion of sub-100 nm filaments 
with length up to few cm [HI [HI [l2i|. A remark- 
able feature of wires obtained by this method is that 
typically they are single crystals. A number of super- 
conducting microcylinders (Sn and In ) fabricat ed using 
Taylor-Ulitovski method were studied jl26l . Il27l | . Bigger 
length values of the filaments (X ~ 1 cm ) as compared to 
those for previously studied whiskers {X ~ 0.5 mm ) en- 
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FIG. 7: R(T) dependence for a tin whisker [l^. Solid symbols 
are the experimental data, dashed line represents the fit to the 
TAPS model [3, [13]. 



abled observation of interesting features related to inter- 
action of current-induced phase slips. The cross section 
values for microcylinders did not differ much from those 
for whiskers, while the homogeneity was clearly worse: 
No qualitatively new phenomena related to fluctuation- 
governed ID superconductivity has been observed. 

There were several reports on experimental studies 
of ID superconductivity in lift-off-fabricated supercon- 
ducting nanostructures (wires and loops) with the line- 
width much smaller tha n of whiskers or microcylinders 
[m [ili, [133, [ml [m . The advantage of the method 
is the ability to fabricate samples with a complicated 
shape in a reproducible way. Typically, the metal is de- 
posited through PMMA mask using thermal or e-beam 
evaporation, or sputtering. The quality of the struc- 
tures is far from perfect: the majority of superconducting 
thin-film structures fabricated at room temperatures are 
polycrystalline with the grain size of about few tens of 
nm. Deposition of metal (in particular low melting ones, 
such as tin or indium) on a cryogenically cooled sub- 
strate might reduce the grain size. In contrast to whiskers 
or microcylinders, the lift-off-fabricated superconducting 
nanostructures studied so far were all in the dirty limit 
I <^ S,. Attempts to quantitatively describe the shape 
of superconducting transition of these quasi-lD systems 
using the TAPS model failed: Experimental curves for 
R{T) were always significantly broader than theoretical 
predictions. Nevertheless, it was believed that with a 
certain adjustment of fit parameters (e. g., re duction of 
the effective size of the phase slip center [l3ll |) a reason- 
able agreement between the experiments and the TAPS 
model predictions could be achieved. Later it was shown 
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that inevitable inhomogeneity of lift-off- fabricated nanos- 
tructures, the presence of node regions and finite size 
effects can dramatically broaden the experimentally ob- 
served dependencies R(T) making any comparison with 
the TAPS theory inconclusive or even impossible (Fig. 

Very thin superconducting nanowires have been stud- 
ied using the template decoration technique: Deposition 
of a metal film on to p of a suspended ins ulating molecule 
used as a template [IM IHl EM Ell (Fig- ©■ Only 
pseudo-four probe configuration is available with the de- 
posited metal film from both sides of the trench serving 
as electrodes. The approach enabled fabrication of su- 
perconducting nanowires (MoGe, Nb and a:InO) with 
effective diameters down to ~3 nm and lengths up to ~1 
/j,m. TEM study revealed amorphous structure for MoGe 
and Q;:InO. 

The R{T) dependencies in a:InO appeared to be sig- 
nificantly wider than it is predicted by the TAPS model, 
though no claims abo ut the impact of the QPS mecha- 
nism were made jlOOl |. Possibly, broadening of experi- 
mental curves R{T) in arInO can be associated with in- 
evitable inhomogeneity of the samples introduced during 
their fabrication. It is well known that physical proper- 
ties of a'.lnO low dimensional systems strongly vary by 
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FIG. 9: R{T) obtained from a:InO nanowire with width 100 
nm together with the data from a similarly prepared 500 fim 
wide film. Inse t: A scanning electron micrograph (SEM) of a 
typical device [lOOl ]. 



annealing 

Interpretation of some experiments on MoGe 
nanowires remains under debates. In earlier reports 
[12, El] very broad R{T) dependencies were observed and 
associated wit h quant um phase slips. However, in later 
works [3, [q^, mi llOSf broadening of a superconducting 
transition in similar MoGe structures was interpreted 
within the TAPS model with no QPS contribu tion (Fig . 
[TU)) . On the other hand, later it was argued [l37l Il38j | 
that fits were actually produced outside the applicability 
range of the LAMH theory (1^51) . M oreo ver, for the 
experimental parameters [H, [90, [Sll) llOS] the strong 
inequalities (|95|) are not satisfied, i.e. this theory has 
almost no applicability range at all, see also Fig. [TT] 
Thus, employing the LAMH model in order to fit the 
resistance curves R(T) for these samples is problematic. 
In addition, even if one performs such fits one is bound 
to use fit values for the electron mean free path much 
larger than the wire diameter (90j which is rather 
unrealistic for such structures (Fig. [10]). At this point 
let us recall that estimates for the co herence le ngth 
C(0) i^oiy^^ both in MoGe thin films [Hi, EH and 
wires ,:32,] typically yield values ~ 5 to 7 nm which 
translates into mean free path values / of^only few nm, 
i.e. much shorter than used in the fits 90]. We will 
return to possible interpretation of these experiments 
below in connection with QPS effects. 

Transmission electron microscopy (TEM) study of Nb 
nanowires, fabricated using template decoration method, 
revealed polycrystallinc structure with an average grain 
size about few nm. In order to fit the data obtained in 
these Nb structures it has been proposed [i^ to represent 
the measured wire conductance 1/ R(T) as a sum of con- 
tributions from normal electrons above the gap and from 
TAPS, 1/R{T) = 1/Rn + 1/Rtaps- The experimental 
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FIG. 11: Fits of the data points for the sample No. 6 of 
Ref. [90I ] to the LAMH theory. The fit parameters are the 
coherence length, ^ — 6.7 nm, and the critical temperature, 
Tc = 4.7 K, marked by an arrow. The LAMH theor y can only 
be appUed within the strip 4.23 K < T < 4.55 K [HI, Mm 
(shown in grey) which is clearly outside the region where R 
strongly depends on temperature. 



FIG. 10: Resistance versus temperature plots for "short" 
MoGe nanowires: insulating (top) and superconducting (bot- 
tom) [13]. Solid lines are fits to the TAPS model. The best fit 
values of the coherence lengths ^(0) are 70.0, 19.0, 11.5, 9.4, 
5.6, and 6.7 nm, and lengths L are 177, 43, 63, 93, 187, 99 
nm for samples 1-6, respectively. Double-step shape of -R(r) 
transitions in the top and bottom plots comes from the su- 
perconducting transition of the contact regions contributing 
to the 2-probe measurement configuration. 



curves R{T) could be reasonably well fitted [lOl [l0| to 
this phenomenological model (Fig. [T^]). Still, the abil- 
ity to detect TAPS contribution in Nb is rather unex- 
pected since the energy gap Aq is very large in this ma- 
terial (see Table I). It is also quite surprising that the 
inhomogeneities in granular Nb nanowires do not seem 
to broaden the R{T) curves, which can be fitted by the 
TAPS model developed for homogeneous quasi-lD sys- 
tems. 



C. Experimental observations of QPS effects in 
superconducting nanowires 

Let us now turn to the experiments detecting QPS ef- 
fects in superconducting nanowires. Mooij and coworkers 
[l4ll | discussed the possibility to observe quantum phase 
slips experimentally and attempted to do so as early as 
in 1987. According to our theory, the crossover between 
TAPS and QPS regimes can be expected at temperatures 
T ~ Ao(r). For superconductors with typical material 




FIG. 12: Temperature dependence of the resistance of super- 
conducting Nb nanowires obtained by template decoration 
method using carbon nanotube as a substrate. The follow- 
ing fit parameters for the samples Nbl, Nb2, Nb3, Nb5, and 
Nb6 are used: Transition temperatures (K) are T'^"'<= = 5.8^ 
5.6, 2.7, 2.5, and 1.9 K; lengths (nm) are X = 137, 120, 172, 
177 and 113; normal state resistance (f2) are _Rjv — 470, 650, 
1600, 4250 and 9500; coherence lengths (nm) are ^(0)= 8.5, 
8.1, 18, 16, and 16.5, respectively. Solid lines show fits to the 
phenomenological TAPS model ^4^. The dashed and dotted 
lines are some theoretical curves that include QPS contribu- 
tions. 
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FIG. 13: Resistive transitions for a set of MoGe wires with 
diffe rent width w upon a single film of thickness t = 5 nm 

[i3i. 



parameters and Tc > IK (see Table I) this condition im- 
plies that QPS effects may become important already at 
Tq — T > 100 mK. As we already discussed, in order to be 
able to observe QPS it is necessary to fabricate nanowires 
with effective diameters in the 10 nm range, -y/s ~ 10 
nm. The wires should be sufficiently uniform and ho- 
mogeneous in order not to override fluctuation effects by 
trivial broadening of R(T) dependencies due to wire im- 
perfections and inhomogeneities [osl ]. Thus, proper fab- 
rication technology is vitally important for experimental 
studies of QPS in quasi-lD superconductors. 

Perhaps the first experimental indication of the ef- 
fect of quantum fluctuations was obtained in amorphous 
MoGe nano wires with effective diameter down to ^/s ~ 
30 nm [139| |. The samples in 4-probe configuration with 
length X up to 20 /im were fabri cated using e-beam 
lithography. Though the paper [l3!^ was mainly focused 
on the effect of disorder in low-dimensional superconduc- 
tors, it was clearly stated that for the narrowest samples 
significantly broader curves R{T) were observed than it 
is predicted by the TAPS model (see Fig. [T3]). The effect 
of quantum fluctuations was pointed out as a possible 
reason for this disagreement. 

Detailed experimental studies of transport properties 
of superconducting nanowires have been carried out by 
Giordano and co-workers [13, [13, IH, H [HI- In these 
experiments In and In-Pb wires with triangular cross- 
section were fabricated using step-edge lithographic tech- 
nique |l43l . Il44l |. This method utilizes the shadow ef- 
fect produced by extremely shallow steps in substrates 
and the corresponding tilted ion beam milling of the de- 
posited material (Fig. [T4|) . Polycrystalline wires were 
fabricated with effective diameters ^/s in the range 40 - 
100 nm with the grain size from 10 to 20 nm. The unifor- 
mity of the samples was controlled with SEM and with 
the accuracy of ~ 10 nm. 

For wider structures a reasonable agreement with the 




FIG. 15: Resistance (normalized by its normal state value) 
as a function of temperature for three In wires; the sample 
diameters were 41 nm (•), 50.5 nm (-I-) and 72 nm (o). The 
solid curves are fits to the TAPS model, while the dashed 
curves indicate fits to the phenomenological model [2C|. The 
sample lengths are 80, 150 and 150 nm, and the normal state 
resistance values are 5.7, 7.1 and 1.2 kfl, respectively [20l |. 



TAPS model was observed. On the other hand, for thin- 
ner wires with ^/s < 50 nm clear deviations from TAPS 
predictions were demonstrated (see Fig. [T5t . The dis- 
crepancy was interpreted as a manifestation of quantum 
phase slippage. In order to explain their observations 
Giordano and co-workers proposed phenomenological de- 
scription based on the Caldeira-Leggett model for macro- 
scopic quantum tunneling with dissipation [25|. Though 
qualitative agreement with this simple phenomenological 
model has been obtained (Fig. llSp . quantitative inter- 
pretation of the data [20| is problematic due to poor uni- 
forrnity of the samples: ± 10 nm for wires with y/s < 50 
nm [20|. 

Recently, Pai et al. 8l| performed a fit of the same 
experimental data [20j to the power-law dependence 
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FIG. 16: Cross-sectional view of the stencil structure. Metal- 
lic films were evaporated th roug h the shadow mask structure 
fabricated on the substrate [l03l ]. 



R{T) oc T^^ ^ derived from the result (fT45)) assumme; 
Gaussian fluctuations of the wire thickness. Unfortu- 
nately the temperature dependence of the QPS fugac- 
ity (|136p was not taken into account in Ref. [8l|. This 
temperature dependence (predominantly determined by 

1 /2 

that of the gap, Score oc Aq' (r)) should, however, still 
be significant at temperatures not very far from Tc and 
will alter the fits. 

It appears that strong granularity of the wires was 
most likely a very important factor in the experiments 
by Giordano and co-workers. For instance, estimating 
S'core for the thinnest wires |4S] with diameters 16 and 
25 nm, for the experimental parameters we obtain re- 
spectively S'core ~ 280^1 and 700A. For uniform wires 
such huge values of the QPS action would totally pro- 
hibit any signature of quantum phase slips. Since QPS 
effects were very clearly observed, it appears inevitable 
that these samples contained constrictions with diame- 
ters ~ 3 10 times smaller than the average thickness 
values presented in [4^. 

Another set of experiments was performed by Dynes 
and co-workers |l03l . Il45l Il46l | . Suspended stencil tech- 
nique has been developed (Fig. [11] ) enabling fabrica- 
tion of qucnch-condensed granular nanowires with cross- 
section area s down to 15 nm^ and length X ranging 
from 1 to 2 /xm from various materials including several 
superconductors, such as Pb, Sn, Pb-Bi. The samples 
edge roughness was claimed to be about 3 nm. A re- 
markable feature of the method is the ability to vary the 
wire thickness t and, hence, its cross section s at a con- 
stant width w in situ inside the cryostat in-between the 
sessions of truly 4-probe R{T) measurements. 

The experiments (l03j | clearly indicated systematic de- 
viations of the experimental data points R{T) from the 
TAPS model predictions (Fig. [T7]). This discrepancy in- 
creases as the wires become narrower. The width of the 
superconducting transition was found to scale with the 
normal state resistance Rn (Fig. fT8|) . It should be noted 
that in lead nanowires as narrow as 15 nm and as thin as 
10 nm no low temperature resistance tails were observed. 
Instead, a less dramatic but systematic broadening of the 
superc ondu c ting transition beyond the TAPS limit was 
noted [l03l Il45| . On the other hand, long resistance 
tails were always present in tin structures fabricated and 
measured using similar technique [l46| |. Very probably. 
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FIG. 17: R{T) dependencies for the same Pb nanowire with 
width w — 22 nm and various thicknesses. Solid lines at the 
top of the transitions are fits to the Aslamazov - Larkin the- 
ory. Dashed lines are fits to the TAPS model, which clearly 
yields steeper R{T) dependencies compared to the experimen- 
tal data for the thinnest wires (103.] . 
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FIG. 18: Width of the transition (defined as a temperature 
inter val between 20% and 80% of the normal state resistance 
[103| ) as a function of the normal state resistance Rn- 



the discrepancy comes from the difference in the heights 
of the potential barrier SF between these two materials 
(see Table I) : Lead wires should be significantly narrower 
than tin ones in order to obtain a similar magnitude of 

fluctuation effects. 

In add ition, we point out that the experiments |l03l . 
Il45l Il46l | demonstrated clear superonducting transitions 
reaching experimental "zero" in nanowires with normal 
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state resistance Rn ^ 10 kfl. This observation is not in 
line with the conjecture [H, HI] that superconductivity 
in MoGe nanowires could only be possible for wire re- 
sistances below the quantum resistance unit Rq ~ 6.45 

kn. 

A clear manifestation of a crucial role of QPS effects 
in superconducting nanowires was provided in the ex- 
periments by Tinkham and co-workers [3^ . l43j . These 
authors developed a novel technique which allowed them 
to fabricate sufficiently uniform superconducting wires 
considerably thinner than 10 nm with lengths ranging 
between ~ 100 nm and 1 fim. This was achieved by sput- 
tering a superconducting alloy of amorphous Mo7QGe2i 
over a free-standing carbon nanotube or bundle of tubes 
laid down over a narrow and deep slit etched in the sub- 
strate (Fig. E). T hree out of eight samples studied in the 
experiments [32| demonstrated no sign of superconduc- 
tivity even well below the bulk critical temperature Tc- 
Furthermore, in the low temperature limit the resistance 
of these samples was found to show a slight upturn with 
decreasing T. In view of that one can conjecture that 
these samples may actually become insulating at T — > 0. 
The resistance of other five samples [HI decreased with 
decreasing T. Also for these five samples no sharp super- 
conducting phase transition was observed. 

All three non-superconducting wires (il,i2 and i3) [111 
had the normal state resistance above the quantum unit 
Rq, while Rn for the remaining five "superconducting" 
samples was lower than Rq. This observation allowed 
the authors |3l| to suggest that a dramatic difference in 
the behavior of these two groups of samples (otherwise 
having similar parameters) can be due to the dissipa- 
tive quantum phase transition (QPT) [2^ [sl] (see also 
Chapter 5) analogous to Schmid phase t rans ition (2ll.l27| 
observed earlier in Josephson junctions Jl47'|. Already at 
this stage we would like to emphasize that an important 
pre-requisite for this QPT is the existence of a source for 
linear Ohmic dissipation at low energies which typically 
requires the presence of some normal shunt resistance. It 
appears that no such condition was fulfilled in the exper- 
iments (32} . Hence, the interpretation of the data [33 1 in 
terms of a dissipative QPT is problematic. 

And indeed, this interpretation was not confirmed in 
the later experiments of the same group [43i] who ob- 
served superconducting behavior in samples with normal 
resistances as high as 40 kfi ^ Rq. The authors [i^l con- 
cluded that "the relevant parameter controlling the su- 
perconducting transition is not the ratio of Rq/R^, but 
appears to be resistance per unit length, or equivalently, 
the cross-sectional area of a wire" (Fig fTOl) . 

This conclusion clearly favors interpretation of the 
data [si, SI either in terms of a BKT-like QPT [H, ^ 
at « 2 or just as a sharp crossover between the regimes 
of vanishingly small and sufficiently high QPS rates 7qps 
which will correspond respectively to superconducting 
and normal behavior of the nanowires. In both cases 
the crucial parameter is the wire cross section s. Both 
these QPT and crossover are expected to occur for wire 
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FIG. 19: Resistance at 1.5 K normalized to normal state re- 
sistance as a function of L/Rn- (a) Linear plot. The dotted 
line is a guide to the eye. (b) Semilog plot with an exponential 
fit. Slope of the fitted line is 0.39 kSl/ nm [13]. 



diameters ~ 10 nm. 

Let us estimate the QPS core action S'coro defined by 
eq. (|113p for the eight samples studied in the experiments 
[3^]. With the density of states Nq = 1.86 x 10^^ s/m^, 
the superconducting critical temperature Tc ~ 5.5 K and 
the measured resistivity p = 1.8 iiQ/m we obtain the co- 
herence length ^ ~ 7 nm in agreement with the estimate 
[3^ . Our estimates for the action S'coro are summarized 
in the following Table II: 



sample 


R/X, kn/nm 


'S'corc 


if 


0.122 


7.8 A 


i2 


0.110 


8.7 A 


iS 


0.079 


12.7 A 


si 


0.038 


25.1 A 


s2 


0.028 


33.7 A 


s3 


0.039 


22.6 A 


ssl 


0.054 


15.4 A 


ss2 


0.044 


19.6 A 
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We observe that ior A ^ 1 the QPS rate jqps 
exp(— S'coro) is expected to be much higher in the normal 
samples il, 12 and 13 than in five remaining wires which 
demonstrated the superconducting behavior in the low 
temperature limit. This observation is consistent with 
the interpretation in terms of the crossover between the 
regimes with low and high QPS rates. 

Further support for this Interpretation comes from the 
data [1^ obtained for more than 20 different nanowires. 
The resistance of these wires measured at T = 1.5 K 
Is presented in Fig. [12] versus the inverse normal resis- 
tance per unit length X/Rm oc s. In Fig. [12] (a) one 
observes a sharp crossover between normal and supercon- 
ducting behavior at wire diameters ^/s ^ 10 nm. This 
crossover could be Interpreted as an indication to the 
QPS-bindlng-unbindlng QPT [H, H^] controlled by the 
parameter fi cx ^/s. Note, however, that all wires studied 
in [4^ were quite short. Hence, this QPT should in- 
evitably be significantly broadened by finite size effects. 

Re-plotting the same data on a semilog scale (Fig. 
[TOl (b)) we Indeed observe a rather broad distribution 
of measured resistances which - despite some scatter - 
can be fitted to the linear dependence on the wire cross 
section s. This fit is highly suggestive of the crossover 
scenario although it cannot yet rule out the (broadened 
by size effects) QPT either. 

In order to finally discriminate between these two op- 
tions it is necessary to analyze the temperature depen- 
dence of the resistance R{T). According to our theory, in 
the linear regime the (caused by QPS) wire resistance is 
defined by eq. (|145p . The dependence i? cx comes from 
pairs of QPS events where the fugaclty (|136p depends on 

1 /2 

temperature via Score oc Ag (T). An additional (weak) 
power-law dependence cx T^t^~^ enters because of inter- 
actions between quantum phase slips in space-time. In 
order to fit the data [43[ we will ignore this power-law 
dependence and use the simplified formula 

R{T) = 6i^^^^|^2i£^exp(-25eorc), (151) 

with S'corc defined in pi3p and 61 being an unimportant 
constant factor. As this formula can only be applied to 
samples with vanis hing low T resistances R{T ^ 0) — > 0, 
we select the data [4J] for the samples 3 to 8 seemingly 
demonstrating such a behavior. The results of the fits to 
eq. (|15ip are presented in Fig. [201 We observe a very 
good agreement between our QPS theory and experiment 
which was obtained using the same value ^ ~ 0.7 for all 
six samples. 

We also note that the fits remain essentially unchanged 
if we take into account an additional power law depen- 
dence cx T^t^~^ (with /i estimated from the wire param- 
eters). This observation proves that Inter-QPS interac- 
tion is Indeed insignificant for the interpretation of the 
data [isj. Thus, similarly to Lau et al. [i^, we con- 
clude that the temperature dependence of the resistance 
of their ultra-thin MoGe wires is determined by QPS ef- 
fects and the observed sharp transition between normal 




FIG. 20; Superconducting transitions of "long" MoGe 
nanowires on top of insulating carbon nanotube used as the 
substrate [s^l, Double-step shape of the R{T) depen- 

dences comes from the superconducting transition of the con- 
tact regions contributing to the 2-probe measurement config- 
uration. The samples' normal state resistances and lengths 
are 1: 14.8 kQ, 135 nm; 2: 10.7 kfi, 135 nm; 3: 47 kfi, 745 
nm; 4: 17.3 kO, 310 nm; 5: 32 kQ, 730 nm; 6: 40 kQ, 1050 
nm; 7: 10 kQ, 310 nm; 8: 4.5 kSl, 165 nm. Symbols stand for 
calculations using eq. 11511 with the single numerical coeffi- 
cient A— 0.7. The critical temperature Tc and the dirty-limit 
coherence length ^(0) used as fitting parameters for samples 
3-8 are 3: 5.0 K, 8 nm; 4: 6.4 K, 8.5 nm; 5: 4.6 K, 8.9 nm; 
6: 4.8 K, 8.9 nm; 7: 5.6 K, 11.9 nm; 8: 4.8 K, 8.5 nm. Data 
for samples 1 and 2 cannot be fitted by any reasonable set of 
parameters. 



and superconducting behavior is most likely a thickness- 
governed crossover between the regimes of respectively 
large and small QPS rates Jqps ■ 

A bulk of recent experimental data [o^, [9l|, Il05j | accu- 
mulated on MoGe nanowires fabricated using the molec- 
ular template method was recently reviewed in detail by 
Bezryadin [4^ . The main observations can be summa- 
rized as follows: (i) "shorter" nanowires {X < 200 nm) 
demonstrate either "weakly insulating" behavior with 
clear features of weak Coulomb blockade [H, HI, 113, , 
or relatively steep superconducting transition R{T) with 
virtually no samples showing an intermediate regime, (ii) 
" longer" samples (200 nm ^ X <1 /im) typically showed 
the behavior which - similarly to the data [43] - can 
reasonably well be interpreted in terms of a crossover 
between the regimes of small and large QPS rates (cor- 
responding to respectively superconducting and normal 
behavior). 

The R(T) curves of longer wires in the regime (ii) 
showed a decrease of the resistance with cooling no mat- 
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ter whether their normal state resistance Rn was smaller 
or bigger than Rg. The crossover between normal and 
superconducting behavior of the samples was controlled 
by the wire cross section s or, equivalently, by the ratio 
Rn/X. As in Ref. the overall picture is consistent 
with the QPS scenario [30l |. 

Turnin g to the regime of shorter wires (i), it was ar- 
gued 0, HH that the transition between superconduct- 
ing and normal behavior for these samples is most likely 
controlled by the total normal wire resistance Rn rather 
than by the resistance per unit length Rn/X, see Fig. 
3 of Ref. [ilj. This observation implies that a signifi- 
cant fraction of short wires with lengths X < 200 nm 
showed superconducting-like behavior even though the 
corresponding estimated values of the QPS core action 
could be as high as ^core ~ 1 (in which case for longer 
wires one expects strong proliferation of QPS leading 
to total destruction of superconductivity). One should, 
therefore, explain why these samples stay superconduct- 
ing instead of being normal or even insulating in the limit 
of low enough T. 

An empirical fact that short samples [3, [9l| stayed 
normal for Rn > Rq and turned superconducting as soon 
BiS Rn ^ Rq allowed the authors to assume the presence 
of some kind of QPT that could, for instance, be simi- 
lar to Schmid dissipative phase transition in Joscphson 
junctions [U [l^, [H, [6l|. Thus, in some sense the au- 



thors [9l| revived an earlier idea S^] , now in application 
to short w ires only. Supporting this idea Meidan et al. 
[l37l Il38j | suggested a phenomenological scheme essen- 
tially employing the well known RG analysis developed 
for resistively shunted Joscphson junctions l2ll. 27l. 54|. 



Unfortunately, as was also noticed in Refs. [1371 Il38l |. 
the experimental parameters correspond to large fugac- 
ity values y_~ 1, i.e . the RG approach ^53] was employed 
in Ref. |137l Il38l | well beyond its applicability range. 
Another - even more significant problem - is that in the 
structures [3 IsH no source for linear Ohmic dissipation 
seems to exist at low energies. Meidan et al. phenomeno- 
logically argued that such dissipation could occur due to 
proliferation of QPS inside the wire (which would make 
the wire to a large extent normal). However, even in- 
trinsically normal wires in contact with superconducting 
electrodes develop a proximity-induced g ap in its energy 
spectrum [HI [mI [lH [HI [HI [l5|| of order Thou- 
less energy exh = D/X'^. Hence, Ohmic dissipation can 
hardly occur in such SNS structures at energies < eTh- 
Since for short samples [H, [9l| with X ~ 40 200 nm 
the energy eTh can easily reach few K, it is hard to expect 
that a dissipative QPT can occur in such systems. 

The same arguments suggest that the proximity effect 
may actually play an important role in the interpreta- 
tion of the experimental data [3, [9l| for short samples. 
Even if the wire is intrinsically normal, supercurrent can 
flow through such a weak link due to the proximity effect. 
Accordingly, superconductivity in the short samples with 
Rn Rq can simply be associated with the onset of dc 
Joscphson current in an SNS junction. The correspond- 



ing scenario is as follows. 

In the interesting limit X £ t he Joscphson energy 
Ej{T) of an SNS junction [Hi, [ill [Hll is exponentially 
small at high temperatures T ^ exh- Hence, at such T 
we have Ej{T) <C T and the weak Joscphson current 
is fully suppressed by thermal fluctuations. Upon de- 
creasing temperature one eventually reaches the regime 
T < eTh, in which case Ej becomes [156] 



Ej ~ erhRq/ Rn- 



(152) 



Thus, for i?7v < i?g at T < eTh we have Ej > T, i.e. 
thermal fluctuations are negligible and the SNS junction 
becomes superconducting. In the vicinity of T ^ eTh this 
argument also allows to predict exponential decrease of 
the wire resistance R{T) cc e'}q){eThRq/ RnT) in a quali- 
tative agreement with exponential dependencies observed 
in di, [qO]. On the other hand, for Rn > i?, at T ~ eTh 
we still have Ej < T, i.e. supercurrent is disrupted by 
thermal fluctuations down to lower T and signiflcantly 
broader R{T) curves are expected. 

At lower T thermal fluctuations become unimportant 
but, on the other hand, quantum fluctuations take over 
destroying the Joscphson current for Ec ^ Ej, where 
Ec — e^/2C is an effective charging energy of an SNS 
junction. Thus, in the low temperature limit one would 
expect to see the crossover at Ej ~ Ec or, equivalently, 
at 



Rq/RN ^ Ec/eT\i- 



(153) 



Now, similarly to the case of Joscphson junctions [21| . 
one can demonstrate that the capacitance C is given 
by the expression C = Cg -\- SC, where Cg is the junc- 
tion geometric capacitance and SC is the renormaliza- 
tion term SC h/RN£Th- Since in our case geometric 
capacitance is most likely very small Cg <C SC, one has 
Ec ^ Rn^tii/ Rq- Substituting this estimate into eq. 
(|153p . one immediately obtains the crossover condition 

Rn ^ Rq 

in agreement with experimental findings. We believe 
that these simple arguments are sufficient to understand 
a superconducting-to-normal crossover observed in short 
wires [13, [13, l9ll |. These arguments also demonstrate 
that for such samples Joscphson physics (both classical 
and quantum) appears to be more important than that 
of TAPS and QPS. Hence, although short wires can also 
demonstrate interesting phenomena, longer wires appear 
to be more suitable for experimental investigations of 
QPS effects. 

Now let us turn to other experiments where QPS in 
superconducting nanowires have been observed. 

The original m etho d of forcing t he m olten alloy into 
a porous media [l57j or capillary |l58l | enables fabri- 
cation of metallic wires with diameters >1 /Ltm. At 
these early stages of research no studies of the shape 
of superconducting transition were performed in such 
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FIG. 21: R(T) dependence for four Pb nanowires (A) - (D) 
with diameters 40, 55, 55 and 70 nm, c orre spondingly. The 
transition is broader for narrower wires II6III. 
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FIG. 22: R{T) curves for 20, 40, 60, 70, and 100 nm wide and 
6 /xm long Sn nanowire arrays containing respectively 18, 1, 
8, 15 and 53 wires in the bundle. The solid lines for 20, 40 
and 60 nm wires are the results based on the TAPS model 
near Tc and QPS model below Tc [TOTj . 



structures. Chemical or electrochemical methods devel- 
oped later allowed filling of very narrow pores with a 
wide variety of ma, terial s inc l uding superconduc tors: Pb 
15^,'T6^,|lM[i6|, Sn [lOSl . flM llM [l6^ . IT63 | and Zn 
102, 165, il66l |. The most widely used media are track- 
etched membranes and alumina nanoporous films |l09l |. 
The diameter of pores in these materials can be rather 
uniform ranging from ^15 nm to few /im. One side of 
the membrane is coated with a thick metal film serving 
as a cathode. The "sandwich" is immersed into a corre- 
sponding electrolyte. Standard three-terminal reduction 
of ions in the solution leads to growth of nanowires from 
the cathode. Typically the nanowires are polycrystalline. 
However, with properly adjusted electrodeposition con- 
ditions nanowires can be made crystalline [107|. The 



method enables fabrication of large amounts of almost 
vertically aligned nanowires of desired morphology and 
diameters. Electric contacts to a bundle of nanowires 
can be made easily connecting the thick bottom cath- 
ode (from which the filaments start to grow) and the top 
metal plate on top of the nanopore membrane. The num- 
ber of nanowires contacting metallic electrodes depends 
on squeezing pressure of t he t op electrode and the filling 
factor of the membrane jl65j |. For unambiguous inter- 
pretation of the data it is desirable to perform electric 
measurements on a single nanowire. If the density of the 
nanopores is sufficiently high, a possible solution is to dis- 
solve the hosting membrane, isolate a single nanowire and 
fabric ate contac ts using conventional lithographic tech- 
nique [l09l Il67j | . A nanoindentation-based method has 
been proposed enabling electric measurements on a single 
nanowire inside the hosting membrane jl68j |. 

Experiments with Zn nanowires [I0i,|i6i,[l6^ lead to 
observations of unexpected "antiproximity" effect. The 
shape of the supeconducting transition R{T) of a bun- 



dle of Zn nanowires with length up to few /im depends 
on the material of the electrodes. It was found that su- 
perconductivity is completely or partially suppressed in 
samples with superconducting electrodes, while clear su- 
perconducting transition has been observed in systems 
with normal electrodes. The origin of the phenomenon 
is not yet clear, while a model has been proposed [169]. 

Experiments with individual Pb nanowires grown us- 
ing nanopore electrochemical method show clear super- 
conducting transitions, which get broadened with reduc- 
tion of the sample diameter (Fig. [21]) . Unfortunately no 
theory fits was provided by the authors, while the linear 
scale of the reported R{T) dependenc ies c omplicates an 
independent quantitative comparison |l6l| . 

Experiments on bundles of several Sn nanowires show 
clear broadening of the R{T) transitions in the narrowest 
samp les (Fig. [22|) , which was associated with QPS effects 
[lOTj l . As the number of wires in the measured bundle is 
not known precisely, quantitative comparison with theo- 
retical predictions appears complicated. In addition, the 
experiments were presumably performed outside the lin- 
ear regime: Only rather high dc excitation current densi- 
ties jbias 10"* A/cm^ were used being just shghtly below 
the experimentally measured critical current density at 
low T, jc{T < Tc) ~ 10^ A/cm^. 

An advanced method of nanowire fabrication uses 
MBE grown InP layer on a cleaved InxGai-.xAs/InP 
substrate as a support for thermally evaporated metal 
(Fig. [23]). The applicabili ty of the approach was demon- 
strated on AuPd and Al |l70 lll7ll |. The method allows 
fabrication of very long (up to 100 /im) superconducting 
{Al) nanowires with effective diameters down to 7 nm. 
The approach enables pseudo-4-terminal measurements 
(Fig. [M)) . The observed broad R{T) dependencies were 
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(a) (b) (c) (d) (e) 



FIG. 23: Schematics of the sample fabrication, (a) The sample is cleaved in small strips which are cut in half and glued together 
with the two (001) plane facing each other, (b) PMMA is spun on the (110) crystallographic plane of the two pieces and a 
pattern is written using standard EBL and then developed, (c) thermal evaporation is used to deposit a film of Ti/Au (the 
portion of the film deposited on the PMMA is not shown) and then the two halves are separated, (d) after lift-off and oxygen 
plasma etching, wet etching is used to define the InP ridge, (e) the wire is formed through the final evapor ation appropriately 
masking the substrate of the sample (top: side view of the evaporation arrangement, bottom: final result) [l70j . 



associated with QPS. Though rather scarce experimental 
data make quantitative conclusions on the QPS mech- 
anism difficult, the experiments on long Al nanowires 
[sol [l72| have a clear message: No correlations between 



the total normal state wire resistance Rn (compared to 
the quantum resistance unit Rq ~ 6.45 hQ) and super- 
conductivity in such wires was found. 



A convenient material to study the phenomena asso- 
ciated with phase slips is aluminum. Its bulk critical 
temperature T^""^ ^ 1.2 K is relatively low (see Table 
I), hence, the QPS rate should be comparatively high 
enabling pronounced manifestation of fluctuation effects. 
An additional useful feature of aluminum is its peculiar 
size dependence of Tc- Although the origin of this ef- 
fect remains unclear, an increase of Tc with reduction 
of the characteristic dimension of aluminum structures 
(wire diameter or film thickness) is a well-k nown exper- 
imental fact and can be taken as granted This 
effect does not allow to interpret broadening of R{T) de- 
pendencies at temperatures T < Tp""^ in terms of sample 
inhomogeneities, such as constrictions. 

Aluminium was chosen for investigations of ID super- 
conductivity in Refs. [i^, H^. It was demonstrated 
that low energy Ar~^ ion sputtering can progressively and 
non-destructively reduce dime nsions of various nanos- 
tructures including nanowires [l73l . Il74| . The penetra- 
tion depth of Ar~^ ions into Al matrix at acceleration 
voltages of ~ 500 eV is about 2 nm and is compara- 
ble to the thickness of naturally formed oxide. The ac- 
curacy of the effective diameter determination from the 
normal state resistance by SEM and SPM measurements 
is about ± 2 nm. Only those samples which showed no 
obvious geometrical imperfections were used for further 
experiments. To a large extend the method allows one to 
study the evolution of the size phenomenon, eliminating 
artifacts related to uniqueness of samples fabricated in 
independent processing runs. The ion beam treatment 
polishes the surface of the sampl^ removing inevitable 
roughness just after fabrication p^4i (Fig. [^5)1 . If there 



were no detectable geometrical imperfections in the orig- 
inal (thick) wires, they could not be introduced in the 
course of diameter reduction by low energy ion sputter- 
ing. 

After a sequence of sputterings (alternated with R{T) 
measurements) the wire diameter was reduced from -y/s ~ 
100 nm down to ^/s < 10 nm. Experiments were per- 
formed on several sets of aluminum nanowires with length 
X equal to 1, 5 and 10 /xm. For larger diameters > 20 
nm the shape of the R{T) dependence is rather "sharp" 
and can be qualitatively described by the TAPS mecha- 
nism. Note that the abovementioned size dependent vari- 
ation of Tc in aluminum nanowires results in broadening 
of the R{T) transition and significantly reduces applica- 
bifity of the TAPS model [93,] (Fig. 01). When the wire 
diameter is further reduced, deviations from the TAPS 
behavior become obvious (Fig. [26|) . Fits to the TAPS 
model fail to provide any reasonable quantitative agree- 
ment with experiment for diameter values below s^^^ < 
20 nm even if one hypothetically assumes the existence 
of unrealistically narrow constrictions not observed by 
SPM. And, as we already discussed, broadening of the 
R{T) dependencies in aluminum nanowires at T < T^""^' 
can hardly be ascribed to geometrical imperfections, such 
as constrictions. We conclude that the most natural in- 
terpretation of our observations is associated with quan- 
tum fiuctuations. 

Let us now perform a detailed comparison of the data 
with theoretical predictions [H, [U HO] discussed in 
Chapter 5. To begin with, we should select a proper 
formula for the current-voltage characteristics of the wire 
which can be applied to the samples (46j . One can employ 
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FIG. 24; i?(r) dependencies for Al nanowires. At tempera- 
tures <1 K, samples s2 100 ptm long, i?jv =86 kS7 (A) and si: 
10 /.im long, iijv = 8.3 kf2 (t) became superconducting. The 
insets show SEM images and schematic of a typical device 
containing a 10 /xm long, nominally 8 nm wide Al nanowire 
(similar to samples si and s2): the Al layer does not entirely 
cover the Au/Ti pads which were measured in series with the 
superconducting wire (in t he main panel this series resistance 
has been subtracted) [172| ]. 



the general formula (|144p derived for sufficiently long 
wires. It remains not completely clear if this assumption 
is well satisfied for the wires used in experiments [i^ . 
though it appears that at least the longest wires with 
AT = 10 /im are already in this regime [Xji^ > 100 in this 
case). In addition, one can recall that in the experimen- 
tal temperature range not far from Tc there can still be 
sufficiently many quasiparticles above the gap which can 
produce additional dissipation. Although it would not be 
fully justified to use eq. (|147p in this case, still dissipation 
can result in additional inter-QPS interaction which can 
slightly modify the power-law dependencies (|145p . (|146p . 

Fortunately, this effect, even if exists, does not change 
our fitting procedure. In order to describe the QPS con- 
tribution to the wire resistance we will use the formula 



i?(T) = b: 



^°(^)^^°-^cxp(-2.. 



(154) 



while for the non-linear (current-dependent) wire resis- 
tance we will use 



R{I)(xI''. 



(155) 



Here 62 is an unimportant constant which remains the 

same for all samples. Let us again remind the reader that 

the dependence R oa comes from pairs of QPS events. 

As we already discussed, the fugacity p36p depends on 

1/2 

temperature via Score cx; Aq (T). An additional power- 
law dependence enters because of inter-QPS interaction. 
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FIG. 25: SPM images showing evolution of the shape of the 
same aluminium nanowire after several sessions of ion beam 
sputtering. Bright color above the horizonatal plane (initial 
level of the substrate) corresponds to metal, dark color below 
indicates sputtered Si substrate. Note the reduction of the 
initial surface roughness of the nanowire 46]. 
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FIG. 26: R{T) curves for the thinnest samples obtained 
by progressive diameter reduction for the same aluminium 
nanowire with length X = 10 /xm. The TAPS model fitting 
is shown with dashed lines for 11 and 15 nm samples with 
the best fit mean free path Z = 3 and 10 nm, correspondingly, 
Tc = 1.46 K and critical magnetic field 5^(0) = 10 mT. Fits 
to eq. (|151|l are shown by solid lines. For 11, 12, 13 and 15 
nm wires the fit parameters are: A ~ 0.1; Tc'. 1.64 K, 1.52 K, 
1.47 K, and 1.47 K; mean free path I : 7.5 nm, 8.2 nm, 9.5 
nm and 9.5 nm; normal state resistance Rm '■ 7200 kf2, 5300 
kn, 4200 kft and 2700 k^ 
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The parameter k, can be equal to k = 2/x — 3 (as in eqs. 
(|145p . (|146p ) or take somewhat different values due to 
dissipative contribution. Below we will just use it as a fit 
parameter instead of the parameter 

We first determine this parameter by fitting the I — V 
curves at a given temperature to the dependence (|155p . 
see Fig. [27l From these fits we obtain the values for 
parameter k (and, correspondingly, for fi) for each sam- 
ple. Now we perform the fits of the resistance data R{T) 
taken in the hnear regime I < Iq = kBTc/^o ^ Ic 
to the dependence (|154p . There are four parameters 
to specify for our fits: critical temperature Tq, nor- 
mal state resistance Rn, electron elastic mean free path 
I (used to re-calculate the dirty limit coherence length 
^ — 0.85(^00""^^^ )i the numerical factor A of order 
unity. The critical temperature and the normal state re- 
sistance can be trivially deduced from the experimental 
data for R{T). The electron mean free path I can be 
roughly estimated from the normal state resistivity pN 
as the product p^l — 5 x 10~^^ilm^ is a well-tabulated 
value for dirty aluminum. Since the cross section for 
ultra-narrow nanowires is known from SPM measure- 
ments within ± 2 nm accuracy, there remains a certain 
room to choose particular values of the mean free path. 
As a rule of thumb, for all Al nanowires with effective 
diameter ^/s < 20 nm the best-fitted mean free path (at 
low temperatures) was found to be roughly equal to one- 
half of the diameter (Fig. [211 caption). This estimate 
appears quite reasonable taking into consideration that 
at these scales and temperatures electron scattering oc- 
curs merely at the sample boundaries. 

Our experimental data and the fits to the theoretical 
dependence are presented in Fig. [511 We observe that 
as the sample diameter decreases the R{T) curves be- 
come progressively broader, exactly as it is predicted by 
the QPS theory. Our fits demonstrate a good agreement 
between theory and experiment, thus confirming the im- 
portant role of QPS effects in our thinnest wires. For all 
aluminum nanowires the best-fit value for the parameter 
A was determined to be A ~ 0.1. This value is smaller 
than that extracted from the fits to the experimental data 
by Lau et al. [i^ . We believe that this difference can be 
attributed to different geometry and degree of inhomo- 
geneity of samples used in these two experiments. 

We also point out that the temperature dependence of 
all fitted curves R{T) is merely determined by that of the 

1/2 

QPS action Score oc Ag (T) which enters the exponent 
in the expression for the QPS fugacity (I136P whereas an 
additional power-law dependence oc T** turns out to be 
insignificant: The fits remain essentially unchanged if we 
set K ~* 0. This is quite natural since 5corc 3> n and, 
on top of that, the temperature dependence of Scorc{T) 
dominates over that emerging from QPS interaction in 
the temperature interval under consideration. 

To complete our discussion of QPS experiments we 
point out that very recently an additional experimental 
evidence for quantu m p hase slip effects in short MoGe 
wires was reported [175|. In this experiment the super- 
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FIG. 27: V{J) dependencies of the same samples as in Fig. 
26 taken at close temperatures stabilized with accuracy ±0.1 
mK. Solid lines correspond to the dependence i? ~ J" with 
K ~ 0.58, 1.31 and 1.49 (from top to bottom) 46]. 



conducting wires were biased by a large current (just 
slightly below the critical one) in order to decrease an 
effective potential barrier for QPS. Similarly to macro- 
scopic quantum tunneling experiments (see, e.g., Ref. 
[53j and further references therein) the contribution of 
the QPS mechanism was identified in Ref. jl75l | by ob- 
serving low temperature broadening of the distribution 
for switching currents at which the nanowire is driven to 
the normal state by quantum fluctuations of the order 
parameter. Further quantitative investigations of this ef- 
fect as well as the analysis of the role of possibly existing 
external noise and some other extrinsic factors would be 
highly desirable. 

Summarizing our analysis of the experimental data on 
superconducting nanowires, we can conclude that homo- 
geneity of these wires is one of the central issues which 
can influence the data interpretation. The methods en- 
abli ng experiment s on the same nanowire while increas- 
ing [l03l Il45l . Il46l | or decreasing 0, [46l| its characteristic 
dimension to a large extent eliminate artifacts related to 
the uniqueness of samples fabricated in independent pro- 
cessing runs. On the other hand, series of experiments 
on M oGe wires [H, H^] provide large statistics of in- 
dependent samples. Another conclusion is that longer 
nanowires appear to be more suitable for experimental 
investigations of QPS effects while the behavior of very 
short wires could be to a large extent determined by the 
Josephson physics (classical or quantum) which can es- 
sentially mask the QPS one. 

The whole scope of the analyzed experimental data 
obtained in different groups within different fabrication 
methods employing different materials clearly demon- 
strates that at sufficiently low temperatures quantum 
phase slips provide an important mechanism which may 
cause non-zero resistance of superconducting nanowires. 
A vast majority of the data on wires with diameters in 
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the 10 nm range confirms that the TAPS scenario cannot 
explain the broad R{T) dependencies observed in these 
structures and favors QPS [2l,[2^,[3^ as the only possible 
scenario. Although in general this scenario can include a 
number of rather sophisticated physical effects, such as 
quantum phase transitions (see Chapter 5), it appears 
that in many cases the rapid change from superconduct- 
ing to normal behavior can simply be interpreted as a 
crossover between the regimes of low and high QPS rates 
[30| controlled by the wire cross section s or, equivalently, 
by the the wire normal state resistance per unit length 
Rn/X. Also the dominating contribution to the tem- 
perature dependence observed, e.g., in sub-15-nm MoGe 
[3^ I43I I and aluminium 4t | nanowires is well described 



by that of the QPS rate 
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1291, |30|, see Figs. [20] and [ 



Although more experiments would be highly desirable, 
already at this stage we can conclude that the existing ex- 
perimental data unambiguously confirm our understand- 
ing of basic features of QPS physics in superconducting 
nanowires. 



D. Open questions and related topics 

1. Negative magneotresistance 

Application of a weak perpendicular magnetic field 
to thin superconducting strips governed by fluctuations 
revealed an unusual effect: negative magnetoresistance 
(nMR) (Figs. [lS]and[2ni). The phenomenon has been ob- 
served in lead [145,] and aluminium [46] nanowires. The 
effect is seen only in very narrow quasi-lD superconduct- 
ing channels in the resistive state sufficiently far from the 
critical temperature Tc- A trivial explanation related to 
the Kondo mechanism is problematic, as the presence of 
magnetic impurities in those experiments is believed to 
be negligible and - in any case ~ has not been indepen- 
dently proven. Even if we were to assume the presence 
of Kondo impurities, several features of nMR would still 
remain unclear. First, the corresponding magnetic fields 
are too small to polarize magnetic moments of any im- 
purities at such values of T. In addition, it is known that 
aluminum is immune to creation of localize d ma gnetic 
moments with concentration up to few at % |l76j . Sec- 
ond, there is a pronounced diameter dependence, making 
nMR observable only in thinnest wires. Third, the on- 
set of superconductivity is not affected by weak magnetic 
fields: nMR is observed only at the bottom part of the 
R{T) transition. Perhaps we should also mention another 
effect - enhancement of the critical current by magnetic 
field - which was observed in MoGe and Nb nanowires 
[l77j | and attributed to interplay between spin-exchange 
scattering from residual magn etic impurities and the or- 
bital and Zeeman effects |178| . However, it remains un- 
clear if this effect is related to nMR in any way. On top of 
that, presently there exists no evidence for the presence 
of magnetic i mpu rities in experiments with aluminium 
[3 and lead ^U^ . 
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FIG. 28: (a) Resistive transitions in zero magnetic field for 
58 nm wide lead wires of difi'erent thicknesses. Superconduc- 
tivity is suppressed with reduction of the wire thickness. The 
symbols represent the points at which the magnetoresistance 
was measured, (b)-(e) Magnetoresistance for the wire at dif- 
ferent t hick nesses. The numbers indicate the temperature in 
Kelvin 



NMR has been observed in tin microbridges [l79| | and 
aluminium nanostructures dem onstr ating the resistive 
transition anomaly (see below) [l28j |. However, the di- 
mensions of those structures were significa ntly larger as 
compared to ultrathin nanowires [46l Il45| . Presumably 
the origin of nMR in these earlier reports can be ascribed 
to ch arge imbalance jl79l | or to sample inhomogeneity 
[l80j |. NMR h as b een predicted in disordered supercon- 
ducting wires |l8l| . However, the contribution responsi- 
ble for nMR within this model is exponentially small as 
compared to the effective resistance produced by T APS. 
This is clearly not the case in experiments [i^ Il45l | . 

A plausible explanation of the nMR effect could be re- 
lated to the reduction of the energy gap Aq in the mag- 
netic field. It was argued [2^ that this reduction leads 
to the following trade-off. On one hand, the core action 
•Score fx Aq^^ decreases with increasing magnetic field 
and, hence, the QPS rate becomes bigger. On the other 



hand, the quasiparticle resistance i?„ 



exp(Ao/T) 
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FIG. 29: Slowly (~1 h) recorded dependence R(T,Be^t = 
const) for the 11-nm sample from Fig. 1261 The perpen- 
dicular magnetic field Bext =19.6 mT was switched on and 
off several times while sweeping the temperature. The top 
branch corresponds to zero field, while the lower one to the 
field "on". Inset: resistance versus perpendicular magnetic 
field R{Bext,T = const) for the same sample measured at 
constant temperature and small ac current [46j. 



decreases too implying increasing dissipation which, in 
turn, suppresses the QPS rate. Provided the latter trend 
dominates over the former one, the contribution of quan- 
tum fluctuations to the wire resistance gets reduced with 
increasing magnetic field, thus leading to nMR [2^. Yet 
another idea employs possible formation of a charg e im - 
balance region accompanying each phase slip event jl83j | . 
This non-equilibrium region, if it exists, would provide 
dissipation outside the core of a phase slip. The cor- 
responding Ohmic contribution can be effectively sup- 
pressed by the magnetic field, resulting in nMR. How- 
ever, so far the validity of the charge imbalance concept 
was only demonstrated at temperatures sufficiently close 
to Tc and its applicability to QPS is by no means ob- 
vious. This mechanism still requires a solid theoretical 
justification. 



2. Step-like R(T) dependencies in crystalline ID structures 

Rather unusual step-like shape of R{T) transition has 
been reported in tin whiskers 184]. In relatively long 
structures X ^ 0.8 mm with "manually" fixed e lectr odes 
using conducting epoxy and/or Wood metal [l85l | the 
steps were observed at relatively high measuring cur- 
rent densities j ~ 10^ A/cm^ (Fig- EO])- In hybrid 
whisker-ba sed s tructures with lithographically-fabricated 
electrodes |ll6l | (Fig. [6]) the steps on i?(T) dependencies 
were observed between closely located probes {X < 10 
Hm.) and at current densities as low as j ~ 10^ A/crn^ 
(Fig. I3ip . For comparison, the experimental critical 
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FIG. 30: Step-like VJT) dependence in tin whisker at various 
excitation currents Il85ll. 



current density in pure tin whiskers at sufficiently small 
temperatures reaches its theoretical value jc(0) ^ 10^ 
A/cm^. The origin of the phenomenon is not clear. The 
step-like R{T) dependencies were not reported in similar 
experiments using tin whiskers [H, U^. Note that the 
steps on R{T) transition should not be confused with 
step-like I — V characteristics observed in ID supercon- 
ducting wires very close to the critical temperature and at 
bias currents exceeding the temperature-dependent crit- 
ical value / > Ic{T). The origin of the I — V steps is 
usually associated with current-induced non-equilibrium 
effects [mil. 



3. Resistive transition anomaly 

Intensive studies of various lift-off fabricated quasi-lD 
nanostructures revealed an unusual "resistive transition 
anomaly" : Resistance increase above the normal state 
value Rn at the top of the superconducting transition 
(Fig . [321) i29, 130, 131, 180, .186, 187, 188, 189, 190, wj, 
Il92j |. No special correlation has been observed between 
the absolute value of the sample length X (distance be- 
tween the probes) and the magnitude of the anomaly. 
Pronounced "bumps" were observed only in relatively 
short structures with width w > X/10. No dependence 
upon the film thickness within the range 30-90 nm has 
been noticed either. The same statement holds for the 
sample topology: both single-connected and non-single- 
connected structures displayed the bumps. For the same 
multiterminal structure, the anomaly could be clearly ob- 
served for a pair of voltage contacts, while the neighbor- 
ing segment showed no signs of the effect. For a given 
set of contacts, the magnitude of the resistive anomaly 
sometimes depends on the cooling history. Heating up to 
^--^50 K could eliminate the effect. Variation of the mea- 
suring current or application of the magnetic field can 
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FIG. 31; An example of a pronounced step-like V{T) depen- 
dence in a hybrid whisker-based structure similar to that of 
Fig. [6l Arrows indicate the direction of the temperature 
sweep. Small hysteresis is a consequence of the temperature 
measurement finite response time. The inset shows an en- 
larged view of several R{T) steps measured at two different 
currents. The curves are shifted slightly due to z ero offset 
drift of the front-end battery powered preamplifier [116| ]. 



modify the magnitude of the bump (Fig. 

Let us point out that experimentalists frequently use 
lock-in technique in order to increase the signal-to-noise 
ratio. Should a dc component occasionally occur (e.g., 
due to improper grounding), experimental dependencies 
could become proportional to dV/dI{T) rather than to 
the desired R{T) = Vac/lac resulting in "veird" -looking 
R(T) curves. Utilization of rf filters, such as widely 
used commercial Tr-filters, can also contribute to non- 
monotonous R{T) transitions due to re-arrangcmcnt of 
the measuring ac current through the ground loop even 
at rather low frequencies '--^lO Hz. 

Various physical reasons for the resistive transition 
anomaly in superconducting nanostructures have been 
proposed, such as the fluctuation-governed resistive state 
[l93l Il94l . Il95| | nonequilibrium quasiparticlecharge im- 
balance within the locus of N /S boundary [l87l llSSt or 
phase-slip centers [l29l Il95| [. It should be noted that 
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FIG. 32: Normalized resistance as a function of tempera- 
ture for five aluminium samples A-E showing the resistance 
anomaly. Longer wires display a smaller and broader peak . 
Inset: The sample configuration for three shortest wires [1291 ]. 



the resistive transition anomaly has only been reported 
in relatively inhomogeneous lift-off fabricated nanostruc- 
tures with the width of the R{T) transition strongly ex- 
ceeding that predicted by the fluctuation models. It has 
been shown that the phenomenon can be qualitatively 
explained by simple geometric considerations concerning 
the shape of the N/S phase bo unda ry of an inhomoge- 
neous wire in t he re sistive state [l96j | . The approach has 
been extended jl80| demonstrating that the combination 
of the two effects provides reasonable agreement with ex- 
periments: (i) a geometry effect of the formation of non- 
perpendicular (to the wire axis) N/S boundaries and the 
corresponding rearrangement of current across the wire, 
and (ii) the existence of nonzero and strongly anisotropic 
effective resistance of the nonequilibrium superconduct- 
ing region close to the N/S interface. Considered sepa- 
rately, neither of these two contributions could provide 
quantitative agreement with the experimental data. 

To summarize, it is possible that the resistive transi- 
tion anomaly originates from "dirty physics" related to 
inhomogeneity of the finite length lift-off fabricated wires. 
This is supported by the fact that no si gns of th e ano maly 
were reported in pure ID systems [l8l.ll9l. IIIGL I185| . 



4- Tc dependence on wire diameter. 

Since early days of low-dimensional superconductivity 
it is known that the critical temperature of thin films 
T^^ and wires T^^ differ from the corresponding bulk 
value r^^"""'. In indium, aluminum and zink nanowires 
Tc increase s wit h decreasing the characteristic dimension 
20, 45, 50. Il02j |. In lead, niobium and MoGe an oppo- 
site tendency has been observed (s^, H ^M, flo3 . ri05| . 
No noticea ble variations of Tn have been reported in tin 
nanowires [iM [ifrTi [iM [T09l | . 

The early models considered size-dependent modifica- 
tion of the electron-phonon spectrum and found the cor- 
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FIG. 33; Top part of the resistive transition R{T) for alu- 
minium multiterminal nanostructure at different measuring 
ac currents (top) and at various exte rnal magnetic fields nor- 
mal to the sample surface (bottom) [l80l ]. 



respo nding var i ation of the critical temperature of thin 
fihus [l97l. llQSl . Il99l |. However, later it was noticed that 
the phonon spectrum renormalization can significantly 
contribute to variations of Tc only in ultra-thin struc- 
tures, while experimentally the effect is observed also in 
rather thick films (e.g. up to t ~ 50 nm in the case of 
aluminum), where the phonon spectrum is almost iden- 
tical to that of bulk samples. Additionally, it was found 
that the critical temperature Tq^ of ultra-thin quench- 
condensed films depends not only on thickness, but also 
on annealing made in situ in the measuring chamber at 



cryogenic temperatures [200|, indicating that the mor- 
phology of the film is also important. Another signifi- 
cant effect is an interplay between disorder and electron- 
electron interactions which can lead to substantial sup- 
pression of Tc in reduced dimensions ^llOl] . On the o ther 
hand, the so-called shape resonance effect l20l| can 
be responsible for the enhancement of Tc- To conclude, 
after years of intensive experimental and theoretical stud- 
ies, the origin of the size dependence of the critical tem- 
perature Tc remains not fully understood and in many 
cases it can be determined by a number of different fac- 
tors. 



VII. PERSISTENT CURRENTS IN 
SUPERCONDUCTING NANORINGS 

A. Persistent currents and quantum phase slips 

It is well known that superconducting rings pierced by 
external magnetic flux develop circulating persistent 
currents (PC) which never vanish. This phenomenon is 
a fundamental consequence of macroscopic phase coher- 
ence of Cooper pair wave functions. In the case of bulk 
metallic rings fluctuations of the phase ip of the order 
parameter can be neglected, i.e. ip can be considered as 
a purely classical variable. In this case the total phase 
difference ip{X) — ip{0) accumulated along the ring cir- 
cumference X = 27ri? is linked to the external flux 
inside the ring by the well known relation 



(^(A) - ^(0) = + 



(156) 



where p is an integer number, (j)^ = ^x/^o and, as be- 
fore, $0 is the superconducting flux quantum. This rela- 
tion implies the existence of a phase gradient along the 
ring which, in turn, means the presence of discrete set of 
current and energy states labelled by the number p. At 
sufficiently low temperatures T <^ Ag quasiparticles are 
practically irrelevant and the grand partition function of 
the ring takes the form 



where 



En 



E 

— oo 



eM~Ep{(^x)/T), 



R 



(157) 



(158) 



defines flux-dependent energy states of a diffusive super- 
conducting ring with radius R and cross-section s. The 
ground state energy E{(j)x) =iniiip Ep((f>x) is a periodic 
function of the flux with the period $0. The deriva- 
tive of the ground state energy with respect to the flux 
defines the persistent current 



(159) 
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(a) £(0)A 




(b) 






















FIG. 34: The energy (a) and persistent current (b) for a su- 
perconducting nanoring as functions of magnetic flux without 
fluctuations, with weak and with strong QPS effects shown 
respectively by solid, dashed and dotted lines [33l ]. 



Combining eqs. (|158p and (|159p one finds 

/ = 2TremVsNoDAos, (160) 
where m is the electron mass and 

vA(h^) = minri \v-\ ^ 1 (161) 

' 2mR $0/ 

is the superconducting velocity. Both the current / and 
the velocity are periodic functions of the magnetic flux 
showing the familiar sawtooth behavior, see also Fig. 

m 

This picture remains applicable as long as fluctuations 
of the superconducting phase (f can be neglected, i.e. as 
long as the ring remains sufficiently thick. Upon decreas- 
ing the thickness of a superconducting wire ^/s down to 
values ~ 10 nm range one eventually reaches the new 
regime in which quantum fluctuations of the phase gain 
importance and, as it will be demonstrated below, essen- 
tially modify the low temperature behavior of supercon- 
ducting nanorings. 

As we have already learned, at low T the most impor- 
tant fluctuations in superconducting nanowires are quan- 
tum phase slips. Each phase slip event implies transfer of 
one flux quantum $o through the wire out of or into the 
ring and, hence, yields the change of the total flux inside 
the ring by exactly the same amount. In other words, 
each QPS yields a jump between two neighboring energy 
states Ep (|158p (i.e. between two neighboring parabolas 
in Fig. 34). As a result of such jumps the flux inside the 
ring fluctuates, its average value ($($3;)) decreases and 
so does the persistent current /. One can also anticipate 
that the magnitude of this effect should increase with the 
ring perimeter X. This is because the QPS rate jqps 



increases linearly with X, i.e. the bigger the ring perime- 
ter the higher the probability for QPS to occur anywhere 
along the ring. 

A quantitative theory of QPS effects in nanorings was 
proposed by Matveev et al. .33J . These authors employed 
a model of a closed chain of Josephson junctions in which 
case QPS can occur only across such junctions. Here will 
extend our theory of QPS (Chapter 5) and analyze the 
effect of quantum fluctuations on persistent currents in 
uniform superconducting nanorings. The key conclusions 
remain the same for both models. 

Before turning to technical details it is instructive to 
point out a formal equivalence between the phenomenon 
discussed here and charging effects in ultrasmall Joseph- 
son junctions (or Cooper pair boxes) [2l|. Indeed, the 
energy states of our ring (|158p driven by the (normal- 
ized) external magnetic flux (j)^ (see Fig. |34| are fully 
analogous to such states of a capacitor in a Cooper pair 
box driven by the (normalized) gate charge Qx/'ie. Fur- 
thermore, we will demonstrate that there exists a direct 
analogy between QPS events changing the flux inside the 
ring by $0 and tunneling events of single Cooper pairs 
changing the capacitor charge by 2e. This equivalence is 
reminiscent of the well known duality between phase and 
charge representations [2l| and will be exploited in our 
consideration. 

We start from the case of sufficiently thick supercon- 
ducting rings where QPS effects can be neglected. In this 
case by virtue of Poisson's resummation theorem (see, 
e.g.. Sec. 3.3.2 in Ref. |2l'|) one can identically trans- 
form the ring partition function (|157p to the following 
expression 

Z^^= exp(i27rfc0,) dOo / VO exp{- So[0]) , 

(162) 

where 

and 1 /En plays the role of a " mass" for a " particle" with 
the coordinate 9. This coordinate represents an effective 
angle. It is formally analogous to the Josephson phase 
variable in which case En just coincides with the charging 
energy. 

Let us now consider thinner rings where QPS effects 
gain importance. In this case we should include all tran- 
sitions between different energy states (|158p labelled by 
the number p. As we already discussed these transitions 
are just QPS events with the rate Jqps defined in eq. 
(|130p . Let us fix p = and take into account virtual 
transitions to all other energy states and back. Perform- 
ing summation over all such contributions we arrive at 
the series in powers of the rate 7qps (or fugacity y) sim- 
ilar to eq. (|135p . Assuming that both the ring perimeter 
X and its thickness are not too large (the latter condition 
restricts the parameter fi) one can neglect weak logarith- 
mic interaction (jl32p between different phase slips. Then 
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spatial QPS coordinates Xn become exact zero modes and 
can be trivially integrated out. As a result, we arrive at 
the following contribution to the partition function 



n=0 



2n 



dri. 



dT2n 



: J2 '5(</.(T)-0„(T))exp (- / 



where, as before, the summation is carried out over neu- 
tral charge configurations ()137p . 



3 = i 



(165) 



and 0(t) is the theta-function. We note that the term in 
the exponent describes virtual energy changes occurring 
due to QPS events. What remains is to add up similar 
contributions from all other parabolas with p ^ 0. In 
this way we arrive at the final expression for the grand 
partition function of the nanoring 



(166) 



Together with eq. (|164p this result fully accounts for all 
QPS events in our system. 

In the absence of QPS (i.e. for "fqps 0) cqs. ()166p . 
()164p obviously reduce to (I157|) . For non-zero 7qps one 
can rewrite the partition function p64p . (I166P in the 
equivalent form of the following path integral 



Z^^ = exp(i27rfc(/)^) / ddo / P61 exp(-5[6l]), 

7 Oq 

(167) 



k— — oo 



where 



s[e] = So[e]~-fQPs / dT cos 0{t) 



(168) 



In order to demonstrate that the partition function (|167p , 
(fT68| is identical to that defined in eqs. (fT64| . (fT66l) it 
suffices to perform the Hubbard-Stratonovich transfor- 
mation of the kinetic term (|163p (which amounts to in- 
troducing additional path integral over ^(t)), to formally 
expand exp(—S[9]) in powers of jqps and then to inte- 
grate out the 0- variable in all terms of these series. After 
these straightforward steps we arrive back at eqs. (|164p . 
(|166p . We also note that this procedure is described in 
Sec. 3.3.5 of Ref. '21], therefore we can avoid further 
details here. 

Eqs. (|167p . (|168p define the grand partition func- 
tion for a quantum particle in a cosine periodic poten- 
tial which is, in turn, equivalent to that for a Josephson 
junction in the presence of charging effects [2lJ. This 



partition function can be identically rewritten as (cf. eq. 
(3.94) in Ref. [21j]) 

OO 

Z^^ = J2 ^M-EM-)IT), (169) 

p— — oo 

where Ep{(j)x) are the energy bands of the problem de- 
fined by the solutions of the well-known Mathieu equa- 
tion. For 7gps <C Ep one has 



-—r arcsm 
27r^ 



T (^) ) ''"^"'^^^ 



(170) 

i.e. the energy bands remain nearly parabolic Ep{(t)) ^ 
Ep((f>) except in the vicinity of the crossing points where 
gaps open due to level repulsion (see Fig. [M]) . In this 
limit the value of the gap between the lowest and the first 
excited energy bands just coincides with the QPS rate 
SEqi — jQps- For larger jqps the bandwidth shrinks 
while the gaps become bigger. In the limit jqps ^ Ep 
the lowest band coincides with 



Eo{<l>) = -^7Q''ps(^^fl)'^'e-«V^^^^^(l - cos(2^0,)). 

(171) 

The gap between the two lowest bands is SEqi = 
\/iqpsEr. 

These results are sufficient to evaluate PC in super- 
conducting nanorings in the presence of quantum phase 
slips. As before, taking the derivative of the ground state 
energy with respect to the flux and making use of eq. 
(|130p we find that for 7qps <C Ep and outside immedi- 
ate vicinity of the points (j)x = 1/2 PC is again defined 
by eqs. (|160p . (|16ip . In the opposite limit Jqps ^ Ep 
we find 

/ = /osin(27r0j,), 

/o = iE|7S''p5(^^i^)'/'e-«VW^. (172) 



We observe that in the latter limit PC is exponentially 
suppressed as 

io oc exp(-i?/i?e), (173) 

where neglecting a numerical prefactor we estimate 

i?,-eexp(S'core/2). (174) 

This simple formula sets the size scale beyond which one 
would expect PC to be exponentially suppressed in su- 
perconducting nanorings with ^/s < 10 nm. E.g. for ^ of 
order 100 nm and S'corc ~ 10 eq. (|174p yields Rc of order 
of a micron. 

We would like to emphasize that during our analysis we 
employed only one approximation neglecting logarithmic 
inter-QPS interaction effects. If needed, such effects can 
also be included into our consideration and may only lead 
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to unimportant modifications in our results for rings with 
very large perimeters X. 

The above results demonstrate practically the same 
qualitative features as those found within a different ap- 
proach for the model of nanorings formed by Josephson 
chains [s^. In particular, both for granular and for ho- 
mogeneous rings the dependence of PC on gradually 
changes from sawtooth to sinusoidal as the ring perimeter 
X increases. This change is accompanied by suppression 
of PC amplitude which eventually becomes an exponen- 
tially decaying function of the ring radius R (|173p in the 
limit R > Re- In other words, eq. (|173p sets the length 
scale Rc beyond which persistent currents in supercon- 
ducting nanorings should be essentially washed out by 
quantum fluctuations. In some sense the scale 2'kRc plays 
the role of a dephasing length in our problem demonstrat- 
ing that zero-point fluctuations can destroy macroscopic 
phase coherence down to T — s- even in superconducting 
systems. 

This conclusion is qualitatively consistent wi th t h e re- 
sults derived for normal metallic conductors |202l . l203l 
|204, 205| and single channel rings coupled to dissipative 
baths [2O5I l20d |207| . In fact, a close similarity between 
the pa rtition fu nction (|167p . (|168p and those studied in 
Refs. |206L l207l | exists on a formal level as well: These 
problems are described by similar path integrals except 
in the case of normal rings coupled to dissipative envi- 
ronments one should include the corresponding non-local 
in time term into the action instead of the last term 
in eq. (|168p . At sufficiently large X these path inte- 
grals can also be handled in exactly the same manner: 
In both problems the dominating contribution comes 
from instantons (kinks) describing quantum tunneling of 
the angle variable 9 between different topological sec- 
tors of the problem. In all cases this procedure yields 
low temperature PC in the form (|172|) where one finds 
Iq cx exp(— (X/L y)^) in the case of Caldeira-Leggett en- 
vironments [205l . l206l . \20T\ and Iq cx exp(-X/L^) in the 
case of so-called dirty electron gas environments j207']. 

An important qualitative difference between our 
present problem and those of normal nanorings with dis- 
sipation lies in the fact that in the latter case dissipation 
explicitly violates time-reversal symmetry (thus causing 
genuine decoherence of electrons) , while no such symme- 
try is violated in eqs. (I167II. (11681) . Hence, in our case 
- in contrast to [205l . l206l l207t - exponential suppres- 
sion of PC (|172p can also be interpreted as a non-trivial 
coordinate-dependent renormalization effect. This dif- 
ference is just the same as that between dissipativeless 
Cooper pair tunneling (which, just like QPS, opens inter- 
band gaps in the spectrum of Josephson junctions) and 
dissipative single electron tunneling (which opens no such 
gaps). For more details on the latter subject we refer the 
reader to the review [2l|. 

Finally we should mention that one can also consider 
a slightly modified situation of rings consisting of thicker 
and thinner parts, as it is shown in Fig. 1351 Assum- 
ing that QPS effects are negligible in a thicker part of 




FIG. 35: Superconducting ring with a "quantum phase slip 
junction" . 



the ring and they can only occur in its thinner part (of 
length d) we arrive at exactly the same results as above 
in which one should only replace the ring perimeter X 
by the length d. In particular, it follows from our analy- 
sis that such rings also exhibit the property of an exact 
duality to mesoscopic Josephson junctions if we identify 
Eji with the junction charging energy Eq and the QPS 
rate Jqps with the Jose phso n energy Ej. For this rea- 
son Mooij and Nazarov [208j | suggested to call systems 
depicted in Fig. [35] "quantum phase slip junctions" and 
argued that any known result on electron transport in 
circuits containing Josephson junctions can be exactly 
mapped onto a dual result for a QPS junction in a dual 
circuit. This observation can be used in metrology, e.g., 
for practical implementation of the electric current stan- 
dard in the above structures. Mooij and Harmans |209l | 
proposed to use rings with QPS junctions for experimen- 
tal realization of quantum phase slip fiux qubits. 

It is also worth pointing out that previously a sim- 
ilar exact duality between phase and c harge var iables 
in Josephson junctions was discussed [2ll . l210l l21l| . see, 
e.g., eq. (5.76) in Ref. [2l|. According to this property 
the Josephson junction itself can also be a QPS junction 
if we interchange the canonically conjugate phase and 
(quasi)charge variables (f irq/e. In this sense one can 
identify the angle variable 6 with the (quasi)charge Trg/e 
passing through the QPS junction. 

B. Parity effect and persistent currents 

In our previous analysis of QPS effects in supercon- 
ducting nanorings we followed the standard technique de- 
veloped for grand canonical ensembles, i.e. we implicitly 
assumed that the total number of electrons in the system 
N may fluctuate and the chemical potential fl is fixed. 
Obviously, this assumption is not correct for rings which 
are disconnected from any external circuit. In that case 
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electrons cannot enter or leave the ring and, hence, the 
number TV is strictly fixed, but the chemical potential fl, 
on the contrary, fluctuates. It turns out that novel effects 
emerge in this physical situation. These effects will be 
discussed below in this section. To a l arge extent we wil l 
follow the analysis developed in Refs. (sH, ll53L I212L l213l | . 

It is well known that thermodynamic properties of iso- 
lated superconducting systems a re sensitive to the p arity 
of the total number of electrons MB, HH, [Hi even 
though this number is macroscopically large. This par- 
ity effect is a fundamental property of a superconduct- 
ing ground state described by the condensate of Cooper 
pairs. The number of electrons in the condensate is neces- 
sarily even, hence, for odd N at least one electron remains 
unpaired having an extra energy equal to the supercon- 
ducting energy gap Aq. This effect makes thermody- 
namic properties of the ground state with even and odd 
N different. Clear evidence for such parity effect was 
demons trated experim entally in small superconducting 
islands pTsl . [218L l219t . 

At the first sight, this parity effect should have little 
impact on the supercurrent because of the fundamen- 
tal uncertainty relation SN5(p > 1. Should the electron 
number N be fixed, fluctuations of the superconducting 
phase (fi become large disrupting the supercurrent. On 
the other hand, in transport experiments with fluctua- 
tions of If being suppressed the parity effect cannot be 
observed because of large fluctuations of N. 

Consider now isolated superconducting rings pierced 
by the magnetic flux ■ In accordance with the number- 
phase uncertainty relation the global superconducting 
phase of the ring fluctuates strongly in this case, how- 
ever these fluctuations are decoupled from the supercur- 
rent and therefore cannot influence the latter. In this 
situation the parity effect may substantially modify PC 
in superconducting nanorings at sufficiently low temper- 
atures. In particular, we will show that changing the 
electron parity number from even to odd results in spon- 
taneous supercurrent in the ground state of such rings 
without any externally applied magnetic flux. In other 
words, our fundamental conclusion will be that the BCS 
ground state of a canonical ensemble with odd number of 
electrons is the state with spontaneous supercurrent. 



1. Parity projection formalism 

In order to systematically investigate interplay be- 
tween the parity effect and persistent currents in super- 
conducting nanorin gs we will empl o y th e parity projec- 
tion formalism [H, [HI, [HO, [22l which we will 
briefly outline here. 



is linked to the canonical one Z{T, N) as 

Z(r,A)= £^(T,iV)exp('^ 



N=0 



(176) 



Here and below TL is the system Hamiltonian. Inverting 
this relation and defining the canonical partition func- 
tions Zf. and Zo respectively for even (iV = N^) and odd 
{N = No) ensembles, one gets 



1 

2^ 



due 



Ze/oiT) 

where 

Ze/oiT, ~^) - iTr{[l±(-l)^] 



Z,/oiT,iTu), (177) 



= -{Z{T,f^)±Z{Trti + i^T)) (178) 

are the parity projected grand canonical partition func- 
tions. For 3> 1 it is sufficient to evaluate the integral 
in (|177p within the saddle point approximation 



(179) 



where ^e/o = ^T\uZf,/o{T,fi) are the parity projected 
thermodynamic potentials. 



^e/o = ri/ - Tin 



-[l± e-z^lf^b-n/) 



The grand canonical partition function 



and VLf/i, = -Tin [Tr{(±l)^e-^(^-'^^)}]. "Chemical 
potentials" /ie/o are defined by the saddle point condition 

Ne/o = -d^e/o{T,fle/o)/dfLe/o- 

The main advantage of the above formalism is that it 
allows to express the canonical partition functions and 
thermodynamical potentials in terms of the parity pro- 
jected grand canonical ones thereby enormously simpli- 
fying the whole calculation. We further note that 0/ is 
just the standard grand canonical thermodynamic poten- 
tial and f2b represents the corresponding potential linked 
to th e partition function Z(T, p, + zttT). It is easy to see 
[22 1| that in order to recover this function one can eval- 
uate the true grand canonical partition function Z(T^ fi), 
express the result as a sum over the Fermi Matsubara 
frequencies usf — 2ttT{1 + 1/2) and then substitute the 
Bose Matsubara frequencies ujb — 2ttTI {I — 0, ±1,...) 
instead of w/. This procedure automatically yields the 
correct expression for Z(T, fi -f i-nT) and, hence, for S7;,. 

Having found the thermodynamic potentials for the 
even and odd ensembles one can easily determine the 
equilibrium current /g/o- Consider, as before, isolated 
superconducting rings pierced by the magnetic flux $2,. 
Making use of the above expressions one finds PC circu- 
lating inside the ring: 



Z(r,^) = Tre-'3(W-AiV) 



(175) 



4/o - If 



± 



,J3{Q.i-Q.;) _|_ 



(180) 
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where the upper/lower sign corresponds to the even/odd 
ensemble and we have defined 



/o 



dQe/o 



f/b 



(181) 

Eqs. (|180p . (|18ip represent a direct generaUzation of the 
grand canonical formula p59p to canonical ensembles. 



2. Homogeneous superconducting rings 



N 





FIG. 36: Superconducting ring with embedded SNS junction 
of length d i34. 



Let us first consider homogeneous nanorings with cross 
section s and perimeter X = 2'kR. As before, rings will 
be assumed sufficiently thin, ^/s < X^. On the other 
hand, below we will neglect QPS effects, i.e. describe su- 
perconducting properties of such rings within the parity 
projected mean field BCS theory. As we have already 
learned, this description is justified provided the condi- 
tion ^ 1 is satisfied. Hence, the ring should not be 
too thin and the total number of conducting channels 
should remain large Af ^ 1. In addition, the perimeter 
X should not exceed the scale 2ttRc (|173p . Finally, we 
will neglect the difference between the mean field values 
of the B CS order parameter for the even and odd en- 
sembles |22Cll . l22ll | . This is legitimate provided the ring 
volume is large enough, V = Xs ^ I/NqAq, where, as 
before, Nq is the density of states at the Fermi level. 

Evaluating thermodynamic potentials ^f/b and ex- 
pressing the result in terms of the exci tatio n energies Sk 
and the order parameter Aq one finds 221 



Qf =n- 2T^ In (2 cosh |^ 



fifc = fi-2r^ln('2sinh 



2TJ ' 



(182) 
(183) 



Here we introduced the electron density in the case of 
even/odd total number of electrons n^/o = N^/o/V. 

We observe that in the case of odd ensembles there 
exists an additional term which modifies the flux depen- 
dence of PC and, as we will demonstrate below, leads to 
a number of fundamentally important effects. Unfortu- 
nately, in the case of homogeneous rings the difference 
between PC in even and odd ensembles turns out to be 
hardly observable. Indeed, these difference is inversely 
proportional to the total number of conducting channels, 
(/e — Io)/Ie ^ 1/M. For this reason the parity effect 
remains vanishingly small in generic metallic rings with 
^ 10"^ • Oil the other hand, for ultra-thin nanorings 
with A/" < 10 PC is essentially wiped out due to prolif- 
eration of quantum phase slips. Estimating ^ A/'^/^, 
we conclude that for '-^ I (i.e. when the QPS fugacity 
is already large and, hence, quantum suppression of PC 
becomes very strong) the number of conducting channels 
yet remains parametrically large A/" ~ ^/^ 3> 1. This esti- 
mate demonstrates that the parity effect on PC is never 
important in the case of homogeneous superconducting 
nanorings and appears to be practically unobservable in 
such systems. 



where Q. = | Ao| VA -f Trj^}, 

i-l^{-^^l~-Air))\uir)^~,^ (184) 

is the single-particle energy operator, — pVg + 

Y^'^p -|- A§, where p is the quasiparticle momentum, = 

{p^ — ji{<^x)/2ra and the superconducting velocity Vs is 
defined in eq. (|16ip . 

The above equations allow to fully determine PC in 
superconducting nanorings with even and odd number 
of electrons /g/o. The parity effect beco mes observ- 
able at sufficiently low temperatures |215l | T < T* k, 
Ao ln(A'oVVAoT*). Here we consider the most interest- 
ing limit T < hvp/X. From eqs. P^ - (fT551) we find 
that in this limit the current le exactly coincides with 
the standard grand canonical result / = engVsS with 
n.s = Up , whi l e the current for the odd ensemble reads 
[3ll2T3ll223l.[22l 



eUoVsS ■ 



e^sgn?;^ 



3. SNS rings 

We now turn to the situations in which parity effect 
gains importance and can be directly probed in modern 
experiments. Let us slightly modify our system and con- 
sider a relatively thick superconducting ring with large 
(7^ S> 1 interrupted by a thin wire of length d with only 
few conducting channels 



(186) 



(185) 



thus forming a weak link inside the superconducting ring, 
see Fig. [351 Without loss of generality this wire can be 
considered normal no matter if it is made of a normal 
or a superconducting material. In the latter case quan- 
tum ffuctuations would fully suppress the order parame- 
ter inside such a wire bringing it into the normal state. 
In contrast, QPS effects in superconducting parts of the 
ring can be neglected thus making the mean field BCS 
description applicable. A clear advantage of these struc- 
tures in comparison to homogeneous rings is that in the 



former case the effect of the electron parity number on 
PC can be large due to the condition (|186p . 

In order to evaluate PC in such structures we again 
employ the parity projection formalism. According to 
eqs. (|180p . (|18ip we need to evaluate both the currents 
///f, and the difference between the "Fermi" and "Bose" 
thermodynamic potentials fib — ilf = ^bf ■ The currents 
If jb can be conveniently expressed via the phase differ- 
ence across the weak link Lp ~ 27r(/)j; by means of the 
general formula [22 5i] 
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^ 2e smy 
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where the sum runs over conducting channels of the nor- 
mal wire and the function Wi{Lo) was evaluated in Ref. 
[225j . The difference of thermodynamic potentials f2b/ 
is defined as a sum of the contributions from supercon- 



ducting iS^b]^ normal iS^b]') parts of the ring, 
former is evaluated with t he aid of eqs. (|182p . (|183p which 
yield the standard result [215j | 
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FIG. 37; The ratio between canonical and grand canonical 
values of PC le/o/Ig.c. (represented by the term in the square 
brackets in Eq. (|190|l versus in a single mode QPC at differ- 
ent temperatures for even (three upper curves) and odd (three 
lower curves) ensembles. Here we have chosen the channel 
transmission T = 0.99 Q. 
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while the latter is found by integrating If /biv) (|187p over 
the phase tp. 

We now consider several important limiting cases. The 
first limit is that of a very short normal wire d —^ 0. This 
is essentially the limit of a quantum point contact. In 
practice this limit is realized for d <ti ^. Even smaller 
values of d are required provided the contact transmission 
is^mall. In the limit of a quantum point contact one finds 



(189) 
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In 



where % define transmissions of the conducting channels. 
Substituting this function int o th e above equations one 
arrives at the final result [s^, |213| | 



FIG. 38: Andreev levels inside a quantum point contact and 
their occupation at T = for even (a) and odd (b) ensembles 
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Here £^{(p) = Ay 1 - 7^ sin^(.^/2) are Andreev energy 
levels in a quantum point contact. 

The first line of eg . (I190p defines the standard grand 
canonical result 12261 . |227] while the term in the square 
brackets accounts for the parity effect in our system. For 
Mn = 1 and at T = this term reduces to unity for even 
ensembles and to zero for odd ones, i.e. PC turns out to 



be totally blocked in the case of odd number of electrons 



De totailyoio 

P^IIiilli 



The physics of this blocking effect is rather transparent 
and can be understood as follows. We first recall that in 
the limit c? — > the Josephson current can be expressed 
only via the contributi ons from discrete Andreev energy 
states E± = ±e{'p) as [229l . l230t 
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Using the Fermi filling factors for these states f± {E± ) — 
[1 + exp(±e(p) / T)]~^ one arrives at the standard grand 
canonical results |226l . l227j | . In the case of superconduct- 
ing rings with fixed number of electrons A'^ these filling 
factors should be modified. Let us set T ^ 0. For even 
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N all electrons are paired occupying available states with 
energies below the Fermi level (see Fig. [55k ) . In this case 
one has f-{E^) = 1, f+{E+) = 0, the current is entirely 
determined by the contribution of the quasiparticle state 
E- and eq. (|19ip yields the same result as one for the 
grand canonical ensemble. By contrast, in the case of odd 
N one electron always remains unpaired and occupies the 
lowest available energy state - in our case E+ - above the 
Fermi level. Hence, for odd N one has f±{E±) = 1 (Fig. 
l38b). the contributions of the two Andreev states in eq. 
()19ip exactly cancel each other, and the current remains 
zero for any tp or the magnetic flux <i>x. This is just the ef- 
fect of parity-induced blocking of PC derived above from 
formal considerations. 

For T > Eq. (fTM]) demonstrates that both for 
even and especially for odd N the current-phase rela- 
tion for QPC may substantially devi ate fr om that de- 
rived for the grand canonical ensemble [226j |. see Fig. [37l 
For even ensembles the supercurrent increases above its 
grand canonical value. This effect is mainly pronounced 
for phases ip not very far from = tt and - at suffi- 
ciently low T - it becomes progressively more important 
with increasing temperature. On the contrary, for odd 
ensembles the supercurrent is always suppressed below 
its grand canonical value. This suppression is gradually 
lifted with increasing temperature, though at phases ip 
in the vicinity of the point <p = tt blocking of PC may 
persist up to sufRciently high T. Eq. (|190p also shows 
that in quantum point contacts with several conducting 
channels and at T — > the current through the most 
transparent channel will be blocked by the odd electron. 
Hence, though blocking of PC remains incomplete in this 
case, it may nevertheless be important also for quantum 
point contacts with Mn > 1- 

Let us now turn to another important limit of su- 
perconducting rings containing a normal wire of length 
d > ^0 ^ hvp / ^a- In contrast to the case d ^ con- 
sidered above, the Josephson current in SNS structures 
cannot anymore be attributed only to the discrete An- 
dreev states inside a weak link, and an additional con- 
tribution from the states in the continuum should also 
be taken into account. Furthermore, for any non-zero d 
there are always more than two discrete Andreev levels 
in the system. Accordingly, significant modifications in 
the physical picture of the parity effect in such SNS rings 
can be expected. 

The key difference can be understood already by com- 
paring the typical structure of discrete Andreev levels in 
SNS junctions (Fig. I39p with that of a quantum point 
contact (Fig. [55)1 . As before, in the limit T all states 
below (above) the Fermi level are occupied (empty) pro- 
vided the total number of electrons in the system is even 
(Fig. [39k). If, on the other hand, this number is odd 
the lowest Andreev state above the Fermi energy is oc- 
cupied as well (Fig. [39b) thus providing an additional 
contribution to the Josephson current. This contribu- 
tion, however, cancels only that of a symmetric Andreev 
level below the Fermi energy, while the contributions of 
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FIG. 39: Andreev levels in a single mode SNS junction with 
d = 6hvF I Ao and their occupation at T = for even (a) and 
odd (b) ensembles [3^ . 
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FIG. 40: The zero temperature current-phase dependence 
([T93)) for SNS rings with = 1: I^{}p) (dashed) line and 
Io(}P^ (solid) line [34| . 



all other occupied Andreev levels and of the continuum 
states remain uncompensated. Hence, unlike in the case 
d — > 0, in SNS rings one should not anymore expect the 
effect of PC blocking by the odd electron but rather some 
other non-trivial features of the parity effect. 

This conclusio n is fully confirmed by our quantita- 
tive analysis [H, |213| |. Let us restrict our attention to 
transparent SNS junctions in which case the function 
Wi{(M) = WJuj) is the same for all transmission channels 
and reads |225l | 

iy(.)^(^^ + ljcosh^— j 

Substituting this function into (|187p and repeating the 
who le ca lculation as above, we arrive at the final result 
[13, l213l | which takes a particularly simple form in the 
limit T -> and d > ^qi 

evpAfn J evpUn ( 7rsgn</?\ 
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FIG. 41: The same as in Fig. 1401 only for the odd ensembles 
(second Eq. (IT93J) and for 7V„ = 2, 3 and 4 [H]. 
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FIG. 42: The zero temperature current-phase relation lo(}p) 
(— TT < ^ < tt) for A/'n = 1 and different values of the param- 
eter y = dl^/hvF [34| . 



This result applies for — tt < < tt and should be 27r- 
periodically continued otherwise. We observe that at 
T = the current /g again coincides with that for the 
grand canonical ensembles 231, 232], while in the case of 
odd ensembles the current-phase relation is shifted by the 
value 7r/-/V. This shift has a simple interpretation being 
related to the odd electron contribution {2e/h)dEi/dip 
from the lowest (above the Fermi level) Andreev state 
Ei{(f) inside the SNS junction. As we have expected, this 
contribution indeed does not compensate for the current 
from other electron states. Rather it provides a possi- 
bility for a parity-induced 7r-junction state ^233j in our 
system: According to Eq. (|193p for single mode SNS 
junctions the "sawtooth" current-phase relation will be 
shifted exactly by tt, see Fig. |40l More generally, we 
can talk about a novel 7r/-/V-j unction state, because in 
the odd case the minimum Josephson energy (zero cur- 
rent) state is reached a.t (p = ±Tr/N, see Fig. |40l For 
any J\fn > 1 this is a twofold degenerate state within the 
interval — tt < (p < n. In the particular case Afn = 2 the 
current-phase relation loi'p) turns 7r-periodic, see Fig. 
41. 

Let us recall that the 7r-junction state can be realized in 
SNS structures by driving the electron distribution func- 
tion in the contact area out of equilibrium |234l . |235| . |236| | . 
Here, in contrast, the situation of a tt- or Tr/A'^-j unction is 
achieved in thermodynamic equilibrium. Along with this 
important difference, there also exists a certain physical 
similarity betw een the effects discussed here and in Refs. 
l234l |235|. l236ll: In both cases the electron distribution 



function in the weak link deviates substantially from the 
Fermi function. It is this deviation which is responsible 
for the appearance of the 7r-junction state in both phys- 
ical situations. 

Perhaps the most spectacular physical consequence of 
the parity effect in SNS rings is the presence of sponta- 
neous supercurrents in the ground state of such rings with 
odd number of el ectro ns. Similarly to the case of stan- 
dard TT-junctions [233j | such spontaneous supercurrents 
should flow even in the absence of a n externally applied 
magnetic flux. Unlike in Ref. [233j . however, here the 
spontaneous current state occurs for any inductance of 



the ring because of the non-sinusoidal dependence loif)- 
Consider, for instance, the limit d ^ Hvp/Aq. In the 
case of odd number of electrons the ground state energy 
of an SNS ring can be written in a simple form 



$2 TrhvpJ^n 
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(194) 



where the first term is the magnetic energy of the ring (£ 
is the ring inductance) while the second term represents 
the Josephson energy of an SNS junction. Minimizing 
(|194p with respect to the flux $ one immediately con- 
cludes that the ground state of the ring is a twofold de- 
generate state with a non-vanishing spontaneous current 



/ = ± 



1 + 



2ewi?A/'„ C 
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flowing either clockwise or counterclockwise. In the limit 
of small inductances C ^ this current does not vanish 
and its amplitude just reduces to that of the odd electron 
current at (/? — > 0. One finds [34j, 213]: 



Isp = eAld/h^vp, 
hp = evp/i^d, 



d > hvp/Ao. 



(196) 
(197) 



At d hvp/ Aq the amplitude of the current hp can be 
evaluated numerically, see Fig. [43l One observes that - 
in agreement with eq. (|196p - Isp increases linearly with d 
at small d, reaches its maximum value /max ~ OAeAo/h 
at d ~ ^ and then decreases with further increase of d 
approaching the dependence (|197|) in the limit of large 
d. For generic BCS superconductors the magnitude of 
this maximum current can be estimated as /max '^-^ 10 -f- 
100 nA. These values might be considered as surprisingly 
large ones having in mind that this current is associated 
with only one Andreev electron state. 

Note that in the above analysis we merely assumed 
that the normal wire is sufficiently clean and, on top of 
that, is in a good electric contact with superconductors. 
Since both these assumptions can be violated in a realis- 
tic experiment it is important to discuss the correspond- 
ing modifications of our results. 
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FIG. 43: The spontaneous current amplitude Isp as a function 
of the parameter y at T = 0. In the inset, the same function 
is shown on the log — log scale. Dashed lines indicate the 
asymptotic behavior of Iap{y) in the limits of small and large 



Assume first that the transmissions of both NS inter- 
faces are small Ti^2 ^ 1- In this limit electron trans- 
port across the junction is mainly due to resonant tun- 
neling through discrete energy levels inside the normal 
metal. For simplicity we will restrict our analysis to a 
single channel junction J\fn ~ 1. The most interesting 
physical situation is realized in the limit of short junc- 
tion d <^ hvp/Ao- In the case of a one dimensional 
metal of length d, the level spacing in the vicinity of the 
Fermi energy is Se ~ Hvp/d. Hence, the condition for the 
short junction regime can also be represented in the form 
i5e 3> Aq. Electron tunneling causes a non-zero linewidth 
of the energy levels which is proportional to Ti^2^e- This 
value is much smaller than Se, hence, the resonances re- 
main well separated. In this situation it suffices to take 
into account only the closest level to the Fermi energy 
inside the normal metal. 

As before, making use of eqs. (|187|1 combined wit h the 
proper expression for the M^-function [225] one finds [213l | 
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A^T sin ip 
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eicp) = A^l-rsin2(^/2), T = Fi + F2, Ti^^/h = 

'?l,2W_F/2(i are the tunneling rates, T„^ax ~ 4FiF2/F^ and 
the total transmission probability at the Fermi energy T 
is given by the Breit-Wigner formula 
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FIG. 44: Andreev levels in a single mode SNS junction with 
d = hvp/A: (a) Ti 2 = 1 and (b) Ti,2 <C 1 and A/T = 0.5. In 
both cases T = 1 EM- 



where en is the energy of a resonant level. 

It follows from Eqs. pM)) . that - although 

the transparencies of both barriers are low - the total 
transmission T and, hence, the Josephson current shows 
sharp peaks provided the Fermi energy becomes close to 
a bound state inside the junction. On the other hand, 
eqs. (|198l) . (|199p demonstrate that even in the vicinity 
of resonances the behavior of the Josephson current as 
a function of the phase difference ip and temperature T 
can substantially deviate from that for transparent SNS 
junctions. 

In order to understand the physical reasons for such 
a difference it is instructive to compare the structure of 
discrete Andreev levels for ballistic (71. 2 = 1) SNS junc- 
tions with that for junctions with weakly transmitting 
NS interfaces 71 _2 ^ 1, see Fig. The spectrum of 
the latter system consists only of a single non-degenerate 
state £o(<^) in the interval < eo < Aq (Fig. [jib). As 
a result, the behavior of eo{ip) at small ip is smooth and 
the derivative of Sq with respect to (p has no jump at 
p = 0. In contrast, in the case of ballistic SNS junctions 
discrete levels become split at arbitrary small values of 
the phase p (Fig. Hlk ) and the derivative of the lowest 
Andreev level with respect to p> acquires a jump ai p> — 0. 
As this feature is absent in resonant SNS junctions the 
spontaneous current in the ground state of such systems 
can only develop at not very small ring inductances. The 
results for PC in SNS rings with resonant transmission 
are presented in Fig. [55l They clearly demonstrate that 
at sufficiently low temperatures the "7r-junction" state 
should be realized in the case of odd number of electrons. 

Finally let us turn to the case of a disordered normal 
wire. The difference between PC values for odd (Jo) and 
even (/e) ensembles is rel a.ted t o the min i gap value egjip ) 
inside the normal metal ^M, [HI IHO, Ii5l|, EM, \lm- 
This relation acquires a particularly simple form in the 
limit T ^ 0: 
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FIG. 45: Zero temperature current-phase relations I^/^i'fi) for 
T = 0. 9, Fi = F2 = Fo and different values of the parameter 
Ao/Fo [213^ . 
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FIG. 46: Phase dependence of the Josephson current at T = 
for the odd and even number of electrons in the ring |153i] . 



The last term in this equation describes the contribu- 
tion to the current from the "odd" electron occupying 
the lowest available state above the minigap Sg{ip) in the 
density of states of the normal metal. 

Let us first evaluate PC for the even ensemble /g. As 
before, at T = this current identically coincides with 
one calculated for the grand canonical ensemble. The 
results for the current-phase relation lei^) are displayed 
in Fig. US] at various impurity concentrations. PC in the 
odd ensemble /q at T = can be evaluated from eq. 
(I200p . Combining our results f or Ip jip) with those for the 
minigap Sg{ip) derived in Ref. [153l | we arrive at a typical 
dependence loi^fi) displayed in Fig. 351 We observe that 
at sufficiently large values of (/? < tt the absolute value of 
the odd electron contribution to PC 2edeg/dif exceeds 
the term lei^fi) and the total current Iq changes the sign. 
This non-trivial parity-affected current-phase relation is 
specific for SNS rings with disorder and it substantially 
differs from the current-phase relations derived above for 
SNS rings with ballistic and resonant transmissions. 

At the same time, as in the previous cases, in the odd 
ensemble there exists a possibility both for a 7r-junction 
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FIG. 47: Josephson energy E{lp) of an SNS ring as a function 
of the phase difference if for the even and odd ense mbles. The 
solid curve corresponds to a vr -junction state [l53l ]. 



state and for spontaneous currents in the ground state 
of the system without any externally applied magnetic 
fiux. Let us evaluate the ground state energy of the SNS 
junction by integrating Eq. ()200|1 with respect to the 
phase if. One finds 



(201) 

where E^/g{ip) are the ground state energies of SNS 
junction for even and odd number of electrons in the 
ring. While the energy Ee{(p) is always non- negative and 
reaches its minimum at 1^9 = 0, in the odd case the ground 
state energy Eo{(p) can become negative reaching its ab- 
solute minimum at ip — tt. This physical situation of a 
TT-junction is illustrated in Fig. [171 

It is easy to find out under which conditions the tt- 
junction state becomes possible. For that purpose it is 
sufficient to observe that for any impurity concentration 
£'e(7r) = a/c/e, where Ic is the grand canonical critical 
current at T = and the prefactor a ~ 1 depends on 
the form of the current-phase relation. The 7r-junction 
condition Eoiir) < is equivalent to the inequality 



£g(0) > alc/e. 



(202) 



It is obvious from Fig. Hz] that in the many channel limit 
the inequality (|202p cannot be satisfied for sufficiently 
large I. On the other hand, for sufficiently short mean 
free paths Ic oc P decays faster with decreasing I as com- 
pared to the minigap Sg(0) oc I, and the 7r-junction state 
beco mes p ossible. In particular, in the diffusive limit one 
finds [HI Ic ^ 10.82eTh/ei?Ar = l.bSeAfnVpP /d^ and 
a ~ 1.05, where i?Ar is the Drude resistance of a normal 
metal. 

Combi ning thes e results with the expression for the 
minigap [l49l Il52j £g(0) ~ 3.12eTh, from the condition 
(|202p we observe that in the odd case the 7r-junction state 
is realized provided the number of conducting channels 
in the junction Mn is smaller than 



7V„ < 0.65d/L 



(203) 
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This condition is not very restrictive and it can certainly 
be achieved in various experiments. For sufficiently dirty 
junctions it allows for a formation of the 7r-junction state 
even in the many channel limit. The condition (|203p can 
also be rewritten as qn < 1.73, where — 8Afnl/Sd is 
dimensionless conductance of a normal wire. 

The condition for the presence of spontaneous currents 
in the ground state of SNS rings with an odd number of 
electrons is established analogously, one should only take 
into account an additional energy of the magnetic field 
produced by PC circulating inside the ring. The ground 
state with spontaneous currents is possible provided the 
total energy of the ring -Etot(7'') becomes negative, i.e. 

Etotin) = 1.8£Th [9N - 1.73] + ^^j^ < 0. (204) 

This condition is more stringent than that for the tt- 
junction state, but can also be satisfied provided C ex- 
ceeds a certain threshold value which can roughly be es- 
timated as ~ 0.1<i>g/eTh- 

We conclude that in the diffusive limit the current- 
phase relation in the odd case is entirely different from 
that in the ballistic case. Also the restriction on the num- 
ber of conducting channels Afn in the normal metal (|203p 
is less stringent. This feature of diffusive SNS rings is 
rather advantageous for possible experimental observa- 
tion of the effects discussed here. 

Still, in practice it would be necessary to fabricate SNS 
rings with few conducting channels in the normal wire 
Afn ^10. This condition can hardly be met for conven- 
tional normal metals. It appears, therefore, that most 
promising candidates for practical realization of such 
structures are junctions with N-layers formed by carbon 
nanotubes or organic molecules. In this respect it is im- 
portant to point out that observations of clear signatures 
of dc Josephson effect in superconducting junctions with 
a weak link formed by carbo n nanotubes were reported by 
several experimental groups [237l . [238l [239l [240l . l24ll . [24a | 
Also intrinsic superconductivity in carbon was claimed 
[243l I244I l245l | , though at this stage it appears that more 
experimental support for these claims would be desirable. 
Regardless of this latter issue, SNS junctions with carbon 
nanotubes are most likely objects in which it would be 
possible to observe the influence of the parity effect on 
persistent currents in superconducting nanorings. 

VIII. SUMMARY 

It is well established that fluctuations play an impor- 
tant role in structures with reduced dimensionality. In 
superconducting materials clear signatures of fluctuation 
effects exist already above the critical temperature Tc 
Slim. Also Tc itself can be reduced due to fluctu- 
ation effects. This reduction is particularly important in 
disordered 2D films and quasi-lD wires [llOj |. In this pa- 
per we addressed fluctuation effects which occur in ultra- 
thin superconducting wires at temperatures below the 



mean field BCS critical temperature. Superconducting 
properties of such systems have been intensively studied 
- both theoretically and experimentally - during past 
years. The key conclusions of these investigations can be 
summarized as follows. 

Thicker superconducting wires are characterized by 
very small Ginzburg numbers Giio ^1 and diameters 
typically > 100 nm. In such systems the superconduct- 
ing transition is broadened due to thermally activated 
phase shps (TAPS) [M IH which cause non-zero resis- 
tance R{T) at temperatures close enough to the critical 
temperature Tc — T ^ Tc- Upon decreasing tempera- 
ture TAPS events become less likely and quantum fluc- 
tuations of the order parameter take over. This is the 
regime of quantum phase slips (QPS) which sets in at 
T < Ao(T). As long as the wire is sufficiently thick and 
the Ginzburg number remains very small, Gim <?SCl, 
QPS events are rare and typically do not lead to any 
measurable consequences. Hence, the behavior of thicker 
wires remains essentially superconducting outside an im- 
mediate vicinity of the critical temperature. Upon re- 
duction of the wire diameter below ^50 nm the QPS 
rate increases drastically. In this regime the wire resis- 
tance R{T) still decreases with temperature but may re- 
main well in the measurable range down to very low T. 
Provided the wire diameter is decreased further the di- 
mensionless conductance of a wire segment of length ^ 
eventually becomes of order ~10 or smaller. In such 
wires QPS proliferate causing a sharp crossover from a 
superconducting to a normal behavior. For generic pa- 
rameters this crossover is expected for wire diameters in 
the 10 nm range [2^, [2^, [s^l . This crossover was indeed 
observed in a number of experiments in wires with thick- 
nesses exactly in this range. Thus, intrinsic supercon- 
ductivity in wires with diameters ^10 nm is destroyed 
by quantum fluctuations of the order parameter at any 
temperature down to T — 0. 

Theoretical analysis of nanowires reveals further in- 
teresting effects, like QPS-binding-unbinding quantum 
phase transition (QPT) [11] between superconducting 
and non-superconducting phases which is predicted to 
occur as the impedance of a superconducting wire be- 
comes of order of the quantum resistance unit Rq « 6.5 
kfi. For typical parameters this condition is also achieved 
for wire diameters in the 10 nm range. To the best of 
our knowledge, no clear experimental evidence for this 
phase transition exists so far. This can be due to rather 
stringent requirements: QPT can only be observed in 
long nanowires at sufficiently low temperatures, ideally 
at r ^ 0. Interesting effects may occur also in short 
nanowires forming weak links between superconducting 
electrodes. In many respects such systems can behave 
similarly to Josephson nanojunctions and weak links. 

Novel effects are also expected in superconducting 
nanorings. While rings formed by thicker wires demon- 
strate the standard behavior familiar from the bulk sam- 
ples, the situation changes drastically as soon as the wire 
thickness gets reduced down to w 10 nm or below this 
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value. In this case QPS effects become important leading 
to strong fluctuations of the magnetic flux inside the ring. 
As a result, the amplitude of persistent current decreases 
and its flux dependence changes from the sawtooth-like 
to a smoother one j33j | . For such rings even at T = per- 
sistent current gets exponentially suppressed by quantum 
fluctuations provided the ring radius R exceeds the crit- 
ical scale Rc ^ i exp(S'corc/2) where (Score ~ 3? G^iil^ 
is the action of the QPS core. For typical wire parameters 
the length R^ can be of order of a micron or even smaller. 
This length constitutes another fundamental scale asso- 
ciated with quantum phase slips. 

Yet another important factor which may influence per- 
sistent currents in isolated superconducting nanorings is 
the electron parity number. This influence is particu- 
larly strong in nanorings containing a weak link with few 
conducting modes. Changing the electron parity number 
from even to odd may result in spontaneous supercurrent 
in the ground state of such rings without any externally 
applied magnetic flux [s^]. At T = this current is 
produced by the only unpaired electron which occupies 
the lowest available Andreev state. Under certain con- 
ditions this spontaneous supercurrent can reach remark- 
ably large values up to ^ eAo/h ^ 10 100 nA which 
can be reliably detected in modern experiments. 



equations: 
G 



F 



G 



1 / (r^T(Xi)^|(X2)> ~{^J\{X2)^^^{x,)) \ 

I \ (^t(Xi)7a|(X2)) (r-V3T(^i)^|(^2)) ; ' 

1 / (rV'T(^i)^i(^2)) -(^i(X2)^T(^i)> ^ 

i \ (^t(Xi)v-i(X2)> (r-iv3T(^i)^i(^2)> ; ' 

1 / (r^j(Xi)^|(X2)) -(^|(X2)^j(Xi)) \ 

I \ {^{x,)^{x^)) (r-v3I(Xi)^j(X2)) ) ' 

1 / (rv'|(Xi)^^(X2)) -(^KX2)^I(Xi)) \ 

I \ (^|(Xi)V'i(X2)> (r-iv3I(^i)^i(^2)> j ■ 



Here t/i^ , -0^ are the electron annihilation and creation 
operators, the symbol T (T^^) stands for the time (anti- 
time) ordering. 

Let us introduce the complete basis of the eigenfunc- 
tions x„(r) of the single electron Hamiltonian iJp: 



(Al) 



Here are the electron eneries in the normal state. 
Defining the quasiparticle energies in the superconduct- 
ing state E„ — S,^ + Aq as well as the BCS coherence 
factors 



i(l 



En 
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iG{t, ri,r2) = ^ Xn(ri)xn(r"2) ulAn{t) + vlBn{t) 
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iG{t,ri,r2) = \vlA„{t) + ulBnit) 

n 

iF{t,ri,r2) = iF{t,ri,r2) = ^ XTi(r"i)x„(r2)M„w„ 



where 



An{t) 



Bn{t) 



-iE, 



t [ 0{t)-fn 
[ l-fn 



fn 
-t) - fn 



fn 



fn - 6{t) 



(A3) 



(A4) 



APPENDIX A: EQUILIBRIUM 
GREEN-KELDYSH FUNCTIONS 



We formally define the sub-blocks G, F, G, F in the 
Green-Keldysh matrix (fTB|) by means of the following 



Here we introduced the occupation probabilities of the 
quasiparticle states /„. In thermodynamic equilibrium 
these filling factors take the universal form 



fn 



1 



=s„/r ■ 



(A5) 
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APPENDIX B: WARD IDENTITIES 

One can demonstrate that the matrix Green-Keldysh 
functions satisfy the following Ward identities 



iF0 + ipF = Ga, (-(3+—{V,f3}]F 



+ Fa, |^/3+_{V,/3}jG 
+ 2iAoGaJd - 2iAoFajF, 
-iF$ - 0F = Fa, (-$ + ^{V, p}\ G 



+ 2iAoFajF - 2iAodajG, (Bl) 



Gp-I3G = iGa, ( -P + —{V, /3} ) G 



iF^A(3+—{V,p}]F, 
G$- $G = iPa, ( -/§ + ^{V, /3} ) F 

- iFa,[(3+^{V,P}]F, (B2) 



where 



pF{t,r) 



/3s(t,r) ; 
is an arbitrary diagonal 2x2 matrix. 

APPENDIX C: KERNELS 

Explicit expressions for the functions xj I'ead 



Xa(0) 



2iVn N< 



X 



+ 



Dq^ 



(B3) 



(a-6)'+£>v 



66-A2\ {E,+E2){l-h-h) 



+ 1- 



E1E2 J {oj + iOf - {El + E2Y 
ii^2-Al\ {E^-EM- f2) 



E1E2 J {w + iQf - {El - £-2)2 J ' 

(CI) 



Dq^ 



^0 



+ 



6)2 + r'V^ii^2 

(£1 + ^2) (1^/1-/2) 

(w + i0)2 - {El + £^2)2 

{Ei-E2){h-h) 



{uj + i0)2 - {El - £2)2 



(C2) 



Xl{Q) 



2m'^NoD 



(^1-6)' 



A2 
^0 



(£l+£2)(l-/l-/2) 



+ 



(w + iO)2 - {El + £2)2 

(£i-£;2)(/i-/2) 

(W + «0)2 - (£1 - £2)2 



(C3) 



Xd{Q) 



e^NoD 



d^id^2 



(6-6)2 



{El~E2){h- f2) 

E2 J {u> + i0)2 - {El - E2Y 
iii2 + Al\ {El + E2){1 - fl - f2) 



£1^2 ; (w + io)2 - (£1 + £2)2 



(C4) 



Here we defined £1,2 = yJ(,i,2 + ^% /i,2 = 1/(1 + 
exp[£i,2/T]). 

In a number of limiting cases the kernels Xj can be 
evaluated exactly. First let us set T = and consider the 
limit of small wave vectors Dq^ <C Aq. This regime 
is relevant, e.g., in the context of microwave absorption. 
At \uj\ < 2Ao we find 



XA 

Xj 

XL 

Xd 



■ arctan 



UJ 



8e2jVoA2 
2m2cr_DAo 



: arctan 



v/4A^ 



'-K 



2ar.Ao 
a;2 



K 



2Ao 
V2A0 



E 



2Ar 



while at |a;| > 2Ao one finds 



XA 



Xj 



XL = 



+ i 



8e2jVoA2 

2m?aDAQ 



2An 



4A2 



(C5) 



TT 
'2 



In 



2An 



4A2 



TT 
'2 



e 

4A, 







-K 



2Ao^ / 2Ao 
UJ - 2Ao 



a; + 2Ao V'^ + 2Ao 



Xd 



1R 

to 



K 



E 



+ t- 



.ao 



4A, 







-is: 



/2A( 



w - 2Ao 



a;-|-2Ao V'^ + 2Ao 



ui J \a;-|-2Ao 

oj - 2Ao 



(C6) 



Here and below K{x) and E{x) arc complete elliptic in- 
tegrals of the first and second kind respectively. Utilizing 
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the relation p6)) one can verify that the expressions for 
XL and xd agree with the well known results [s^, [36j . 

At temperatures close to Tc one can derive analytic 
expressions for the kernels in the limit Z?^^, Aq <C 
2itTc. At < 2Ao these expressions read 



7C(3) N^Al No \Lu\Dq^ 



2n^ T2 2T 2Ao 



V2Ao 
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X n 
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4A2 



(C7) 
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Here n(x, y) is the complete elliptic integral of the third 
kind. At higher frequencies |a;| > 2Ao we find 
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APPENDIX D: RELATIONS BETWEEN THE PHASE AND THE ELECTROMAGNETIC POTENTIALS 



Integrating out the electromagnetic fields in eq. ([7T|) we arrive at the action ((72|) where the function J^{uj,q) is 
defined by the following general expression 



( XJ , ,2 I XL „2\ ( C 



XD 



I XjXl 



(Dl) 



The electromagnetic potentials are linked to the fluctuating phase of the order parameter field via the following saddle 
point conditions 
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